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Chapter  1 
Overview 


1.1  Introduction 

The  research  described  in  this  report  concerns  the  following  fundamental  issues  in  control 
system  design: 

1.  Selection  of  meaningful  yet  analytically  tractable  performance  criteria, 

2.  Sensitivity  of  closed  loop  stability  and  performance  to  both  structured  and  unstruc¬ 
tured  plant  uncertainty  and  model  errors,  and 

3.  Approximation  of  ideal  compensators  with  implement  able  ones. 

Since  the  seminal  work  in  this  area  [52],  there  has  been  very  substantial  progress.  [21] 
contains  a  good  survey  and  a  comprehensive  bibliography  for  the  field  up  to  the  time  of  writ¬ 
ing  in  1983.  Since  then  developments  for  linear  lumped  parameter  systems  have  focussed  on 
computational  techniques  and  applications.  Although  [52]  suggests  applicability  to  systems 
with  irrational  transfer  functions,  actual  development  of  a  theory  taking  into  account  the 
characteristics  of  infinite  dimensional  plants  began  with  the  mutually  independent  work  [9] 
and  [19]. 

During  the  last  five  years  there  has  been  significant  progress  on  results  for  distributed 
parameter  systems.  However,  the  progress  has  not  matched  that  of  the  results  lumped 
parameter  systems  over  a  comparable  period.  This  is  partly  because  the  research  related 
to  H°°  control  of  distributed  parameter  systems  has  generally  appeared  in  a  very  abstract, 
highly  mathematical  form,  in  which  the  connection  to  real  problems  is  made  only  in  a  much- 
used  simple  delay  example. 

In  the  work  we  describe  here,  we  have  focussed  our  theoretical  work  on  issues  which  arise 
in  applying  the  theory  to  non-trivial  examples.  These  examples  can  serve  as  archetypes  of 
distributed  parameter  systems  of  interest  to  the  Air  Force. 

This  research  is  intended  to  develop  both  theory  and  algorithms  capable  of  providing 
realistic  control  systems  for  physical  plants  which  are  appropriately  modeled  as  infinite  di¬ 
mensional  linear  systems.  Typical  physical  plants  would  be  large  orbiting  radar  antennas, 


4 


large  robot  arms  such  ?s  the  shuttle  Remote  Manipulator  System,  and  high  precision  or¬ 
biting  optical  systems  such  as  surveillance  satellites.  We  hypothesize  that  maintaining  an 
infinite  dimensioned  model  as  long  as  possible  during  the  design  process  will  lead  to  designs 
which  better  capture  intrinsic  infinite  dimensional  characteristics,  at  lower  computational 
cost  and  with  greater  analytical  certainty,  than  alternative  procedures  which  employ  finite 
dimensional  approximations  as  a  precursor  to  design.  However,  we  must  ultimately  check 
this  hypothesis  by  comparing  designs. 

In  pursuing  this  program,  our  approach  is  twofold:  First,  we  must  examine  sufficiently 
realistic  models  of  typical  prototype  plants  to  see  what  characteristics  models  for  theoretical 
analysis  should  have  in  order  to  yield  useful  results.  Our  work  to  date  along  this  line  indicates 
that  much  existing  theoretical  work  makes  overly  strong  assumptions,  which  prevent  the  work 
from  being  applicable  to  even  basic  plants  of  interest.  Having  distilled  from  such  prototype 
plants  what  we  consider  to  be  essential  characteristics,  we  must  derive  theoretical  procedures 
for  obtaining  a  solution. 

Secondly,  we  recognize  two  facts  about  infinite  dimensional  models,  and  direct  our  re¬ 
search  to  deal  with  these  realities:  One  fact  is  that  even  infinite  dimensional  models  will 
themselves  only  be  approximations,  and  the  other  is  that  infinite  dimensional  models  which 
have  the  fidelity  we  seek  to  any  complex  physical  system  will  probably  not  be  tractable  ana¬ 
lytically  (in  the  sense  of  our  being  able  to  compute  explicit  expressions  as  we  can  for  simple 
delay  systems).  Thus,  having  an  explicit  inner/outer  factorization  as  we  allow  ourselves  to 
assume  in  order  to  develop  theoretical  results,  is  not  a  realistic  assumption,  even  for  some¬ 
thing  so  relatively  simple  as  a  damped  Timoshenko  beam  model  (compared  to  the  shuttle 
Remote  Manipulator  System). 

1.2  Summary 

Our  main  results  are  as  follows: 

1.  (Chapter  2)  We  have  also  calculated  and  analyzed  transfer  functions  for  certain  beam 
models  [11],  [31],  for  use  as  testbeds  of  our  theoretical  results.  We  now  have  five  transfer 
function  models  of  differing  complexity,  derived  from  physically  motivated  plants,  which 
have  led  us  to  the  different  issues  we  are  studying.  These  are:  the  delay  (with  rational 
factor),  the  damped  Euler- Bernoulli  beam,  the  damped  Timoshenko  beam,  the  two-delay 
multivariable  plant  (with  rational  factors),  and  the  multivariable  damped  Euler-Bernoulli 
beam  with  torsion. 

2.  (Chapter  3)  We  have  extended  outer  factor  absorption  results  to  cover  certain  irrational 
outer  factors  [12].  This  ensures  that  this  critical  step  (which  is  frequently  ignored  in  the  lit¬ 
erature)  has  a  sound  basis.  The  method  of  proof  may  yield  a  useful  computational  technique 
for  implementation,  as  well. 

3.  (Chapter  4)  We  have  solved  the  scalar  mixed  sensitivity  problem  with  irrational  outer 
factors  and  rationed  weighting  functions  [13],  [14].  In  [12]  we  explicitly  calculated  a  standard 
transformation  of  this  criterion,  leading  us  to  the  observation  that  in  the  case  of  plants  with 
poles  in  the  closed  left  half  plane,  for  independently  determined  weighting  functions  the 
infimal  norm  is  independent  of  the  outer  factor  of  the  plant.  On  the  other  hand,  we  also 
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observe  that  the  choice  of  weighting  functions  should  depend  on  the  outer  factor  of  the  plant, 
and  in  fact  this  reasoning  leads  to  a  need  to  extend  results  further  to  the  case  of  irrational 
complementary  sensitivity  weighting. 

4.  (Chapter  5)  We  have  successfully  solved  a  mixed  sensitivity  problem  for  a  non-trivial 
multi-input/multi-output  system  [47].  We  believe  this  is  the  first  such  example  solved.  (In 
[34]  results  are  presented  which  appear  quite  general,  but  in  fact  implicitly  assume  the 
satisfaction  of  a  certain  special  commutativity  property  which  does  not  usually  nold,  as  in 
our  example.) 

5.  (Chapter  6)  We  have  developed  a  technique  for  numerically  computing  the  infimal  norm 
of  the  mixed  sensitivity  criterion  without  an  explicit  inner/outer  factorization  of  the  plant 


1.3  Papers  and  Reports 

The  following  is  a  list  of  publications  written  using  the  support  of  this  contract.  The  results 
in  these  documents  have  been  incorporated  into  the  present  report. 

1.  D.  S.  Flamm,  “A  Model  of  a  Damped  Flexible  Beam,”  ISS  Report  No.  54,  June  14, 

1990,  Department  of  Electrical  Engineering,  Princeton  University. 

2.  D.  S.  Flamm,  “Outer  Factor  ‘Absorption’  for  H°°  Control  Problems,”  ISS  Report 
No.  55,  July  31,  1990,  Department  of  Electrical  Engineering,  Princeton  University. 
(Submitted  to  The  International  Journal  of  Robust  and  Nonlinear  Control.) 

3.  D.  S.  Flamm  and  H.  Yang,  “Some  Comments  on  the  H°°-Optimal  Scalar  Mixed  Sen¬ 
sitivity  Problem,”  ISS  Report  No.  56,  August  16,  1990,  Department  of  Electrical 
Engineering,  Princeton  University. 

4.  D.  S.  Flamm  and  H.  Yang,  “H°°-Optimal  Mixed  Sensitivity  for  General  Distributed 
Plants,”  ISS  Report  No.  57,  August  31,  1990,  Department  of  Electrical  Engineering, 
Princeton  University.  (Submitted  to  IEEE  Trans.  Automatic  Control.) 

5.  D.  S.  Flamm  and  H.  Yang,  “H°°-Optimal  Mixed  Sensitivity  for  General  Distributed 
Plants,”  Proc.  1990  Conference  on  Decision  and  Control,  Dec.  5-7,  1990,  pp.  134-139, 
Honolulu,  HI.  (condensed  version  of  ISS  Reports  56  and  57.) 

6.  D.  S.  Flamm,  H.  Yang,  Q.  Ren  and  K.  Klipec,  “Numerical  Computation  of  Inner 
Factors,”  ISS  Report  No.  58,  August  8,  1990,  Department  of  Electrical  Engineering, 
Princeton  University. 

7.  H.  Yang  and  D.S.  Flamm,  “Mixed  Sensitivity  Design  for  a  Real  Multivariable  Delay 
Problem,”  ISS  Report  No.  60,  Sept.  7,  1990,  Department  of  Electrical  Engineering, 
Princeton  University.  To  appear  in  Proc.  American  Control  Conference,  June  26-28, 

1991,  Boston,  MA. 
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1.4  Notation 

H3,  H°°:  Hardy  spaces  on  the  right  half  plane. 

L1,  L3,  L°°:  Lebesgne  spaces  on  the  real  axis  R  or  imaginary  axis  jR. 

H3:  the  orthogonal  complement  of  H3  in  L3,  so  L3  =  H3  ©  H3. 

L^x,  (H3xl)_,  etc.:  matrices  with  entries  in  the  corresponding  spaces. 

C:  tLe  subspace  of  L°°  consisting  of  bounded  continuous  functions. 

IT+:  the  projection  L3  — >  H3. 

Et_:  the  projection  L3  —*  H3. 

11k:  the  projection  L3  — >K,  where  K  is  a  closed  subspace. 

IIk:  the  projection  from  La  onto  a  subspace  K. 

z*(s):  the  involution  of  x(s),  i.e.,  x*(s)  =  x(— I).  Here  the  bar  denotes  the  complex  conju¬ 
gate. 

i,,z0:  inner  and  outer  factors,  respectively,  of  x  6  H°°. 
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Chapter  2 
Models 


In  this  chapter  we  present  four  models  which  we  have  been  examining  in  order  to  motivate 
and  guide  oui  investigation. 

Since  one  of  the  main  applications  in  which  we  are  interested  is  the  control  of  large  space 
structures,  three  of  the  models  are  of  flexible  beams.  We  think  of  these  as  building  blocks 
from  which  we  well  construct  mere  complex  and  more  realistic  space  structure  models  as 
the  theory  of  design  with  them  matures.  The  simplest  model  is  a  Euler- Bernoulli  beam 
with  viscous  damping.  We  also  present  here  a  damped  Timoshenko  beam  model,  and  beam 
model  with  two  degrees  of  freedom.  The  Euler-Bemoulli  beam  is  a  single-input /single-output 
model  wHch  we  h“.v*  substantial  detail,  and  its  characteristic*  have  influenced 

several  of  our  directions  of  research.  The  Timoshenko  beam  is  also  a  SISO  model,  but  it  is 
substantially  more  difficult  to  analyze,  mainly  because  an  explicit  expression  for  a  transfer 
function  does  not  seem  to  be  available.  We  present  it  here  to  illustrate  the  complexity 
involved,  and  as  an  object  of  future  study.  Similarly,  the  transfer  funcion  for  the  two-degree- 
of-freedom  beam  is  intractable  analytically,  but  we  intend  it  a*  an  iniria!  object 
study  for  multi-input /multi-output  space  structure  models. 

The  fourth  model  is  a  multivariable  delay  model  for  an  automotive  system.  This  is  of 
interest  for  several  reasons:  it  is  a  realistic  MIMO  irrational  transfer  function,  it  was  an 
intractable  problem  in  its  exact  form  (prior  to  the  current  work),  and  it  has  been  examined 
by  means  of  approximation  by  other  researchers,  so  that  computational  results  are  available 
for  comparision  purposes. 


2.1  A  Damped  Euler-Bernoulli  Beam 

The  model  we  present  in  this  section  is  relatively  tractable  by  analytical  means.  We  are 
able  to  find  exact  singularities  and  approximate  zeros,  and  to  characterize  the  asymptotic 
behavior  for  high  frequency. 
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2.1.1  Model  Derivation 


With  the  addition  of  two  tractable  types  of  viscous  damping  to  the  classical  Euler-Bernoulli 
beam  model,  one  obtains  for  the  partial  differential  equation  [Clough  and  Penzien,  p.  302], 


d3 


d3v 


:  EI-=~, -  +  c,I 


dx 2  \  dx 3 


d3  v 

iwdt 


d\ 


dv 


+  mw+crn 


(2.i; 


where  v(x,  t)  is  the  lateral  displacement  of  the  beam,  p(x,  t )  is  the  external  load  on  the  beam, 
c  is  the  resistance  to  transverse  velocity  and  c,  is  the  resistance  to  strain  velocity,  c  has  the 
interpretation  of  external  damping,  and  c ,  has  the  interpretation  of  internal  damping.  We 
are  interested  in  the  case  of  no  external  load,  (p(z,  t)  =  0),  and  we  assume  a  beam  of  length 
1  with  free  boundary  conditions  at  both  ends. 

With  the  intent  of  computing  the  transfer  function  from  <x  torque  applied  at  one  end  to 
the  displacement  at  the  other  end,  we  take  Laplace  transforms,  obtaining 


(El  +  c,Is)vw(x,  s)  +  (ms3  +  cs) v(x,  s)  =  0 
The  problem  is  separable,  and  setting 

v(x,s)  =  ifi(x)  ■  w(s) 

we  obtain 


<p(4)(x) 


<p(x) 

The  solution  is  of  the 

form 

t>(x,  s)  =  (k i  cos  (ax] 

+ 

where  w(s)  and  the  constants  ki 

are  d 

conditions  at  each  end  axe  related  to  1 

/  v(0,s)  ) 

/ 

v'(0,s) 

v"(0,s) 

\t/"(0,  s)/ 

V 

ms *  +  cs 
El  +  c,I  s 


(2.2) 


(2.3) 


2  4' 
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0 
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0 

a 

0 


1  0  \ 
0  a 
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\ 
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kt 

k3 

/ 
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and 


(  v{L,s)  \ 
v'(L,s) 

/  cos  (aL) 

sin  (aL) 

cosh  (aL) 

—  asia(aL) 

a  cos  (<lL) 

a  sinh  (aL) 

v"(L,s) 

-a3  cos  ( aL ) 

-a3sin(oL)  a2  cosh  (aL) 

\v"(L,s)j 

\  a3  sin  (a L) 

—a3  cos  (aL)  a3sinh(aZ<) 

sinh(aL)  \ 


(kx\ 

k7 

k3 

\  ki  , 

w(3) 

(2  r 


The  relation  between  the  boundary  conditions  at  either  end  can  be  expressed  as 
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/  ”(£>•»)  \ 
v  (L. s)  _  1 

v"(L,a )  ~2 

V  v  (L,  3)  / 

(  cos  (aL)  +  cosh  (aL) 
a(sinh(aL)  -  sia(ai)) 

aJ(cosh(ai)  -  cos(ai)) 

V  a3($inh  (aL)  +  sin  (aL)) 


■ink  ( aL )  4-  *in  {aL) 
a 

cos  (aL)  4-  cosh  (aL) 
a(sinh  (aL)  -  sin  (aL)) 
a2(cosh  (aL)  -  cos  (aL)) 


co%h(aL)  —  coM  (aL) 

»inh (aL)-f tin  (aL) 
a 


cos  (aL)  +  cosh  (aL) 
a(sinh  (aL)  -  sin  (aL)) 


»inh  (oL)  -tin  (aL) 

--j 

co*k  (aL)-co«  (aL) 

nnk(aL)+»m  (a  L) 
a 

cos  (aL)  -r  cosh  (aL ) 


/  v(0,s)  \ 
v'(0,s) 
v"(0,s) 
\v"'(0,3)J 


Free  boundary  conditons 


at  both  ends,  plus  a  torque  T(s)  at  the  i  =  L  end,  result  in 


(  u"(0,j)  =  0 

u"'(0,s)  =  0 
|  EIv"(L,s)  =  T(s) 
{  v"'(L,  3)  =  0 


(2.9) 


Solving  these  simultaneous  equations  for  v(L,  s)  and  T(j)  gives  us  the  transfer  function  from 
a  torque  T(s)  at  the  i  =  L  end  to  the  displacement  at  i  =  0.  Taking  L  =  1  we  obtain 


where 


P(s) 


u(0, j) 

T(s) 


cos  (o)  —  cosh  (a) 
£/a2(cos  (a)  cosh  (a)  -  1) 


3*  +  -3 


a  =  — 


El 

m 


+  c-t 


(2.10) 


(2.11) 


2.1.2  Essential  Singularities 

The  appearance  is  that  the  functions  appearing  in  the  numerator  and  denominator  of  (2.10) 
may  have  essential  singularities  both  in  the  left  half  plane  and  on  the  imaginary  axis,  because 
the  arguments  to  these  functions  are  roots  of  a  rational  function  of  s. 

We  shall  show  that  we  can  factor  the  transfer  function  as  the  quotient  of  two  functions, 
each  of  which  is  the  composition  of  a  holomorphic  function  with  a  rational  transformation  of 
the  complex  plane.  Still,  numerator  and  denominator  are  not  analytic  functions  of  s.  This 
rational  transformation  has  a  pole  in  the  left  half  plane,  and  we  shall  show  that  the  pole 
becomes  an  essential  singularity.  Although  the  essential  singularities  for  both  the  numerator 
and  denominator  are  at  the  same  point  in  the  complex  plane  there  is  not  cancellation. 

Using  power  series  expansions  we  find  that 

(2  8  10  \ 

T\  +  +  fol  + '  ■  ■)  (212) 
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and 


Thus 


2*  24  26 

cos  (a)  cosh  (a)  -  1  =  -  —  a*  +  —a8  -  —  a1J  +  . .  . 


4!  8!  12! 

cos  (a)  -  cosh  (a) 


(2.13) 

(2.14) 


is  a  holomorphic  function  of  a4,  and  a4  is  in  turn  a  rational  function  of  s,  as  claimed.  Similarly 
for  D(a)  =  EI(cos  (a)  cosh  (a)  —  1). 

This  allows  us  to  express  the  transfer  function  in  the  form 


(2.15) 


Since  numerator  and  denominator  are  holomorphic  functions  of  a4,  an  essential  singularity 
could  only  occur  at  infinity,  and  it  is  easy  to  see  that  this  in  fact  occurs.  Furthermore,  this 
essential  singularity  of  numerator  and  denominator  is  not  a  removable  singularity  of  the 
transfer  function,  i.e.,  there  is  not  cancellation  between  numerator  and  denominator.  One 
can  see  this  by  observing  that  (2.10)  has  no  limit  as  a-*oo,  say  along  the  real  axis.  One 
can  check  that  a/+oo  for  values  of  s  in  the  right  half  plane:  a4  becomes  unbounded  only 
as  3 — ►oo,  and,  if  c,  ^  0,  as  s— *  —  f-  .  These  are  the  essential  singularities  of  the  transfer 
function. 

Figure  2.1  shows  the  behavior  of  (2.10)  as  s— >±oo  along  the  reed  axis. 


-22  I _ i _ i _ I _ i _ i _ i _ I 

-1500  -1000  -500  0  500  1000 


Figure  2.1:  |F(s)|  as  s-* oo  along  the  real  axis,  (m  =  1,  E  =  .1,  /  =  .1,  c  =  .001,  c,  =  .001) 


2.1.3  Analysis  of  Poles  and  Zeros 

The  analysis  continues  by  finding  the  zeros  and  singularities  of  the  numerator  and  denomina¬ 
tor  of  the  plant  transfer  function.  For  this  purpose  we  shall  first  regard  the  transfer  function 
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as  a  function  of  a.  So  we  factor  the  beam  transfer  function  as  P(s)  =  with  A'(a)  and 

D(a )  as  above. 

The  zeros  of  N(a)  axe  given  by  a  =  (l±j)mr,  for  n  =  1, 2, . . .,  or  a4  =  —  4(nir)4.  Notice 
that  the  zero  corresponding  to  n  =  0  is  cancelled  by  the  denominator  in  (2.14).  Taking 
an  =  — 4(n7r)4,  we  note  that  the  values  of  a-*  are  real  and  negative. 

Let  (3n  be  the  values  which  a4  assumes  at  the  zeros  of  D(a).  We  can  show  that  the  values 

/ 3n  are  all  real  and  positive,  and  asymptotically  approach  the  values  Obviously 

a4  =  0  is  also  a  zero. 

In  terms  of  the  Laplace  transform  variable  s,  using  (2.11)  we  find  that  the  numerator 
zeros  are  given  by 


z*  n 
z2n+l 


} 


-( OnC.I  +  c)±\J{<XnC§I  -F  c )3  -  4maw£J 
2m 


(2.16) 


for  n  =  1,2,...  . 

The  denominator  is  zero  at  s  =  0,  and  at  the  points  s  =  -  ^  and 


p2n  1  ~(Pnc,I  +  c)±\J(j3nc,I  +  c)3  -  4 m0nEI 

Pan+i  I  2m 


,  n 


1,2,... 


(2.17) 


which  are  all  in  the  left  half  plane.  Both  numerator  and  denominator  are  unbounded  in  every 
neighborhood  of  the  point  s  =  —  and  have  zeros  dense  at  that  point.  This  is  consistent 
with  the  fact  that  the  point  is  an  essential  singularity  of  the  transfer  function  as  mentioned 
above.  Figure  2.2  shows  a  typical  pattern  of  poles  and  zeros  of  the  transfer  function,  with 


s.  =  100. 

c. 


100 
50 
0 
-50 

-800  -600  -400  -200  0  200  400  600 

Figure  2.2:  Poles  ard  Zeros  of  Beam  Transfer  Function  (m  =  1,  E  =  .1,  /  =  .1,  c  =  .001,  c,  =  .001) 
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2.1.4  Factorization  of  the  Transfer  Function 

Since  ^P(s)  £  H°°  fl  C,  the  inner  factor  must  be  the  product  of  a  Blaschke  product 
and  possibly  an  exponential  of  the  form  t~‘a  with  a>0.  But  e'a^yP(.s)  is  unbounded  on 


the  positive  real  axis,  so  a  =  0.  Thus  the  transfer  function  can  be  factored  as  P(s)  = 
Pu(a)Pl(a)P0(a),  where  Pu(s)  =  i*1, 


P>(*) 


n 


a  -  zn 
a  +  zn 


U  ~  **\ 

1  ~zl 


(2.18) 


and  P0(a)  6  L1  fl  C.  P0  has  essential  singularities  at  the  points  s  =  oo  and  a  =  Thus 
the  outer  part  of  the  beam  transfer  function  is  indeed  irrational,  incorporating  the  essential 
singularities. 

Following  the  technique  in  [Callier  and  Desoer,  p.  655],  we  can  readily  compute  a  coprime 
factorization  of  this  plant  as 


Pi.)  = 

t+ 1 


(2-19) 


and 


p,(0)Po(0) 


P,(s)P0(s)  + 


(3  +  1) 


3  +  1 


P;(0)Po(0) 


=  1 


(2.20) 


2.1.5  Asymptotic  Behavior  of  the  Transfer  Function 

We  now  examine  the  behavior  of  the  transfer  function  at  infinity,  i.e.,  whether  the  transfer 
function  is  strictly  proper.  Our  main  point  here  is  to  show  that  it  decreases  at  infinity  in 
the  closed  right  half  plane  faster  than  any  polynomial. 

The  demonstration  of  this  fact  follows  two  observations: 


1.  The  image  of  the  right  half  plane  under  the  map 


/  ma3  +  ca  \  4 
V  EI  +  c,Ia) 


(2.21) 


satisfies  |Im(a)|<Re(a)  and  approaches  the  cone  |  <|  arg(a)  |<  ^  in  the  right  half 
plane  uniformly  in  arg(s)  for  large  |s|.  As  Is]— >oo,  |o|— »oo. 

2.  Let  a  =  x  +  jy  with  z,  y  €  R.  Then 

(8  cos(z  +  y)  cosh  (z  +  y)  +  8  cos  (z  -  y)  cosh  (z  -  y)  \ 

-  4  cosh  (2y)  -  4  cos  (2y)  -  4  cosh  (2z)  -  4  cos  (2z)  / 

- Lli - - — // 

1  1  /cosh(2z  +  2y)  +  cos  (2z  +  2y)  +  cosh  (2z  -  2y)  +  cos  (2z  -  2y)' 

-  8  cos  (z  -  y)  cosh  (z  +  y)  -  8  cos  (z  +  y)  cosh  (z  -  y) 

^  +  2cos(2y)cosh(2y)  +  2cos(2z)cosh(2z)  +  8  y 

(2.22) 


From  observation  (1)  any  a  in  the  image  of  the  right  half  plane  satisfies  y  =  r0z,  with 
ra  <  1.  Again  by  observation  (1),  for  sufficiently  large  |s|,  ra  is  bounded  from  zero.  Then 
using  observation  (2)  it  is  easy  to  check  that  for  a  sufficiently  large  in  the  right  half  plane 
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|.P(s(a))|3  is  bounded  above  by  a  positive  constant  times  e~3r *  for  some  fixed  r  >  0.  This 
implies  that  any  polynomial  in  |a|  times  )P(s)|  goes  to  zero  at  infinity  in  the  right  half  plane. 
Without  much  more  work  one  can  in  fact  check  that  |.P(jiu;)|  is  asymptotically  equal  to 


i  naz 

E\j  ml\w\ 


.  c-*Vh 


(2.23) 


where  k  =  sin  This  can  be  checked  in  the  plots  in  Figures  2.3  and  2.4  (which  use  the 
same  parameters  for  the  beam  model  as  above). 


Figure  2.3:  | P  (s)|  along  the  imaginary  axis. 


Figure  2.4:  log|F(s)|  along  the  imaginary  axis. 


2.2  A  Damped  Two-Degree-of- Freedom  Beam 

The  model  for  a  flexible  beam  with  a  tip  body,  including  both  bending  and  torsional  vi¬ 
brations  consists  of  a  pair  of  decoupled  partial  differential  equations  with  coupled  boundary 
equations.  The  partial  differential  equations  are  [41]: 


dt3  p  dtdxA  p  dx 4 


(2.24) 
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(2.25) 


d2<f>(x,t)  2SGJ  &<Kxtt)  GJ  d24>{x,t) 


di2 


dx1 


=  0 


pK 2  dtdx 2  pK. 2 

where  y(x,t)  represents  the  transverse  displacement  of  the  beam  and  <j>(x,t )  represents  the 
angle  of  twist  of  the  beam.  The  angle  of  rotation  of  the  motor  is  given  by  0(t);  consequently 
9"{t )  represents  the  angular  acceleration  of  the  motor. 

Taking  Laplace  transforms  of  each  equation  and  using  zero  initial  conditions  (since  the 
transfer  function  is  derived  from  a  zero  resting  state)  yields: 


F  T 

S2Y{x,  3)  4  —  (26s  4  l)rW(*f  j)  =  -X33Q(3) 

P 

GJ 

s2$(x,s) - j(2£j  4  l)$^(x,.s)  =  0 

p  K 

Both  equations  are  separable  and  therefore  have  solutions  of  the  form: 

F(x,s)  =  to(a)v(x) 

*(*.«)  =  /3(a)a(x) 

Using  these  relationships  in  equations  2.26  and  2.27  yields: 


(2.26) 

(2.27) 


v^(x) 

-s2p 

(2.28) 

v(x) 

EI(26s  +  l) 

cS2\x) 

s2pn2 

(2.29) 

a(x) 

GJ(26s  4  1) 

Note:  the  solution  given  is  for  the  associated  homogeneous  equation  to  equation  2.26;  equation 
2. 27  is  already  homogeneous 

Solutions  to  the  equations  2.26  and  2.27  we  therefore  of  the  form: 

Y(x,  s)  =  (k\  cos(\\x)  4-  fcj  sin(Aii)  4  k3  cosh(Aix)  4  sinh(A1x))u>(s)  4  yp(x,  s)  (2.30) 

$(x,  s )  =  (aieA,x  +  a2e~Xl*)(3(s)  (2.31) 

where  yp(x,s)  is  a  particular  solution  of  the  nonhomogeneous  equation  2.26.  A  solution  to 
equation  2.26  is  given  by  yp(x,s)  =  —x©(s)  so  we  have: 

Y(x,s)  =  ( ki  cos(Aix)  4  fcjsin(Aix)  4  &3Cosh(Aix)  4  fc4  sinh(Aix))to(.s)  -  i0(s)  (2.32) 

The  boundary  conditions  are  related  to  the  ki  and  a*  by: 


/  ^(0,3)  \ 

/ 

Y'(0,s) 

Y"(0,s) 

\  Y"'(0,  s)  J 

\ 

1  0 

0  Ai  0  \i 

-A2  0  A2  0 

0  -A*  0  \\  ) 


1  0  \  (  kx\ 

k2 
k3 

\kAJ 


w(s)  - 


(  o\ 
1 
0 

V  o  } 


©(-) 
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(2.33) 


and 


/  Y(L,s)  \ 

/  cos(Ai  L) 

sin(A1L) 

cosh(A1I,) 

sinh(AiL)  \ 

Y‘(L,3 ) 

— Ai  sin( Ax  Z/) 

Ai  cos(AiL) 

\i  sinh(AiZ,) 

Ai  cosh(A1L) 

Y"{L,s) 

— Aj  cos^AjL) 

—  A*  sin^Ii) 

A*  cosh(A!.L) 

A*  sinh(AiL) 

\  Y'"(L,s)  j 

^  A^sin(AiL) 

—  A*  cos(A1L) 

Aj  sinh(AiL) 

AiCosh(Ai  L)  j 

f  k'  \ 

k-i 

k$ 

w(a)  - 

(  L  ) 
1 

0 

Q(s) 

(2.34) 

\  k4  ) 

l  0  j 

for  the  transverse  displacement  and 


and 


(  *(0,3) 

V  *'(0,3) 


(  *(*,-)  ^  _ 
V  )  ~ 


gAjL  e -AjL 

AaeAli  — Aae-A,i 


(2.35) 

(2.36) 

for  the  angle  of  twist. 

Next  we  combine  equations  2.33  and  2.34  to  eliminate  the  fc,  and  w(s),  resulting  in: 


/  Y(L,a)  \ 

( 

Y'{L,  s) 

Y"(L,s) 

\  Y"'(L,3 )  yi 

i 

«inh(Ai.&)+«ui(Ai&)  r  \ 

eo«h(AiD+co«(AiI»)  ^ 

iinh( 

,  2 

Aj  coah(A1I.)-co*(AiL) 

2  / 


1  coHAtLJ  +  co^Aii) 

Al(unh(Al£.)  -  iin(AiL))  coa(Ai&)  +  co*h(AtI>) 

Aj(co«h(AiL)  -  cotiXiL))  Ai(iinii(AiIi)  -  linfAtL)) 
\  Aj(imh^Att>)  +  iin(AiIi))  A^cod^A,  £,)  -  coa(AiL)) 


coth(  t)~coi(  L) 
«iah(  AiL^4»aia(  Xil>) 

CO«(AlL)  +  CO«h(AxIr) 

Xi{amb(X\L)  -  tin(Ai L)) 


/  y(o,j)  \ 
Ho,  3) 
no,  3) 

v  no,  3) 


»inh(  Xy  L)—sin(  Ax  L)  \ 

ij  \ 

coah(  A i  L)^-co>(  Xi  L) 

x\ 

*iah(  XiLjt-f  tia(  Xx  L) 

cot(\iL)  +  coth(XiL)  ) 


(2.37) 


Similarly,  combining  equations  2.35  and  2.36  to  eliminate  the  Oi  and  (3(a) 
(  #(Z,j)  \  _  1  /  eXlL  +  e~x*L  ±(e*’L-e-x’L  \  f  #(0,3)  ' 

\  #'(£,3)  J  2  \  A3(eA,i  -  e"Aji)  ex*L  +  e~x*L  )  \  *'(0,j)  , 

The  simple  boundary  conditions  are  given  by: 


we  obtain: 

|  (2.38) 


no,*)  =  o 

HO,  a)  =  0 
#(0,s)  =  0 


(2.39) 
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while  the  coupled  boundary  conditions  are  [41]: 

0  =  ma2[(L  +  c)0(a)  +  Y(L,  a)  +  cY'(L,  a)  +  e$(Z,  a))  -  £1(1  +  2 6a)Y"\L,  3X2.40) 
0  =  mca2[(L  +  c)Q(a)  +  Y(L,a)  +  cY'{L,a)  +  e$(L,3)} 

+  Jo*2[Q(a)  +  Y‘{L,  j)]  +  £/(l  +  2 Sa)Y"(L,  a)  (2.41) 

0  =  Tnea2[{L  +  c)Q(a)  +  Y(L,a)  +  cY,(L>a)  +  e$(L,a)} 

+  Jea2$(L,  a)  +  GJ(  1  +  2 6a)V(L,  a)  (2.42) 


Now,  using  equations  2.37  and  2.38  with  the  boundary  conditions  given  above,  we  have 
9  equations  in  10  unknowns.  All  that  remains  to  be  done  to  find  the  transfer  function  is  to 
eliminate  8  of  the  variables.  There  are  in  fact  two  transfer  functions  of  interest  in  this  system. 
One  is  the  transfer  function  from  0(s),  the  angle  of  rotation  of  the  motor,  to  Y(L,a),  the 
transverse  displacement  at  the  end  of  the  beam;  the  other  is  the  transfer  function  from  0(s) 
to  $(Z,s),  the  angle  of  twist  at  the  end  of  the  beam. 

Using  Macsyma,  these  transfer  functions  were  found  to  be: 


Y(L,a)  Nx(a) 

ew  m 


(2.43) 


fCM 

e(j)  dm 

where  Ni(a),  N7(a)  and  D(a)  are  given  by: 


(2.44) 


Nx(a)  =  ~[a3J.  sinh(AjZ)  +  X2GJ(26a  +  1)  cosh(A2Z)]{A3£3J3(2Ss  +  l)J[sinh(A1I)  +  sin(AxI)  - 
A1£(cos(A1ii)  cosh(AxZ)  +  1)]  -  2X\EIL{26s  +  l)mcs3(sin(AxZ)  sinh(AxZ))  + 

XlEI(2Sa  +  l)mc33(sinh(A1I)  -  sin(Ai Z))  +  XiEIL(2Sa  +  l)ms3(cos(A1X)  sinh(A1I)  - 
cosh(A1X)  sin(AxZr))  +  X\EI{26a  +  l)s3(mc3  +  /0)(co»h(A1I)  -  cos(A XZ))  - 
X\EIL{26a  +  l)s3(me3  +  J0)(cot(XlL)  sinh(A1£)  +  cosh(A1h)  sin(AxZ))  + 
JoLma*(coa(XiL)  cosh(AxZ)  -  1)}  -  A3£/sinh(Aji)s3(2^s  +  l)me3{s3Jo[cosh(A1h)  - 
cos(A iL)  -  XiL(cot(XiL )  sinh(AxZ)  +  cosh(AxZ) sin(AxZ))]  +  XiEI(2Sa  +  l)[sinh(AxZ)  + 
sin(AiX)  -  A1I(cos(AiI)cosh(AiI)  +  1)]} 

(2.45) 

Nj(s)  =  A3£/sinh(A2Z)mes3(2£j  +  l){/o33(cosh(AxZ)  -  cos(AxZ))  +  (26s  +  1)A \EI  ■ 

[sinh(AiIi)  +  sin(AxZ)  +  A1c(cosh(A1I)  +  cos(AxZ))]} 

(2.46) 

D(a)  =  [s2J.  sinh(A2Z)  +  X2GJ{26a  +  1)  co.h(AaI)]{ms470(cos(A1X)  cosh(AxZ)  -  1)  - 

(2 6a  +  l)3A}£3/3(cos(A1I)cosh(A1£)  +  1)  +  ma3X1EI{26a  +  l)[cos(AxZ)  sinh(AaZ)  - 
cosh(Aji)  sin(Aii)  -  2A1csin(AiI) sinh(AxZ)]  -  (me3  +  Jo)X\EIa3(2Sa  +  1)  • 

(cos(AxL)  sinh(AxZ)  +  cosh(AxZ)  sin(AxZ))}  -  me2a3  sinh(A2Z)Ax-E/(2£3  +  1)  • 

[■ X\EI(26a  +  1)(1  -f  cos(Ax^)cosh(AxZ))  +  Joa3(coa(XxL) sinh(XxL)  +  cosh(XiL) sin(AjZ))J 

(2.47) 
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2.3  A  Damped  Timoshenko  Beam 

The  partial  differential  equation  model,  including  damping  terms,  for  the  Timoshenko  beam 
is  given  by: 

d4v  c,  dsv  /  p  p  \  d*v  dv  p  d3v  p3  d*v  o 

dx*  +  £  dx*dt  ~  \~KG  +  EJ  d^dt3+Cdt  +  ~EK3dt3  +TEG~di*  =°  ^2'48^ 

let 


c, 

E=C" 


(-L. 

\KG 


+  |)  =  C  = 


=  C<, 


EK3  KEG 

By  taking  the  Laplace  transform  of  the  equation  we  obtain: 

(1  +  cia)v^4^  +  c2a3v W  +  (css4  +  c4a3  +  c3a)v  =  0 


=  c5 


(2.49) 


where  v  =  v(x,a)  and  primes  denote  differentiation  with  respect  to  x.  This  is  an  ordinary 
differential  equation.  Making  the  usual  assumption  that  the  solutions  have  the  form  c(s)e*^* 
and  applying  the  differential  operator  we  have: 


(1  +  Cis)A4  +  cajaAa  +  (c5s4  +  c4j2  -f  c3s)  =  0 


where  the  roots  of  this  fourth  order  equation  yield  four  possible  values  of  A(s).  Solving  this, 
we  discover  that  the  solutions  take  the  form: 


v(x,  a)  =  fci(.s)cosh  Ati  +  fca(j)sinh  Axx  +  k3(a)  cosh  Aax  +  k4(a)  sinhAji  (2.50) 


where 

Ai  = 

Aj  = 


-c2a3  +  (~4c1c8j5  +  (cj  -  4c5)s4  -  4(cic«)j3  -  4(c4  +  cic3)sa  -  4c3s) 

2(1  +  cis) 

c2a3  -  (~4cic8s5  4-  (cj  -  4c8)s4  -  4(cxC4)s3  -  4(c«  +  CiC3)a3  -  4c3j)  » 

2(1  +  Ci  j) 


n  * 

2 


(2.51) 


The  boundary  conditions  at  each  end  of  the  beam  (x  =  0,  x  =  L)  are  related  to  kx (s) . . .  k4(a) 

by: 


v{0 ,a) 

■  1  0  1  O' 

'  A3)  ' 

AO,  a) 

0  Ai  0  Aa 

kt(a) 

v"(0  ,a) 

Aa  0  Aa  0 

k3(a) 

.  v'"(0 ,a)  . 

.  0  A?  0  A3  . 

.  k*U) . 

v(L,a) 

cosh  Ail,  sinhAiL  coshAjL  sinhAjL 

’  A3) ' 

v'(L,a) 

AisinhAiL  AicoshAiL  AasinhAaL  AacoshAa L 

k2(a) 

v"(L,a) 

AacoshAiL  AasinhAiZf  AjCoshAjL  AjsinhAjL 

k*(s) 

Lvw(L,s)J 

AasinhAiL  AfcoshAiL  A3sinhA3Z»  AfcoshAjL. 

.  k4(a)  _ 
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The  relationship  between  boundary  conditions  at  either  end  is: 


v(L,a) 

v(0,s) 

v‘{L,a) 

—  T 

v\0,a) 

v"(L,s) 

—  X 

v"(0  ,a) 

.  v"'(L,a )  . 

v"'(0,  a) 

Aj  cosh  Ai  L  —  A*  cosh  Aa  L 
AiAj  sinh  \\L  —  AaAa  sinh  A  aL 
AjAj  cosh  \\L  —  Aj  Aj  cosh  Aa  L 
AjAj  sinh  XxL  —  Aa  Aj  sinh  A ji 


r2  sinh  AiL  —  t2  sinh  A 2L 
Aj  cosh  \\L  —  Aj  cosh  A 2L 
Aj  Aj  sinh  A XL  —  Aa Aa  sinh  A ji 
A^Aj  cosh  XxL  —  A* A*  cosh  A 2L 


cosh  Aj  L  —  cosh  A  XL 
Aa  sinh  Aa L  —  Ai  sinh  Ai  L 
Aj  cosh  AjX  —  Aj  cosh  XxL 
Aj  sinh  Aj  L  —  A®  sinh  A  XL 


sinh  AjL  —  ^  sinh  A XL 
cosh  Aj  L  —  cosh  A XL 
Aj  sinh  Aj L  —  Ai  sinh  XxL 
Aj  cosh  Aj L  —  Aj  cosh  A XL 


The  boundary  conditions  considered  are: 

u"(0,  a)  -  ^s2v(0,s)  =  ^ 

f»"'(0,*)-[l  +  ife],V(0,,)  =  0 

«•(£.»)  =  W 

y(l1.)-[l+i].V(Ill)  =  0 


The  relationship  between  the  boundary  conditions  given  by  the  matrix  T  along  with  the 
boundary  conditions  themselves  define  a  system  of  eight  equations  in  nine  unknowns.  If  we 
solve  these  equations  to  get  u(0,  a)  in  terms  of  the  torque  T(a),  then  we  have  the  transfer 
function  from  the  torque  applied  at  one  end  to  the  displacement  at  the  other  end.  The 
solution  was  performed  on  Macsyma,  and  yielded  the  transfer  function: 


P(,\  -  v(°’3)  _  Mfl 

(  ^  T(»)  D(a) 


where: 


N(a)  =  -GK((K  +  EG){GK  +  £)((  A*  +  A2)  sinh^L)  sinh(AaI)  + 

2AjAj(1  -  cosh(AxL)cosh(Aj L)))pV  +  /f(-(Aj  +  A \)EG{K  +  EG)  ■ 
(Aj  sinh2(AjL)  +  Ai  sinh2(A XL))  -  (2A2A \G2K{K  +  EG)  +  (A2  +  A2)  • 
(Aj  -  A,Aj  +  A \)EGK  -  A1A,(Aa  -  Xx)2E2G2  +  (A,  +  X*)E2)  • 
sinh(A1L)  sinh(Aji)  +  AxAj(A5  +  A2)(G(tf  +  EG){GK  +  2  E)  + 
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E(GK  +  £))(cosh(AiL)cosh(A2L)  -  1)  +  AiAj(A£  -  A \)G(K  +  EG)  ■ 

{GK  +  £)(cosh(AaL)  -  cosh(AiI,)))pJ34  +  EGK2{{X\  +  A \)E  • 

(A2  sinh2(A2£)  +  A3sinh2(AxL))  +  A2A2((A2  +  X\)G{2K  +  EG)  + 

(Aj  —  A1)2£J)siah(AiI>)sinii(Aa£)  +  AxA2(  —  E{\\  cosh(AjjL)  +  Aj  cosh(AxL))  + 
EG2{X\  cosh(AjL)  +  A4cosh(A XL))  +  A2A2£(1  -  G2)(cosh(A2L)  +  cosh(Ai£))  + 
(A,  -  X*)GK (cosh(XiL)  -  cosh(A 3L)))  -  AxAa((A2  +  X\)3GK  + 

(Aj  +  A  \)EG2  +  4A2Aj£J)  cosh(AxL)  cosh(AjL)  +  AxAa((A2  +  A  \f{GK  +  E)  + 
2A?A \EG2))ps2  +  \\X2E2G2K\-{\\  +  A4)sinh(AxL)sinh(AaI)  +  AxA2  ■ 

(A2  +  A2)(cosh(A1I)cosh(A2I)  -  1)  +  AiA2(A2  -  A2)(cosh(A2L)  -  cosh(A  1L)))) 

D(s)  =  EI((K  +  EG)(GK  +  £)((A2  +  A2)  sinh(AxL)  sinh(AjL)  + 

2AtA2(l  -  co3h(A1I)cosh(A2I)))p438  +  K(-( A2  +  A \)EG{K  +  EG)  ■ 

(A2  sinh2(A2I)  +  Ax  smh2(AxI))  +  (-(A4  +  A \)E{GK  +  E)  - 
4AJA 2G2K(K  +  EG)  +  AxAa(A2  -  4AXA2  +  A \)EG(K  +  EG)) 
smh(AxZ)sinh(A  3L)  +  AxAa(A2  +  A2)(2G2tf(Jr  +  EG)  + 

4EGK  +  ZE2G2  +  £2)(cosh(AxI)cosh(A2£)  -  1  ))pV  +  Gtf2((A3  +  X3)  • 

E(E( X3  sinh2(A2 L)  +  A2  sinh2(Ax£))  +  X1X3G{K  +  EG)  ■ 

(Ax  sinh2(A2£f)  +  A2  sinh2(Ax£)))  +  X1X3(XlX3(X 2  +  A \)G2K{K  +  EG)  - 
(A4  -  4AiA2(A2  +  A2)  +  X\)EGK  -  (A4  -  2AXA^  -  2A?Aa  +  A \)E2G2 
+2Ax  A2(A2  -  AxA2  +  A \)E2)  sinh(AxL)  sinh(AjL)  +  X1X3(2X2lX2G2K(K  +  EG)  + 
2(A4  +  3A?A2  +  A \)EGK  +  (A4  +  4A2A2  +  A \)E\G2  +  1))  • 

(1  -  cosh(AxL)cosh(A2L)))pV  +  A2A2£G2tf2(-(A3  +  A3)£  • 

(A2  sinh2(AjL)  +  Ax  sinh2(AxL))  -  (2A2A2G(2iC  +  EG)  + 

(A4  -  Ax  A2  +  2AjAj  -  A2Aj  +  A \)E)  sinh(AxL)  sinh(Aj£)  +  AjAj(A2  +  A2)  • 

(G(2K  +  EG)  +  3£)(cosh(AxL)  cosh(AjL)  -  l))ps2  +  X\X\E2G3K*  ■ 

((Aj  -f  A2)sinh(AxL)sinh(AjIr)  +  2AxA2(l  —  cosh(AxL)cosh(AjL)))) 

If  we  use  the  relationship  between  Ax  and  Aj,  some  further  simplification  is  possible.  For 
example,  if  we  write  Ax  =  (u(s)  +  u»(s))»,  we  have  Aj  =  (u(s)  —  ti/(s))*.  Then  we  observe: 

A2  +  A2  =  2u(s) 

A2  —  Aj  =  2u>(j) 

A?Aj  =  u2(s)-w2(s) 

X*  +  Aj  =  u2(s)  +  ^2(s) 

A2Aj  +  AjAj  =  AjA2(/\J  +  A,)  =  2u3(j)  —  2uw2(s) 


These  substitutions  can  be  made  to  simplify  N(s )  and  D(s)  as  functions  of  j.  However,  such 
simplification  would  introduce  complicated  functions  of  s  into  the  equation  so  we  prefer  to 
leave  the  representation  in  terms  of  Ax  and  X3. 
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valve  input 


spark  advance 


torque  load 


engine  speed 


Figure  2.5:  Engine  Model  Structure 

2.4  A  Multivariable  Delay  Problem 

In  [46],  the  idle  speed  control  of  a  typical  V6  fuel- injected  engine  with  computer  controlled 
management  system  waa  considered.  The  model  structure  consists  of  delays  and  second 
order  dynamics.  The  valve- to- rpm  transfer  function  contains  no  zeros.  The  spark-to-rpm 
transfer  function  has  one  zero.  The  system  is  described  by  Figure  2.5. 

From  Figure  2.5  we  see 


P  = 


where 
Pi  = 


0,(s  +  0 ,) 


We  primarily  will  be  interested  in  the  transfer  function  from  d  to  e,  which  we  denote  by 


/  0 

1 

ON 

0 

0 

1 

Pi 

Pi 

Ps 

\P1 

Pi 

P3/ 

A 

_  (  Pll  P\  2  \ 
\Pn  P») 


:e-#r',P3  = 


0t(s  +  9 #) 


■K.e 


-IT, 


U  =  j  =  Pn  +  Px*K(  1  -  PnK)-'Pn 


+ 


Pi 

1  -  Pl*l  -  P3*3 


£  ) 

Pa^i  +  PzkiJ 
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Chapter  3 

Outer  Factor  Absorption 


3.1  Introduction 

A  crucial  step  in  the  derivation  of  the  solution  to  H°°  optimal  control  problems  is  to  show 
that  the  outer  factor  of  the  plant  can  be  ignored  at  a  certain  point  in  the  calculation  of  the 
infimal  norm  of  the  H°°  control  criterion.  Specifically,  given  the  expressions 

mlJW-FMQU 


and 

*=  inf  \\W-MQ\U 

QeH® 

with  W  €  L00,  M  6  H°°  inner  and  F  6  H°°  outer,  the  assumpton  is  made  that  =  u,,  and 
computation  proceeds  to  find  m. 

The  validity  of  this  step  has  been  demonstrated  for  rational  plants  [49]  and  for  plants 
with  irrational  inner  part  but  rational  outer  and  unstable  puts  [10], [18]. 

Our  recent  study  of  an  example  [11]  has  lead  us  to  be  concerned  with  plants  having 
irrational  outer  part.  In  the  present  work  we  extend  the  results  in  [10], [18]  to  the  case  of 
irrational  outer  factors  for  the  plant. 

The  organization  of  this  chapter  is  as  follows: 

In  Section  3.2  we  review  how  the  H°°  optimal  weighted  sensitivity  problem  for  unsta¬ 
ble  plants  can  be  reduced  to  a  minimization  problem  which  is  affine  in  a  free  (functional) 
parameter,  via  the  so-called  “Q-parametrization.”  When  the  outer  and  unstable  factors  of 
the  plant  and  the  weighting  function  are  rational,  is  known  that  under  additional  neces¬ 
sary  assumptions  on  the  weighting  function  a  sequence  of  approximating  solutions  can  be 
constructed  by  first  solving  a  simpler  problem  in  which  an  outer  outer  function  has  been 
“absorbed”  into  the  free  parameter,  and  then  approximately  “extracting”  the  outer  factor. 
Proposition  1  summarizes  these  previously  known  necessary  and  sufficient  conditions.  For 
the  beam  model  of  Section  2.1  the  outer  factor  of  the  plant  is  irrational.  In  this  case  the 
standard  construction  is  [18]  does  not  apply.  In  order  to  be  able  to  treat  still  more  general 
plants,  we  next  define  a  general  zero  of  a  function,  and  proceed  to  prove  corresponding 
necessary  conditions  in  Lemma  1  of  Section  3.3. 


In  Section  3.3  we  present  the  main  result  of  this  chapter,  Proposition  2,  showing  that  the 
conditions  of  Lemma  1  are  also  sufficient  under  reasonable  assumptions  on  the  plant.  The 
key  ideas  behind  this  extension  of  previous  results  are  a  construction  to  approximately  invert 
irrational  outer  functions,  and  a  definition  of  a  generalized  zero  of  a  function  which  allows 
us  to  formulate  satisfactory  necessary  and  sufficient  conditions  for  outer  factor  absorption. 

3.2  Problem  Definition  and  Background 

Our  main  result  here  forms  one  step  in  the  solution  of  H°°  problems.  For  clarity’s  sake  we 
first  define  the  overall  H°°  problem  which  forms  the  setting  of  our  result.  This  is  the  basic 
problem  set  forth  in  [52). 

BASIC  PROBLEM:  Given  a  proper  transfer  function  P(s)  (the  plant),  and  a  weighting 
function  W\(s)  €  H°°,  find  the  infimal  norm  over  all  weighted  stable  disturbance  sensitivity 
transfer  functions  attainable  by  stabilizing  proper  feedback,  t.e., 

i§f||Wk(*)-S(-)|U, 

where  C  —  {stabilizing  and  proper  feedback  compensators}. 

In  Figure  3.1,  feedback  connection  for  the  problem  is  illustrated.  5(s)  =  (1  +  PC)-1  is 
the  transfer  function  from  a  disturbance  d(s)  at  the  output  of  the  plant  to  the  closed  loop 
output  j/(s).  The  transfer  function  of  the  feedback  compensator,  C(s),  is  constrained  to  be 
proper  and  stabilizing  for  the  closed  loop. 


Figure  3.1:  General  Feedback  System 

Remark  1  Since  ||Af  •  ^Ijao  —  Halloo  f0T  anV  inner  function  M  and  X  £  H°°,  we  assume 
( as  usual)  without  loss  of  generality  that  W\  is  outer. 

As  is  well  known,  this  problem  has  been  solved  under  various  additional  assumptions. 
Here  we  separate  the  additional  assumptions  into  two  versions  of  the  basic  problem: 

PROBLEM  R  (rational):  Assume,  in  addition  to  the  assumptions  of  Basic  Problem,  that 
P(s)  and  W\(s)  are  rational  functions, 

PROBLEM  I  (general  inner):  Assume,  in  addition  to  the  assumptions  of  Basic  Problem, 
that  Wi(s)  is  a  rational  function,  and  P(s)  =  i/>(s)v>(j)Pr(s),  where 
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•  V’(j)  is  a  general  inner  function, 

•  Pr(s)  is  rational 

•  <p(s)  is  a  general  outer  function  in  H°°  with  inverse  in  H°° 

Remark  2  In  each  of  these  problems,  as  a  consequence  of  our  assumptions  one  can  factor 
the  plant  as  P(s)  =  P,(s)  •  F0(s)  •  Pu(s),  an  inner- outer-unstable  factorization  of  the  plant, 
where  the  unstable  factor  Pu  is  rational,  analytic  in  the  open  left  half  plane,  and  bounded  from 
zero  in  the  open  right  half  plane.  We  can  assume  without  loss  of  generality  that  Pu  =  D~l , 
with  rational  D  €  H°°.  Analyticity  of  D~l  in  the  open  left  half  plane  implies  that  D  has  no 
zeros  there,  i.e.,  D  =  DXD0,  where  Dx  is  inner  and  Da  is  outer  with  zeros  (as  a  function  on 
C)  only  on  the  (extended)  imaginary  axis. 

The  (standard)  first  step  in  the  solution  involves  the  parametrization  of  all  attainable 
disturbance  sensitivity  functions,  which  is  as  follows: 

Theorem  1  Given  a  plani  P  as  above  ( Problems  R  or  I)  with  factorization  as  in  Remark 
2,  there  exists  an  H°°  function  V(s)  3uch  that,  if  the  compensator  C  is  proper  and  stabilizes 
the  closed  loop,  then  Q  =  —fpc  satisfies 

Q  6  H“ 

and 

(1  -  PC)-1  =  D(V  -  P,P0Q), 

with  Px,  Pa  and  D  as  m  Remark  2  above.  Conversely,  ifQE  H°°  then  (1  —  PC)-1  =  D(V  - 
PXP0Q)  G  H°°,  and  C  =  is  a  proper  stabilizing  compensator.  Defining  C  =  {proper 

stabilizing  feedback  compensators},  we  use  this  parametrization  to  conclude 

inf  ||Wk(l  +PC)-1||0O  =  u£J\WxDV -WxDPiPmQ\\-  (3.1) 

=  inf  ||^1Z)0(V'-P,P0Q)!|00  (3.2) 

Proof:  See,  for  example,  [23,  p.  10].  V  is  the  coefficient  function  in  the  Bezout  identity 
giving  a  coprime  factorization  of  P  :  P,P0U  +  DV  =  1,  and  P  =  PXP0/D.  I 


As  is  well  known,  this  theorem  allows  us  to  simplify  notation,  since  (3.1)  shows  that 
finding  a  solution  is  equivalent  to  finding  a  solution  to  a  related  affine  minimization  problem. 
With  this  simplification,  taking 


W  = 
M  = 
F  = 


WXDV 

DXPX 

WlD0P0 


the  problems  become  equivalent  to 


inf  \\W-FMQl 


geH 


(3.3) 


(3.4) 


where, 
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•  for  Problem  R,  W ,  F  and  M  are  rational; 

•  for  Problem  I,  W  and  M  are  not  necessarily  rational,  and  F  is  the  product  of  a 
rational  outer  function  with  a  stably  invertible  irrational  outer  function.  In  this  case 
we  can  assume  without  loss  of  generality  that  F  is  rational,  since  the  invertibility  of 
the  irrational  outer  factor  of  F  allows  us  to  ignore  it  in  the  minimization. 

Previously  known  results  can  be  summarized  as  follows 


Proposition  1  Let  Wi  be  rationed  and  outer,  and  let  P  =  ipipPT  with  ip  inner,  <p  an  outer 
function  in  H00  invertible  in  H°°,  and  PT  rational.  Define 

Mo  =  inf||^1(s).(l  +  P(s)C'(s))-1||0O 

=  inf  ||^1(s)Z)0(s)(l/(s)  -  Pt(s)P0(s)Q(s)) Hoc 
=  inf||W(s)-F(s)M(s)<?(s)|U 


M,  =  mf!|^(s)-Af(s)Q(a)||0O 

=  g‘nHf-  IW-W-M*)  ~  ftWGWIU 
Let  V(F)  =  {z  €  JR  :  F(z)  =  0}.  Then 

M«  >  W{z)  for  all  z  €  V(F) 


if  and  only  if 

do  =  Ml- 

Proof:  This  is  a  special  case  of  Proposition  2  below. 


(3.5) 

(3.6) 

I 


Remark  3  (3.5)  is  obviously  a  necessary  condition  for  (3.6),  so  the  essence  of  Proposition 
l  is  the  sufficiency  of  (3.5)  for  (3.6).  However,  for  more  general  F  and  W,  the  necessity 
of  (3.5)  is  not  certain:  If  P0  is  discontinuous  we  first  need  to  generalize  the  notion  of  zero. 
Then  one  might  consider  the  case  of  W  and  F  both  beign  discontinuous  at  a  generalized  zero 
of  F  Other  possible  pathological  cases  are  mentioned  in  the  conclusion.  In  Lemma  1  below 
we  proved  the  necessity  for  the  relatively  straightforward  case  of  Pa  continuous  at  its  zeros. 

In  considering  H°°  problems  for  distributed  parameter  systems,  it  has  become  clear  that 
the  assumptions  of  Problem  I  are  too  restrictive. 

EXAMPLE:  For  the  damped  Euler- Bernoulli  beam  model  of  Section  2.1  the  outer  part 
has  an  essential  singularity  at  infinity,  and  the  transfer  function  rolls  off  faster  than 
any  rational  function.  (This  is  not  just  a  result  of  poles  clustering  at  oo,  since,  for 
example,  an  infinite  Blaschke  product  may  have  such  a  cluster  point  of  poles,  but  does 
not  roll  off.) 
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In  mixed  sensitivity  problems  with  irrational  outer  factors  [13],  it  is  also  desirable  to  allow 
irrational  weighting  functions,  which  yields  a  problem  equivalent  to  taking  Wx  irrational  in 
the  present  “pure”  sensitivity  case.  Referring  to  (3.3),  we  see  that  for  this  example,  W,  M 
and  F  are  all  irrational.  If  these  irrational  functions  were  continuous  and  bounded  from  zero 
on  the  imaginary  axis,  there  would  be  no  need  for  further  developments.  However,  irrational 
zeros  require  a  different  construction  to  prove  sufficiency  of  (3.5),  and  the  possibility  of 
discontinuity  of  the  outer  factor  (which  does  not  occur  in  this  beam  example)  makes  it 
necessary  to  use  a  different  condition  than  (3.5),  since  W(z)  will  not  be  defined  at  a  point 
of  discontinuity.  With  this  motivation,  we  shall  consider  the  following  generalization  of 
Problem  I: 

PROBLEM  G  (general  inner  and  outer):  Assume,  in  addition  to  the  assumptions  of  the 
Basic  Problem,  that  P(s)  =  ■tp(s)(fi(s)Pr(s),  where 

•  ip(s)  is  a  general  inner  function. 

•  ip(s)  and  Wi  are  general  outer  functions  which  have  only  finitely  many  “generalized 
zeros”  (defined  below)  on  the  (extended)  imaginary  axis  j'R.  For  technical  reasons,  we 
also  assume  that  v*(s)  and  are  continuous  at  their  own  zeros  on  the  (extended) 
imaginary  axis. 

•  Pr(j)  is  rational. 

In  order  tc  give  the  necessary  conditions  which  are  the  natural  extension  of  (3.5),  we 
define  a  general  zero  of  /  to  be  a  point  on  the  imaginary  axis  for  which  |/|  is  not  essentially 
bounded  from  0  on  any  neighborhood  of  the  point.  We  make  this  precise  as  follows: 
Definition:  Let  /  be  a  Lebesgue  measurable  function  on  jR.  For  z  G  jRU  {oo}  and  e  >  0 

take  Nt(z )  to  be  the  open  interval  (z  -  ie,  z  +  tc).  (  If  z  =  oo,  take  Ne(z )  =  (;,  oo)  U  (-oo, 
—  i).  )  Define 

essran(/,  z)=  fle>o  {  essential  range  of  [/(x)|  restricted  to  jVe(z)} 

(the  essential  range  of  /  at  z).  We  say  that  z  is  a  general  zero  of  /  if 

inf  Rf(z )  =  0, 


and  we  write 

V(/)  =  U  e  jR  U  {oo}|  inf  Rf{z)  =  0}  . 

I 


Remark  4  V(/)  consists  of  zeros  of  finite  order  and  certain  singularities  of  f .  It  must  have 
Lebesgue  measure  zero  for  f  to  be  a  boundary  value  function  of  a  function  in  H°° . 

EXAMPLE:  f{juj)  —  e-v^"  has  a  general  zero  at  oo  (taking  the  square  root  to  have  a 
branch  cut  on  the  negative  real  axis).  /(ju/)  is  the  boundary  value  of  the  H°°  function 
/(s)  =  e~'f~‘ ,  which  is  continuous  on  jR  and  has  an  essential  singularity  at  oo. 
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Remark  5  In  this  third  problem  as  well,  as  a  consequence  of  our  assumptions  one  can  factor 
the  plant  as  P(s)  =  Pi(s)  •  P„( 3)  •  Pu(j)  «  »»  Remark  2. 

As  indicated  in  Proposition  1,  solutions  to  the  first  two  problems  (R  and  I  via  (3.4))  can 
proceed  by  first  neglecting  the  factor  F(s)  and  solving  the  following: 

W=  inf J|W  -  iflJIU.  (3.7) 

The  transformation  of  (3.4)  to  (3.7)  we  call  “outer  factor  absorption”  because  of  the 
construction  used  to  obtain  am  approximate  solution  to  (3.4)  from  a  solution  to  (3.7).  This 
idea  first  appeared  in  [49,  p.  591],  where  a  construction  was  given  for  the  case  of  P  a  rational 
function.  For  the  case  of  P  having  rational  outer  and  unstable  parts  but  general  inner  part 
[10,  p.  69]  gave  the  first  such  construction.  In  the  latter  case  it  is  also  shown  in  [10], [18] 
that  a  sequence  of  approximate  solutions  {^„}  to  (3-4)  can  be  constructed  from  a  solution 
to  (3.7)  such  that  they  are  rational  functions. 

Remark  6  The  infimum  in  (3.7)  is  actually  a  minimum  [24,  p .  135}. 

Now  we  can  proceed  to  develop  parallel  results  for  Problem  G. 

3.3  Results 

We  first  state  the  necessary  condition  for  outer  factor  absorption  in  Problem  G. 

Lemma  1  Let  W  €  L°°(jR),  let  F  be  continuous  at  each  point  ofV(F)  .  Take 

/*.=  inf  \\W  —  FMQ\\oo.  (3.8) 


Then 

Ho  >  sup  Rw{z)  for  every  z  6  V(F). 

Proof:  Let  z  €  V(F).  For  each  e  >  0  by  continuity  of  F  at  z  there  is  am  interval  Nt  =  ( z—j-y , 
z  +  jy)  with  7  >  0  such  that  |F(u;)]  <  e  for  every  w  €  Nt.  Take 

<7  =  sup  Rw{z)- 

Suppose,  contrary  to  the  conclusion  of  the  lemma, 

p.0  =  cr  —  a  for  some  a  >  0,  (3.9) 


and  take  Qa  €  H00  to  satisfy 


\\W  -  FMQaWoo  <  Ho+  <^- 
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Of  course  we  always  have 


inf  ||  W-  FMQWoo  =  Mo  <  \\W  -  FMQa H.. 


Take  e  =  a  /  8||M(?a||oo-  Then  for  u;  6  We, 

IM^O'w)  -  F(ju)M(jw)Qa(jw)\  >  |W(ju;)|  -  |F(jw)M(jw)Qa(jw)| 

>  |Wr(i")l  “  |  ae- 


By  definition  of  a,  there  is  a  set  of  non-zero  Lebesgue  measure  St  C  jV€  such  that 
|  W' (c*^) |  >  cr  —  |  for  each  w  G  St.  Thus  for  u  €  St 

Mo +  7  >  W-FMQoU 

4 

>  \W(jtjj)  —  F(ju)M(ju)Qa(juj)\  a.e. 

.  a 

>  a - a.e., 

4 

and  we  conclude  that 

a 

a  ~  2' 

This  contradiction  with  (3.9)  establishes  that  n0>  <r.  I 


Remark  7  When  W(jlj)  is  continuous  at  z,  |W(z)|  =  sup  ^(2). 

We  next  show  that  the  conditions  of  Lemma  1  are  also  sufficient  for  outer  factor  absorp¬ 
tion  in  Problem  G.  To  motivate  our  development,  we  review  the  ideas  behind  the  proof  of 
sufficiency  in  Proposition  1  as  presented  in  [10], [18]  .  The  proof  is  by  construction,  which 
consists  of  two  steps: 

1.  find  a  sequence  {7k}  C  H°°  such  that 

Mfc  -»  Mo  as  k  -»  00,  (3.10) 

where 

Mk  =  \\W  -  TklU  , 

and 

7k(z)  =  0  for  all  z  G  V(F), 

2.  find  a  sequence  fn  such  that 

(a)  /„  =  Fxn  for  some  i„  G  H°°  (F  divides  fn  in  H°°) 

(b)  H/nlU  ->  1  as  n  ->  oc 

(c)  |1  —  fn\  — ►  0  uniformly  on  finite  intervals  of  the  imaginary  axis  which  exclude 
zeros  of  F. 
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We  leave  the  details  to  the  proof  of  Proposition  2  below,  but  here  we  give  the  ideas 
behind  these  two  steps  to  motivate  the  proof. 

Given  a  solution  Q  to  (3.7),  we  would  like  to  factor  F  out  of  it,  and  still  have  F~1Q  6  H°° 
so  that  F~lQ  gives  a  solution  to  (3.8).  However  F-1  has  poles  at  the  points  of  V(F),  and  Q 
will  not  generally  have  zeros  at  these  points,  so  this  does  not  give  a  solution.  Our  approach 
is  to  find  a  sequence  Qn  such  that 


||  W  -  FMQn\\  — » /i„.  (3.11) 

A  naive  approach  is  multiply  Q  by  another  function  so  that  the  product  is  very  close  to 
Q  over  some  union  of  intervals  which  excludes  V(F),  and  then  has  zeros  of  the  appropriate 
order  at  the  points  of  V(F).  The  natural  way  to  attempt  this  is  the  following:  Assuming 
the  desired  zero  is  at  oo  and  has  order  m,  simply  multiply  Q  by  (7^)  •  Letting  n  — ►  oo, 
the  interval  on  which  the  product  is  close  to  Q  grows.  This  works  for  Problem  R,  but  fails 
for  Problem  I. 

The  difficulty  is  as  follows:  In  Problem  I  {W  —  MQ)(s)  has  constant  magnitude  on  the 
imaginary  axis,  and  so  (W  —  MQ)(j<j)  describes  a  circle.  In  general,  the  locus  will  continue  to 
revolve  indefinitely  as  |u/|  increases.  Assuming  for  simplicity  that  W(jui)  approaches  a  con¬ 
stant  as  |ui|  — *  00,  we  conclude  that  MQ(ju)  approaches  a  circle  with  center  W(ioo).  Simply 
reducing  the  magnitude  of  Q  (multiplying  by  the  “roll-off”  function  (7^)  )  as  frequency 
increases  without  regard  to  phase,  can  result  in  the  sensitivity  locus  leaving  the  asymptotic 
circle.  For  a  roll-off  function  having  only  poles  the  minimal  phase  versus  frequency  charac¬ 
teristic  is  that  of  a  one-pole  roll-off.  To  obtain  less  phase  but  still  get  infinite  attenuation  at 
00  one  must  roll-off  more  slowly.  That  is  the  motivation  for  the  roll-off  functions  hn  below 
which  accomplish  Step  1.  This  shows  us  how  to  modify  Q  so  as  to  introduce  zeros  at  the 
points  of  V(F)  yet  approach  arbitrarily  closely  to  fa. 

The  essential  point  is  that  in  the  rational  plant  case,  at  high  frequency  the  product  PQ 
approaches  a  real  constant.  Roll-off  can  therefore  be  introduced  into  the  feedback  loop  at 
high*  frequency  without  regard  for  the  phase,  whereas  in  the  case  of  general  inner  factors  PQ 
has  unbounded  phase. 

In  order  to  factor  F  out,  the  zero  we  introduce  by  roll-off  must  be  at  least  as  steep  as  the 
zero  of  F.  However  Step  1  introduces  only  very  slow  roll  off,  with  vanishing  phase  change  as 
n  — +  00,  in  order  to  allow  for  the  presence  of  general  inner  factors.  As  a  consequence,  we  will 
have  zeros  of  fractional  order,  and  so  the  order  of  the  zeros  will  still  result  in  F~lTn  £  H°°. 
The  roll-off  could  be  fast  after  a  frequency  at  which  the  loop  gain  has  decreased  sufficiently. 
However  one  must  be  careful  that  the  fast  roll-off  does  not  introduce  excessive  additional 
phase  at  lower  frequencies. 

The  idea  of  Step  2  is  that  after  the  gain  of  Tn  is  sufficiently  small,  we  can  tolerate  the 
introduction  of  large  phase  change  by  the  roll-off  function.  The  problem  then  becomes  one 
of  introducing  the  zero  of  appropriate  order  while  constraining  phase  change  to  be  small 
outside  the  interval  where  Q  has  not  been  rolled  off  sufficiently.  In  Problem  I  this  can  be 
accomplished  using  the  same  functions  (77^)  in  the  naive  approach  above,  but  just  taking 
the  pole  to  be  sufficiently  large.  In  the  more  general  case  of  irrational  F,  it  is  not  necessarily 
obvious  how  to  find  a  function  which  has  zeros  having  order  at  least  that  of  the  zeros  of 
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F  (these  latter  zeros  will  generally  be  singularities),  yet  which  is  close  to  1  on  a  specified 
interval  excluding  zeros.  This  construction  for  the  irrational  case,  and  the  measure-theoretic 
arguments  to  accomodate  general  zeros,  are  the  keys  to  the  generalizations  in  this  report. 

Remark  8  The  sequence  {<?„}  satisfying  (8.11)  directly  gives  a  sequence  of  stabilizing  com¬ 
pensators  via  Cn  =  i-+%nQit  • 

For  simplicity  of  exposition,  we  start  with  the  case  of  a  single  irrational  zero  on  the 
extended  imaginary  axis.  Without  loss  of  generality  we  assume  this  zero  is  at  oo.  (Use  a 
conformal  map  to  move  the  zero  to  any  finite  point.)  In  parallel  to  [18],  our  construction 
here  consists  of  two  steps,  which  we  state  here  as  lemmas.  Step  1  is  accomplished  using  the 
following: 

Lemma  2  Suppose 

H  \\W  -  MQU  =  \\W  -  MQU  (3.12) 

with  Q  6  H°°,  and 

Hi  >  sup  Rw(oo).  (3.13) 

Let  hn^s)  =  ^  a  >  0,  7  >  0  .  Then  the  sequence  Tn  =  hnQ  satisfies  Hn  —►  Hi 

where  Hn  =  \\W  -  MTnWoo. 

Proof:  See  appendix.  A  simpler  version  which  is  suitable  for  treating  Problem  I  is  given  in 
[10,  pp.  69-71]  and  [18,  pp.  519-521].  In  these  references  it  is  shown  that  Tn  can  be  taken 
to  be  rational  functions.  The  parallel  result  is  also  true  in  the  present  case,  but  we  do  not 
prove  it  here.  In  Problem  G  we  do  not  assume  continuity  of  W  at  the  zeros  of  Pa ,  and  so 
the  proof  is  more  involved.  | 


Corollary  1  Assume 

MiC;RU{oo}. 

Suppose  (8.18)  holds,  and  that  Q  solves  (8.12).  Let 

1  l/n 

with  a  >  0,  7  >  0. 


m*)  =  n 

;=i 


3  —  Zi 


3  + 


r+i 


y-*i  J 


Then  the  sequence  Tn  =  hnQ  satisfies  Hn  -»  Hi  where  Hn  =  II W  -  MTn\\ 

Proof:  The  proof  differs  from  that  of  Lemma  2  only  in  the  details  of  the  estimates.  hn  here 
is  the  hn  function  in  Lemma  2  composed  with  conformal  maps  which  take  the  points  of  V 
to  oo.  ■ 


For  Step  2,  in  [10,  pp.  69-71]  and  [18,  pp.  521-522]  we  use  the  obvious  construction 
for  the  rational  case:  we  can  easily  find  the  order  of  the  zero  at  infinity,  and  construct  a 
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sequence  satisfying  the  conditions  in  Step  2.  For  example,  suppose  that  F  =  (t^y)”1-  Then 

we  can  take  fn  =  (^)"\  and  xn  =  F~lfn  =  ■ 

However,  in  the  case  of  general  irrational  P  it  is  not  necessarily  obvious  how  to  construct 
a  sequence  {/„}  to  have  sufficient  roll-off  at  the  irrational  zeros  on  the  imaginary  axis  as  P, 
and  yet  approach  1  in  the  desired  manner. 

From  Lemma  2,  each  T„  goes  to  zero  at  infinity,  as  desired,  but  the  roll-off  of  the  con¬ 
structed  zero  at  infinity  will  not  generally  be  fast  enough  to  make  F-17*  an  element  of  H°°. 
The  issue  addressed  next  is  how  to  construct  /„  (and  xn)  so  that  the  zero  rolls  off  at  least 
as  fast  as  that  of  P,  and  yet  have  Qn  =  xnTn  satisfy  (3.11). 

Lemma  3  For  any  outer  fund  in  p(s)  G  H00  there  exists  a  sequence  {xn}  H°°  such  that: 
(a)  ||pxn||  -►  1  as  n  ♦  oo, 
and 

(b)  |1  —  p(ju)xn(juj)\  — »  0  uniformly  a.e  on  compact  subsets  of  the  imaginary  axis  which 
exclude  zeros  of  p(i.e.,  such  that  Rp^juj)  is  bounded  from  0  at  every  point  on  the  interval). 

Proof:  We  need  the  following  theorem  [24,  p.  85]  which  we  quote  here  transformed  from 
the  unit  disk  to  the  imaginary  axis  and  with  a  change  of  notation: 

Theorem  2  Let  p(s)  be  an  outer  function.  Then  there  are  functions  {xn}  in  H“  such  that 
|xn(s)p(s)|  <  1 

xn(ju)p(ju)  — ♦  1  almost  everywhere.  | 


We  will  not  prove  this  theorem  here,  but  we  will  use  the  fact  that  in  the  proof  of  this 
theorem,  xn  is  defined  as  follows:  Let 

Un(jv)  =  min(A„t  -  log  |p(jw)|) 

where  each  A*,  0  <  An  <  oo,  is  taken  to  be  large  enough  so  that 

One  then  takes 

[1  t 80  u>3  +  i  1 

-  /  ^(ju,)— V  ■  -—du,  . 

ir  J- *  u>  4-  is  1  +  u>3 

This  {xn}  satisfies  the  theorem. 

This  theorem  directly  gives  us  (a),  and  for  (b)  we  need  only  show  the  uniformity  of  the 
convergence  jl  —  fn\  ->0on  finite  intervals  of  the  imaginary  axis  which  exclude  zeros  of  p, 
with  fn  —  xn  ■  p.  For  this,  we  check  that  on  any  closed  finite  interval  Q  C  jR  such  that 
V  fl  fi  =  0,  for  n  sufficiently  large,  |/„|  =  1  and  arg (/„)  <  £  a.e.  on  the  interval. 


Let  Sn  =  {s|s  €  j'R  and  n  <  -  log  |p(s)|  }.  Then  5n+1  C  Sn,  and  we  take  Sn  =  closur e(5n\ 
(fin  Sn)).  Using  the  assumption  that  Rp  is  bounded  from  0  on  Cl,  we  can  take  no  sufficiently 
large  so  that  for  every  n  >  nn,  meas(Cl  D  Sn)  =  0  and  fifl  5n  =  0.  Then  for  s  e  Cl, 

I/Ml  =  |F(s)|.|*n(s)| 

=  -exp^s)] 

=  1  a.e. 


The  phase  of  fn  is  given  a.e.  on  the  imaginary  axis  by: 

Mg(^'"»  =  -;/Z(rh  +  zh)I°slu,‘)l,1‘ 

_l  f°°  (  ujt  +  1 
7T  J- oo  \(1  +  f2)(u>  —  t), 


log 


Then 


=  -i/fc  ((T-^){J_0)  log 


a.e. 


If  u  €  Cl,  we  can  write  a.e 

larg  (/n(;w))|  <  - 


( 

\ 

1 

ess  sup 
t  €  S„ 

ut  +  1 

f  |log|/«(*f)| 

ir 

OJ  —  t 

An  |  1  +  t2 

y 

dt 


<  - 


ess  sup 
*  e  Sn,  1-6(1 


\ 

vt  +  1 

f 

log  |/n(*OI 

V  -  t 

) 

Js„ 

1  +t2 

dt. 


Now,  E  L1  because  /„  €  H”  [24,  p.  66),  and  Cl  and  Sn  are  disjoint  by  choice  of 

n,  so  the  ess  sup  is  finite.  This  provides  a  uniform  bound  on  |  arg(/„(ya;))|  for  u>  €  Q.  Also, 
Js„  1 1°*i{ntat<il  1  dt  — ►  0  as  n  — *■  oo,  since  meas(Sn )  — *  0  and  |  log  |/n(it)||  is  non-increasing  as  a 


function  of  n  for  fixed  t.  Finally,  supt€Siti  „en  is  non- increasing  as  n  increases.  Thus, 
as  n  — ►  oo,  arg(/„)  — ►  0  uniformly  on  any  closed  interval  which  excludes  zeros  of  |p|.  I 


Remark  9  The  key  step  is  separating  Sn  from  Cl  with  closed  sets,  except  for  a  set  of  measure 
zero,  so  that  the  essential  distance  between  the  two  is  bounded  from  zero.  This  is  what  allows 
us  to  bound  the  phase  of  /„. 

Now  we  can  prove  the  following: 

Proposition  2  Let 

Mo=  inf  \\W-FMQW", 
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where  W ,  F  and  M  6  H°°,  and  F  is  outer  and  M  is  inner.  Let 

M.=  inf  \\W-MQW... 

QeH- 

Assume  that  V(F)  is  finite,  and  that  F  is  continuous  at  each  point  ofV(F)  .  Then 

Ht  >  sup  Rw{z)  for  all  z  £  V(F),  (314) 


if  and  only  if 


lia  =  Hi  . 


(3.15) 


Proof:  By  Lemma  1  (3.14)  is  a  necessary  condition  for  (3.15).  For  the  converse,  it  is  obvious 
that  Ho  >  Hi-  We  shall  construct  a  sequence  Qn  such  that  Ho  <  ||W  ~  FMQn\\  — *  Hi,  and 
therefore  conclude  that  Ho  5-  Mi •  Suppose  first  that  V(F)  =  {oo}.  Let  W, =  sup  fZw'(oo), 

and  assume  Ht  >  Woo.  Take  hn(s)  =  [7^]  ^  with  a  >  0,  7  >  0  and  Tn  =  hnQ  as  in 
Lemma  2,  where  Q  solves  (9)  .  Let  {xn}  he  the  sequence  resulting  from  substituting  F  for 
p  in  Lemma  3,  and  we  take  Qn  =  xnTn.  We  shall  show  that  \\W  —  FMQn\\  —*  Ht  by  finding 
conditions  on  n  so  that 

\\W-FMQn\\<l*i  +  e.  (3.16) 

For  €  >  0,  for  each  n  we  claim  we  can  take  wn  >  0  to  satisfy  for  |w|  >  u/n  both 


\W(ju>)\  <  Woo  +  e  a.e. 


(3.17) 


and 


l^n(;w)|  < 


H%-\W(ju)\  +  e 
H*  +  \W(jv)\  °'e' 


(3.18) 


To  verify  this  claim,  note  that  (3.17)  will  hold  by  definition  of  W*  and  by  hypothesis  on 
Hi  for  |u>|  sufficiently  large  independently  of  n.  For  (3.18)  to  hold  it  is  easy  to  check  that  a 
sufficient  condition  is 


l^nl  >  01 


Hi  +  \W(ju)\ 

Hi  ~  \W{j<jj)\  +  e 


a.e.  on  (u>„,  00). 


This  will  always  hold  for  large  enough  n  because  \W(jw)\  is  essentially  bounded  at  00,  and 
\W(juj)\  —  e  is  essentially  bounded  from  Hi  at  00. 

Let  {/„}  be  a  subsequence  of  {Fxn}  such  that 


|1  -  /n(;w) |  <  — for  |u/|  <  u /„.  (3.19) 

n 

The  way  we  pick  n  will  be  to  pick  it  large  enough  so  that  for  |w|  <  u/„(13)  holds.  This 
works  because  for  |cj|  >  u>„ ,  our  definition  of  w„  will  ensure  (13)  holds.  We  see  this  as  follows: 
We  know  that  Q  satisfies  the  condition  that  (W  —  P,Q)  is  all-pass,  say 


Hxeia{u)  =  (W-  MQ)(jw)a.e. 


(3.20) 
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Recalling  the  definition  Xn(jw)  =  W{jiu)  —  Pt(ju)Tn(ju),  assuming  we  have  chosen  un 
to  satisfy  (3.17)  and  (3.18),  and  taking  Hi  and  a(u>)  as  in  (3.20),  we  have 

Hi  ~  \W(ju)\  +  €  >  +  \W{ju)\)a.e. 

\fn{jv)\  <  1  a.e.,  so 

-  W(ju) I  <  \hn(ju>)\  •  -  W{ju>) I 

<  IMiw)l(^  +  l^0'w)l) 

<  —  |W(ja;)|  +  e  a.e. 


But 

(^(jw)  -  F{juj)M{ju)xn{ju)Tn{juj))  =  W{joj)  +  /n(jw)hn(ju;)(/z,e,“(l*,)  -  W(ju;))  a.e., 


so 

\{W(ju)  -  F(jW)Af(;u;)xn(ju;)rn(ju;))|  <  I^MI  +  |/n(jw)An(ju;)(/ile*a(1*')  -  W(jw))| 

<  Hi  +  e  a-e- 

So  now  we  treat  the  case  \lj\  <  un,  essentially  repeating  the  argument  in  Lemma  2.  Let 

Yn(jw)  ^  W(jW)  -  M(ju>)fn(ju)Tn{ju>) 

and 

5n(jw)  =  fnijv'lhnijv). 


\Yn(ju>)\3  =  \W{Ju)  +  pn(ju;)[/i,eM*^  -  W{ju>)]\' a.e. 
Set  ujm  e  [0,  u/n]  to  a  frequency  at  which 


(  \ 


ess  sup  |Kn(jw)l 
I  w  6  [0,  «»] 

V 


€  Ry(wm). 


Now  let  7  =  arg(yn)  and  define 

g  =  sup  Rjn(ujm),  6  =  sup  Ry(ujm),  Wm  =  sup  Rw{um),  and  a  =  a(u/m). 

These  are  functions  of  n,  as  is  o»m.  One  can  see  from  (3.19)  that  |  arg[/n(<7u;m)]|  <  1/n  , 
and  from  the  definition  of  h„(s)  that  0  <  arg (A„(ju;m)l  <  27r/n.  Therefore  -l/r.2  <  f>  < 
2tt /n  +  l/nJ. 

Now  we  show  that  we  can  pick  n  as  a  function  of  e  >  0  such  that  ||Tn||  <  H*  "h  e-  ’n 
the  proof  of  Lemma  2,  |y^fyw)|S  <  [iym(l  —  y)  +  <7/*t]3  +  2 gniWm6  a.e.  Given  n,  there  are 
two  possibilities: 
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case  (i).  (u/m  >  n)  Exactly  as  in  the  proof  of  Lemma  2  we  find  \Yn{j<jj)\2  <  +  e  a.e. 

for  sufficiently  large  n. 

case  (ii).  (o>m  <  n)  Then  since  || W(_;u;)|j00  >  Wm  and  g  <  1,  as  before  we  have 
(sup  RyM)2  <  £  +  (1  -  g)  •  2||^(ju,)||oo/x,  +  +  (l  _  <7)2||VU(ju,)||^0  a.e. 

71 

and 

1  -  9  <  |1  -  5n(;o;m)j  =  |1  -  fn(jwm)hn(jLJm)\ 

Using  the  estimate  of  (1  —  h)  from  the  proof  of  Lemma  2  along  with  (3.18)  we  find 
1  —  g  — *  0,  and  we  conclude  that  —  H^Hoo  — ♦  0  as  n  — ►  oo. 

To  handle  the  case  of  V(F)  finite  with  more  than  one  point,  note  that  Corollary  1  and 
Lemma  3  already  apply  to  this  case.  Use  the  above  to  find  conditions  on  n  for  each  element 
of  V  after  transforming  it  to  oo,  as  in  Corollary  1,  and  take  the  maximum  of  the  resulting 
values.  I 


This  proposition  extends  readily  to  the  case  of  mixed  sensitivity,  which  amounts  to  the 
following  problem  [14]:  Solve 


A 


inf 
ze  H<* 


W  -FMZ 
V 


with 

w(s)  ±  w*(s)w1(s)w;(s)-1r-(s)  -  w2(3)r(s)-1n(s)X{s), 
M(s)iiV,(s)Z},(s), 
F{3)^N0(3)D0(3)W3{s)R-1{3), 
V(s)^Wl(s)R(s)t 

and  R  €  H°°  and  outer  satisfies 


(3.21) 


w;{s)w2(s) 

wyWWiW  +  W7(*)W,(3)  h  [  m  h 

Here  Nt,  Dt  6  H°°  are  inner,  Not  Da  and  W2  G  H°°  are  outer,  and  W2  is  outer  (but  not 
necessarily  in  H°°). 

Corollary  2  Let  F,  M,  W,  and  V  be  as  above.  Assume  W2  is  continuous  jR,  and  F  is 
continuous  on  V(F),  which  has  finite  cardinality.  Define 

W-  MZ 

V 


n%  =  inf 
zeH- 


Assume  further  that 


>  supOIVIjoo,  [Wi(;u;)|  :  ju>  6  ae(M)}. 


(3.22) 
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Take 


Ho  =  inf 
zeH° 


W-MFZ 

V 


Then 

if  and  only  if 


Ho  =  Hi 


Hi  >  SUP  Rww+vv(z). 

x6  V(f) 


Proof:  The  assumption  (3.22)  allows  one  to  transform  the  mixed  sensitivity  problem  to  a 
sensitivity  problem  [44].  The  calculation  is  as  follows:  For  Z  €  H°° 

H 2  >  (W(ju)  -  M{ju)Z{iu)Y{W{juj)  -  M(ju)Z{ju>))  +  V\ju)V{ju)  a.e. 

From  (3.22)  Hi  ~  V*(ju)V (ju>)  >  0,  so  we  cam  find  an  outer  spectral  factor  T  such  that 

T-{jw)T(iu)  =  -  V(jw)V(jw). 

Then  we  have 

1  >  \T-\ju,)(W(jv)  -  M(j«)Z(i«))|. 

Of  course,  we  must  have 


L  =  z^JT-\W  -  MZ)\\e 


=  inf  ||T-1W  —  MZ\\oo 

Ze-  H-  11  1 


Under  the  assumption  (3.22)  it  follows  that  there  exists  a  unique  Z  £  H°°  satisfying 

l  =  \T~1W-MZ\  a.e. 

From  Proposition  (2)  we  have  that  the  condition 


1  >  sup  Rw-'t(z) 

x€V(F) 


is  necessary  and  sufficient  for 

1  =  zngJT-'W-MFZ\\oa. 

But  we  know  then  that  for  every  e  >  0  we  can  find  Z,  6  H“  such  that 

1  +  e  >  -  M(ju)Zt(ju))\  a.e. 

or 


(1  +  e)\T(ju)\  >  \W(ju)  -  M(jw)Zt(ju>)\  a.e. 


which  implies 


\\T\\>  jniJW-MZ\\. 

Thus,  after  some  more  algebra,  we  obtain  the  desired  conclusion. 


The  assumption  (3.22)  corresponds  to  the  case  in  which  the  norm  of  the  corresponding 
“Hankel  plus  T^eplitz”  operator  [50],  [14]  corresponds  to  an  eigenvalue. 
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3.4  Conclusions 


We  have  established  the  validity  of  “outer  factor  absorption”  as  a  step  in  the  solution  of  H°° 
sensitivity  and  mixed  sensitivity  problems  for  a  wide  class  of  distributed  parameter  plants. 
In  the  process,  we  have  presented  a  constructive  approach  for  “re-inserting”  the  outer  factor 
into,  or  “extracting”  it  from,  a  solution.  Although  we  know  of  no  physical  plants  not  covered 
by  the  class  treated,  it  is  far  from  general:  The  reader  may  wish  to  consider  such  pathological 
cases  as  W  and  F  both  being  discontinuous  at  a  zero  of  F,  F  a  continuous  function  with 
zeros  dense  a  some  point,  or  F  having  zeros  dense  on  some  segment  of  the  imaginary  axis. 
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Chapter  4 

SISO  Mixed  Sensitivity 

4.1  Introduction 

The  objective  in  this  chapter  is  to  solve  as  explicitly  as  possible  the  H°°-infimal  mixed  sen¬ 
sitivity  problem  for  linear  time-invariant  single/input  -  single/output  distributed  parameter 
systems. 

Our  motivation  for  considering  distributed  parameter  H°°  problems  comes  from  work 
intended  to  address  infinite  dimensional  models  (irrational  transfer  function)  of  large  space 
structures.  The  model  motivating  our  departure  from  existing  results  is  the  damped  flexible 
beam  of  Section  2.1. 

In  all  previous  work  ( e.g .  [9] ,[10] ,f  19] ,[51] ,[50] ,[34j ),  the  authors  assume  rational  proper 
weighting  functions,  plant  outer  factors  which  are  invertible,  and  finite  dimensional  instabili¬ 
ties.  Our  assumptions  here  are  more  general  in  each  regard,  although  in  the  mixed  sensitivity 
computations  we  do  not  treat  irrational  weighting  funrtions. 

[50]  assumes  certain  operators  to  be  continuous  on  the  imaginary  axis,  and  this  assump¬ 
tion  was  used  to  compute  the  essential  spectrum  of  an  operator.  In  the  unstable  plant  case 
with  general  inner  factor  this  is  no  longer  true.  In  this  report  we  show  that  the  operator  aris¬ 
ing  from  the  possibly  discontinuous  function  differs  from  continuous  operator  by  a  compact 
operator,  and  therefore  has  the  same  essential  spectrum.  This  idea  is  really  an  application 
of  the  same  technique  which  has  been  used  in  the  distributed-plant  H°°  problem  from  the 
start  [9],  [19],  [10]. 

At  one  point,  computations  of  certain  basis  functions  become  complicated  by  the  pres¬ 
ence  of  irrational  functions,  and  the  dimension  of  a  related  subspace  is  increased  by  the 
introduction  of  an  extra  finite  Blaschke  product.  Below  we  present  calculations  treating  this 
case. 

Finally,  we  present  the  application  of  these  calculations  to  the  model  of  a  damped  flexible 
beam  of  Section  2.1.  This  model  has  an  irrational  and  non-invertible  outer  factor,  and  this 
fact  has  motivated  us  to  develop  the  general  solution  for  that  case. 

We  would  like  to  mention  that  the  paper  [34]  treats  a  similar  problem,  although  only 
rational  outer  factors  are  allowed  there.  Even  in  that  rational  case,  we  believe  that  the 
present  work  indicates  the  connection  with  previous  results  in  a  more  transparent  manner. 
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4.2  The  Mixed  Sensitivity  Problem 

The  general  setup  is  indicated  in  the  block  diagram  below. 


Figure  4.1:  General  Feedback  System 

The  sensitivity  S(s )  is  the  transfer  function  from  the  disturbance  d(s)  to  the  output 
y(s).  The  complementary  sensitivity  1  —  S(s)  is  the  transfer  function  from  the  disturbance 
d'(s)  to  y(s),  which  is  the  same  as  the  transfer  function  from  d(s)  to  y'(s).  Obviously  a 
small  sensitivity  over  a  frequency  range  of  interest  means  good  disturbance  rejection,  and  it 
is  known  that  a  small  complementary  sensitivity  means  good  stability  margin.  [6],  [40]. 

The  H°°  mixed  sensitivity  problem  is  the  following  [32]:  find  a  stabilizing  feedback  com¬ 
pensator  which  minimizes 

sup  (|W1(;w)S(j(i;)|3  +  |Wj(;u;)[l  -  5(ju;)]|3)  . 

U/ 


where  W^s)  and  IVa(s)  are  frequency-dependent  weighting  functions  which  serve  to  empha¬ 
size  or  deemphasize  the  importance  we  attach  to  the  magnitude  of  S  or  (1  -  5)  at  different 
frequencies.  This  problem  is  equivalent  to  the  following:  choose  a  stabilizing  feedback  com¬ 
pensator  C(s )  to  solve 


inf  ||r(j)|U 

C(l) 


(4.1) 


where  we  regard 


T{s) 


Wx(s)5(a) 

W2(s)(l-S{s)) 


as  an  operator  from  the  Hardy  space  H3  to  H3  x  H3. 


4.2.1  Assumptions  about  the  Plant 

We  first  describe  the  structure  of  our  plant  model.  We  assume  that  the  plant  has  trans¬ 
fer  function  P(s),  with  the  following  further  assumptions:  (See  [13]  and  [12]  for  further 
explanation.) 

Assumption  A:  P(s)  has  a  coprime  factorization  P(s)  =  N(s)/D(s)  satisfying 

N(a)X(s)  +  D(s)Y(s)  =  1,  (4.2) 

where  N(s),  D(s),  X(s)  and  Y(s)  are  all  functions  in  H°°. 
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Assumption  B:  Let  N(s)  =  Ni(s)N0(s )  and  D(s)  =  Di(s)D0(s)  be  the  inner-outer  factor¬ 
izations  of  these  functions.  We  assume: 

B.l:  |D0(s)|  is  bounded  away  from  0  at  oo  in  the  right  half  plane  Re(s)  >  0. 

B.2:  Dt  is  a  finite  Blaschke  product. 

B. S:  Let  S0  denote  the  subset  of  the  extended  imaginary  axis  iR  U  {00}  on  which 

WjNoDo  vanishes  essentially.  Then  WiNoD0  is  continuous  at  each  point  of  Sa. 

Assumption  C: 

C. l:  Wx(s)  €  H“  PI  C  is  outer; 

C.2:  Wx(s)W2-x(s)  €  H°°,  W2{s)N(s)X(s)  e  H“; 

C.8:  Wx(s)  and  Wj(j)  are  rational  functions  of  s. 

Remark  10  It  seems  to  have  been  known  that  stabilizabUity  is  equivalent  to  the  existence 
of  a  coprime  factorization  (see,  for  example,  [Zames  and  Francis  P.  590]),  but  the  first 
published  demonstrations  of  this,  for  the  case  of  plants  which  are  ratios  of  H°°  functions, 
appear  to  be  [28]  and  [42].  However  the  plants  considered  by  these  authors  do  not  further 
restrict  the  denominators,  and  so  non-causal  plants  are  allowed.  As  a  consequence  of  this 
setup  ,  non-causal  systems  may  arise  as  “stabilizing”  compensators.  Our  assumptions  B.l 
and  B.2  eliminate  this  possibility. 

Remark  11.  Assumptions  B.l  and  B.2  together  guarantee  that  P(s)  is  analytic  and  bounded 
in  some  right  half  plane.  Thus  P{s)  is  a  shifted  version  of  an  element  of  H°°,  and  so  is  the 
Laplace  transform  of  a  growing  exponential  times  the  (distributional)  impulse  response  of  a 
causal  stable  linear  time  invariant  system.  Assumption  B.2  is  not  necessary  for  this  (we 
could  have  simply  assumed  D(s)  is  bounded  away  from  0  at  00  in  the  right  half  plane),  but 
it  is  satisfied  by  all  physically  motivated  problems  we  have  seen. 

Remark  12  Assumption  B.S  we  use  to  prove  that  we  really  can  “ absorb  N0D0  into  the  free 
parameter,”  as  we  explain  below.  (We  shall  also  need  some  assumptions  on  the  weighting 
functions.) 

Remark  13  We  know  of  no  plants  arising  from  models  of  real  problems  which  violate  any 
of  these  assumptions.  In  the  introduction  we  presented  an  example  consisting  of  a  damped 
flexible  beam  which  has  irrational  Na. 

Remark  14  We  use  Assumption  C.3  to  prove  Lemma  4  *n  Section  4-4-2  and  successive 
results.  We  would  like  to  eliminate  this  assumption  for  the  reason  alluded  to  in  Section  2.2: 
in  order  to  obtain  a  proper  compensator  by  using  an  improper  weighting  function  W-i,  As¬ 
sumption  C.2  would  be  sufficient.  However,  our  proof  of  Lemma  4  requires  that  the  function 
(WiW^/Wy  be  meromorphic  in  the  right  half  plane,  which  will  not  hold  for  some  irrational 
Wj  which  satisfy  C.2.  C.l  and  C.2  imply  Wf1  is  analytic  in  C+  and  outer.) 
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Remark  15  We  allow  W2(s)  to  be  an  improper  function.  In  practice  one  would  pick  W2 
based  upon  one’s  knowledge  of  P(s).  Prom  this  perspective,  we  are  allowing  W2  to  have  poles 
at  the  zeros  of  P{s)  on  the  imaginary  axis,  if  the  designer  wishes.  Actually,  we  not  only  allow 
W2  to  be  improper,  we  recommend  it  for  the  following  reason:  If  W2P  is  not  strictly  proper, 
then  when  we  find  a  solution  to  the  mixed  sensitivity  problem,  the  optimal  complementary 
sensitivity  will  be  strictly  proper,  in  fact  enough  si  to  make  the  resulting  compensator  proper. 
This  can  also  be  accomplished  as  in  [44]  by  considering  W2P~1(  1  —  S)  instead  of  Wj(l  —  S), 
but  the  latter  technique  is  less  general.  (Our  assumptions  encompass  the  case  of  [44]  because 
[44]  assume  W2  and  P  both  have  no  poles  or  zeros  on  the  imaginary  axis.) 


4.2.2  Transformation  to  Standard  Form 

In  this  section  we  review  results  from  [13]  which  show  how  to  find  an  operator 

Tx  .  Ha  -  Ha  x  Ha 

having  structure 


Tx(s)  ± 


G(s)  -  M(s)Z(s) 

F(s) 

such  that  T*(s)T(s)  =  T((s)Ti(s).  These  calculations  may  be  considered  a  clarification  and 
extension  of  the  calculations  in  [8],  where  the  assumptions  about  the  plant  and  weighting 
functions  were  less  general,  and  the  results  slightly  less  explicit. 

Proposition  3  We  can  write 

W}(.)  Wi(») 


with  R  €  H°°  and  outer.  Furthermore,  R  1  G  H°°. 

Proof:  Obviously  is  non-negative.  Thus 

w;{jw)w2{ju)  _  (y  ^ 

VjU <*>)/ 


(4.3) 


—4 

Wdj^Wxiju,)  +  W;{ju>)W7{]u)  \  W:;{ju)W7{y. 


G  L°°(R). 


Now  from  [Hoffman,  p.  53]  (translated  from  HM  of  the  disk  to  the  half-plane)  we  know 
then  that  a  necessary  and  sufficient  condition  to  factor  a  function  f]jw)  €  L°°(R)  with 
/  >  0  a.e.  as  /  =  h*h,  with  h  6  H°°,  is  that  log[/(;w)]/(l  +u;3)  €  L^R).  So  we  check: 


log 


w 


w;0u)w,(iu) 


(jw)Wi(jw)  -i-  w;(ju>)w2(j 


h)I=1o8(i+ 


wdj^Wxjjuy 

w;(jw)w2(ju>)/ 


<  log  (1  +  \wx{ju)wt  1(;^)||^o)  • 


This  gives  (4.3)  by  using  the  result  cited  from  [26].  Taking  an  inner-  outer  factorization, 
it  is  obvious  that  R  can  always  be  chosen  to  be  outer.  That  R  G  H°°  follows  from 


(w* 


0  «L" 
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Remark  16  R(s)  €  H°°  is  outer  since  R  1  6  H°°. 


Remark  17  This  result  does  not  require  Assumption  C.S,  i.e.,  rationality  of  the  weighting 
functions  is  not  required  for  the  factorization  in  (4-8). 


Now  taking 

w(s)iw2(s)irl(s) 

(4.4) 

G(s)^W:(s)W1(s)W(s)~l  -  W(s)N(s)X(s), 

(4-5) 

M(s)  =  N(s)D(s)W2(s)R-\s), 

(4.6) 

and 

F(s)^W1(s)R(s), 

(4.7) 

it  is  easy  to  check  that 

r(s)T(s)  =  (G(s)  -  M(s)Z(s))*(G(s)  -  M{s)Z(s))  +  F*(s)F(s). 


Since  ||T(s)||20  =  ||T’*(a)T’(a)||00  ,  we  can  find  |jT’||0O  by  finding 


G-lfiZ 

F 


Thus 


=  inf 

zeH» 


inf 

zeH« 


G-MZ 
F 

J  HOO 

G  -  NtN0DiD0WiR-xZ 
F 


(4.8) 

(4.9) 


R  is  invertible  in  H°°, 

Z) 


so  we  can  immediately  write  (“absorbing”  R~l  into  the  free  parameter 


=  inf 

ze  H» 


G  -  NiN0DiD0W3Z 
F 


(4.10) 


4.2.3  Absorption  of  the  Outer  Factor 

To  reduce  (4.10)  to  the  desired  form,  we  need  Proposition  2  of  Chapter  3.  The  point  is  that 
we  can  absorb  the  outer  factor  of  NDW3  into  the  “free  parameter”,  just  as  in  the  rational 
[49]  and  rational-with-delay  cases  for  pure  sensitivity  [10],  [18]. 

The  following  is  simply  a  translation  of  Proposition  2  of  Chapter  3  to  the  present  context: 
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Corollary  3  Let  p  be  as  in  (4-8)  above,  and  take 


Mi  = 


G  -  NiDiZ 


F 


Suppose  V(N0D0W2 )  is  a  finite  set  consisting  only  of  ordinary  zeros.  Then  p  =  px  if  and 
only  if 

Pi  ^  sup  U  ess  ran(G*G  +  F*F,  z)  (4.11) 

zeV(NaDaWt) 

is  satisfied.  I 


Remark  18  V{N0D0W2)  u  the  set  S0  defined  in  Assumption  B.8. 

Remark  19  It  is  obviously  the  cast  that  pi  <  p,  so  this  proposition  asserts  necessary  and 
sufficient  conditions  for  p  <  Pi  when  S0  contains  only  ordinary  zeros. 


4.3  Implications  for  Design 


In  the  case  of  stable  plants,  one  can  take  X  =  0  in  (4.2),  and  assuming  the  outer  factor 
can  be  absorbed  into  the  parameter  aZ"  as  described  above,  we  see  that  the  infimal  mixed 
sensitivity  is  independent  of  the  outer  factor  of  the  plant. 

In  greater  generality,  suppose  Di  =  1.  This  means  that  any  unstable  poles  of  the  plant  are 
on  the  imaginary  axis.  Then  using  Corollary  3,  Proposition  2  and  the  fact  that  R~lX  6  H°°, 
we  have 


p  = 


inf 

zeHa 


inf 
ze  H« 


inf 

ze  H« 


inf 
Zg  H“ 


inf 

zeH- 


inf 

ZgH* 


WZWx{W;)-lR*  -  WiR-'NX  -  ndw2z 
WXR 

WiWx{W;)-'R'  -  WjR-'NX  -  NiZ 
WxR 

w*wx{w;)~lR*  -  W2R-'NX  -  w2n0n,z 
WXR 

w;wx{w;  )-1r*  -  w2n(r~1x  -  z) 

WXR 

WfW\(W2)~l R*  -  W2NZ 
WxR 

wfWx{w;ylR*  -  N;Z 
WxR 


Thus  the  infimal  mixed  sensitivity  is  independent  of  the  outer  factor  of  the  plant,  Naj Da. 
An  interesting  application  of  this  observation  in  the  stable  plant  case  is  an  explanation  of 
some  previous  numerical  experiments.  In  the  “Integrated  Structural  Analysis  and  Control” 
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project  at  Aerospace  Corp.  [1],  the  goal  was  to  produce  a  computational  methodology  for 
improving  the  performance  of  large  space  structures  by  simultaneously  designing  the  struc¬ 
ture  and  control  systems.  We  used  a  simple  idealized  satellite  example,  where  we  imposed 
mass  and  cross-sectional  area  constraints  on  structural  members  and  a  complementary  sen¬ 
sitivity  constraint  to  provide  stability  margins,  and  we  applied  a  general-purpose  non-linear 
programming  package  to  minimize  the  weighted  sensitivity  subject  to  these  constraints. 

In  this  case,  we  found  there  to  be  essentially  no  benefit  in  simultaneous  design  over 
separately  computed  structural  and  control  designs.  In  [29]  we  presented  the  hypothesis 
that  this  “separability”  was  due  to  the  absence  of  right  half  plane  zeros  in  the  plant  model. 
The  conjecture  was  heuristically  justified  by  appealing  to  the  well-known  fact  that  for  the 
pure  sensitivity  problem  it  is  the  inner  part  of  the  plant  which  limits  performance.  Thus 
the  idea  was  that  whatever  the  structure  was,  so  long  as  the  plant  had  a  trivial  inner  factor, 
the  limiting  norm  of  the  sensitivity  transfer  function  did  not  depend  upon  the  plant,  and 
amounted  to  whatever  the  complementary  sensitivity  constraint  dictated.  Previously,  it  was 
not  clear  why  the  limit  which  the  complementary  sensitivity  constraint  imposed  should  have 
been  independent  of  the  plant.  The  present  result  provides  the  explanation,  since  for  the 
problem  considered,  the  structural  design  constraints  present  preserved  the  lack  of  inner  and 
unstable  factors. 


4,4  Optimal  H°°  Mixed  Sensitivity 

We  employ  Corollary  3  by  making  the  additional  assumption 

Assumption  D:  n,  >  sup  U  ess  ran (G*G  +  F*F,z )  . 

*es<, 

Accordingly,  we  shall  treat  the  problem 

* -z-SSi.il  (%"*)»-  <*■«> 

where  E,  G,  and  F  are  as  in  (4.3),  (4.5)  and  (4.7),  and  Af  is  given  by 

M  =  NiDi.  (4.13) 

Let 

r = (n+MG*n.M*G  +  n+rF) :  h 2  -  Ha. 

As  pointed  out  in  [44],  /i3  =  ||T||.  This  is  the  starting  point  for  the  development  in  [50]. 

Now  we  parallel  the  developments  in  [50].  The  basic  idea  is  that  if  T  has  an  eigenvalue 
larger  than  the  essential  spectral  radius  of  T,  then  ||T||  is  equal  to  the  largest  eigenvalue  of 
T.  So  one  computes  the  esssential  spectral  radius,  and  searches  for  the  largest  eigenvalue. 
Since 

imi  < 
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we  know  an  upper  bound  for  the  largest  eigenvalue.  Therefore  a  search  can  be  confined  to  a 
known  finite  interval. 

There  are  two  complications  in  applying  [50]  directly.  First,  in  [50]  the  operators  G 
and  F  were  assumed  to  be  continuous  on  the  imaginary  axis.  This  allowed  the  authors  of 
[50]  to  apply  a  theorem  in  [33,  p.  125]  to  the  continuous  (on  j R)  function  G*G  +  F*F  to 
determine  the  essential  spectrum  of  Uk(G*G  +  F*F) |K  (the  compression  of  G*G  +  F*F 
to  K),  where  K  =  H2  0  MHJ  (after  showing  that  T  differs  from  IIk(G;*(?  +  F*F)|k  by  a 
compact  operator).  In  the  present  work  G  is  not  continuous  (in  the  general,  unstable  case), 
due  (at  least)  to  the  term  W(3)N(s)X(s)  in  which  N(s)X(s)  may  be  discontinuous  on  the 
imaginary  axis.  However,  since  we  assume  that  Di  is  a  finite  Blaschke  product  we  are  able 
to  show  below  that  IIk(G!*G+  F’-FJIk  differs  from  the  compression  of  a  continuous  function 
by  a  finite  rank  operator,  and  therefore  the  essential  spectrum  is  the  same  as  that  of  the 
compression  of  the  continuous  function.  Alternatively,  one  can  say  that  T  differs  from  an 
operator  to  which  we  can  directly  apply  the  [50]  computations  by  a  finite  rank  operator,  so 
the  essential  spectra  are  the  same. 

This  is  essentially  the  same  idea  used  in  [50]  (as  mentioned  above)  to  show  that  the 
essential  spectra  of  (the  equivalent  there  of  our)  T  and  IIk(G*G  +  F*.F)|k  are  the  same. 
This,  in  turn,  is  the  same  idea  which  was  used  in  [9,  p.  16]  and  [51]  for  the  pure  sensitivity 
problem. 

The  second  complication  comes  in  computing  the  eigenvalues  of  T,  because  G  can  be  am 
irrational  function.  This  we  resolve  by  means  of  some  additional  computations  as  explained 
below. 

In  [50]  the  idea  is  to  note  that  T  differs  from  the  multiplication  operator  (G*G  +  F*F)\k 
by  a  finite  rank  operator.  Then  the  condition  that  x\  is  an  eigenfunction  with  eigenvalue  A2 
implies  that  the  image  of  i*  under  the  multiplication  operator  A3  —  ( G’G  +  F*F)  lies  in  a 
certain  finite  dimensional  subspace  of  L2.  Computing  residues  of  the  image  of  x\  under  the 
action  of  the  finite  rank  operator  at  certain  points  in  the  complex  plane  allows  one  to  find 
necessary  conditions  for  A2  to  be  an  eigenvalue.  This  computational  procedure  was  also  used 
in  [51]  in  the  pure  sensitivity  case.  Here  we  follow  the  same  approach,  with  differing  details 
to  account  for  both  more  general  assumptions  and  the  particulars  of  the  mixed  sensitivity 
problem. 

4.4.1  Essential  Spectrum  of  T 

By  definition, 

t  =  [n+(A*  +  B*)n_(A  +  b)  +  n+/”F] :  h3  ->  h3  ,  (4.14) 

where 

A  =  N* D* ,  B  =  -  WN0XD *  .  (4.15) 

Note  that  by  Assumption  C,  WtR-1  N0X  6  H°°. 

Now  we  expand  the  first  half  of  the  right-hand  side  of  (4.14)  to  obtain 

n+(A*  +  s*)n_(A  +  b)  =  u+Amn.A  +  n+A*n_s  +  n+s*n_(A  +  b). 
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Flailing  Assumption  B.2  tha.  A  is  <*  finite  eKasenke  product,  we  see  (see  Appendix  A) 
that 

n_S  :  H3  -  Hi  and  II +B*  :  Hi  -+  H2 
are  finite  rank  operators,  so 

n+A*n_S  :  H3  -►  H2 

and 

n+S*II_(A  +  B)  :  H2  -►  H3 

are  finite  rank  operators.  This  means  that  T  is  a  finite  rank  perturbation  of  the  operator 

T0  =  n+A*n_A  +  n+F*F. 


Upon  observing  this  property  of  T,  we  have 
Theorem  3  The  essential  spectrum  ofTis 

<X'(T)  =  {|  W(ju)  |3:  ju>  €  <r.(AA)}  U  {inf  |  F(ju)  |3,sup  |  F(ju)  |3}. 

Proof:  Since  T  is  a  finite  rank  perturbation  of  %,  we  have 

a.(T)  =  a.(T0). 

Furthermore, 

A*  A  +  F*F  =  W’W^W^R^W^W^W;)-1^)*  +  W1R{WlRy 

=  W]WX.  (4.16) 

Wx  €  H°°  is  continuous  on  the  imaginary  axis  since  it  is  rational,  so  the  conclusion  follows 
from  (the  proof  of)  Theorem  3  in  [50].  I 


Remark  20  A  priori,  we  might  wish  to  assume  that  A* A  4-  F*F  €  C  in  order  to  apply  the 
essential  spectral  mapping  theorem  [SS,  p.  125],  as  in  [50].  However,  because  of  (4-16),  we 
only  need  the  continuity  ofWx. 

4.4.2  Eigenvalues  of  T 

Now,  in  parallel  to  the  developments  in  [50],  we  show  how  to  compute  the  eigenvalues  of  T. 
The  idea  in  [50]  translated  to  the  present  context  consists  of  the  following  steps: 

1.  If  #*3  is  an  eigenvalue  of  T,  then  the  eigenspace  KM  of  T  associated  with  the  eigenvalue 
p?  satisfies  KM  C  (4Ha)x,  where  4  is  an  inner  function  to  be  defined  below. 

2.  T  |km=  W*W i  +  AM  where  AM  is  a  finite  rank  operator. 
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(4.17) 


3.  The  condition  for  p  to  oc  an  eigenvalue  is  p2xtl  =  Txp,  or 

(p2  -  W7Wi)*M  =  AMx„. 

4.  Since  AM  has  finite  rank,  we  can  find  a  finite  basis  for  the  range  of  Ap,  which  allows 
us  to  write  a  matrix  representation  of  Ap.  By  identifying  finitely  many  /i-dependent 
points  in  C  at  which  we  know  A^x^,  we  can  find  necessary  and  sufficient  conditions 
for  p  to  be  an  eigenvalue. 

5.  Since  we  are  looking  only  for  eigenvalues  which  may  correspond  to  the  norm  of  T,  we 
need  only  search  for  the  largest  eigenvalue  on  the  interval  (p,(T),  ^  ^  j  ). 

The  following  lemma  is  a  modification  of  [50,  Lemma  1];  the  proof  is  basically  the  same. 
The  only  difference  is  that  here  we  use  the  finite  Blaschke  product  Bw  to  “absorb”  the 
unstable  poles  of  Wf(Wf  )~lR?. 

Remark  21  By  Assumption  C.S,  W\,  Wj  and  (therefore)  R  are  rational.  The  existence  of 
Bw  and  Bu  in  the  following  lemma  depend  upon  this. 

Lemma  4  Let  Bw  be  the  finite  Blaschke  product  which  satisfies  BwWfflVf)"1  R*  €  H°°. 
Also,  let  Bft  be  the  Blaschke  product  whose  zeros  are  those  zeros  of  (p2  —  F‘F )  lying  in 
Re(s)  >  0. 

If  p2  is  an  eigenvalue  of  T,  and  xM  a  corresponding  eigenfunction,  then 

x„€(  B^BwMHDnH3, 


Proof:  Since  p 3x<1  =  Tx„,  we  have 

(, p 2  -  f*f)xh  =  (n+Mcrn.M’c?  -  n  .f*f)xia. 

Multiplying  the  above  equation  by  (p2  —  we  get 

b*b^m*xm  =  (p2  -  F^Fy'BiB^Artn+MG'n.ArG  -  n  _f*f)xm.  (4.is) 

For  any  h  €:  H2, 

<  BwM*(Tl+MG*Tl-M*G  -  H.FmF)xltth  > 

=  <  n+AfCrn_ATG*M,  BwMh  > 

=  <  MCrn.liTGx^  BwMh  > 

=  <  n_Af*Gx„,  M*GBwMh  > 

=  <  U-M*Gxtl,  WxW({W;yxR*Bwh  -  BwWN0XN,h  > 

=  0  . 
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B^M*(U+MG*U.MmG  -  n .F*F)Xti  e  Hi. 

Note  also  that  if  z  is  a  zero  of  (p,3  —  F*F),  so  is  —  z.  By  the  construction  of  B we  see 
that  (p3  —  F’F)-1!?*  is  analytic  on  the  left  half  plane.  Therefore 

(p3  -  F,F)~1B^B^rM*(U+MG*U-M*G  -  n_F’F)xM  €  Hi  , 

and  combining  this  with  (4.18)  we  get 

B^BlyM'x^  6  Hi. 

Therefore 

xM  6  (B^BvirMHl)  n  H3  =  (B^BwMHi)1  . 

I 

The  next  step  is  to  show  that  T  differs  from  the  multiplication  operator  W*Wi  by  a 
finite  rank  operator  when  restricted  to  the  subspace  =  (B^BwM H3)1.  Continuing  the 
parallel  to  [50],  we  present  this  as: 

Lemma  5  T |km  «  a  finite  rank  perturbation  of  the  multiplication  operator  W*Wi  |k„- 

T|Km  =  |Km  - 

where 

A„  A  _n_(AM  +  F*F)  -  A’ll+A  +  II-A'II+A  +  n+A*II_B  +  n+BTI_(A  +  &19) 

and  rank(AM)  =  2[order(Wi)  +  order(Dj)]  +order(S<1)  (see  Appendix  B). 

Proof:  We  first  expand  T : 

r  =  [n+tA’  +  B^n^A  +  sj  +  n+F’F]  i*, 

=  n+A’n.A  +  n+A’n.jB  +  n+fi’n.fA  +  sj  +  n+F’F.  (4.20) 

Now  we  note 

n+A*n_A  =  (A*-n_A*)(A-n+A) 

=  a*a  -  A*n+A  -  ii_a*a  +  n_A*n+A, 

II+F’F  =  F’F  -  II_F*F  (4.21) 

Using  Lemma  4  we  see  that  11+ A  has  finite  rank  on  K^,  and  we  saw  in  the  discussion  be¬ 
fore  Theorem  3  that  U-B  and  11+ B*  have  finite  rank  on  H3  and  H3,  respectively.  Therefore 
we  conclude  that  A-  has  finite  rank  on  KM. 

Combining  (4.20)  with  (4.19),  (4.16)  and  (4.21),  we  can  write 

r  =  ww  +  am.  (4.22) 

To  compute  the  rank  of  AM)  we  find  a  basis  for  the  range  of  AM:  see  Appendix  B.  ■ 
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4.4.3  Computational  Details 

In  this  section,  we  shall  give  further  details  for  computing  the  optimal  performance  /x,  fol¬ 
lowing  the  ideas  of  [50].  We  take  into  account  the  explicit  formulas  above  as  well  as  the 
complications  mentioned  before. 

Applying  (4.22)  to  an  eigenvector  x^,  're  can  write 

(/x2  -  w*w{)x^  =  [n_A*n+A  -  a*u+a  -  n-\v*Wi 

+(A*  -  n_A*)n_B  +  n+£*n_A  +  n +£*n_s]x.  (4.23) 

In  order  to  use  the  condition  (4.17)  to  find  eigenvalues,  we  express  in  terms  cf  a 
basis  for  its  range  when  restricted  to  an  eigenspace  of  T.  From  Lemma  (4)  we  know  that 
xn  G  H2)x,  so  we  need  only  examine  the  image  of  under  AM.  A  direct 

application  of  the  technique  in  [50]  leads  to  expressions  for  a  basis  in  terms  of  rational 
functions  and  M. 

The  calculations  here  are  an  extension  of  those  in  [50],  so  [50]  is  covered  as  a  special  case. 
Note  that  in  the  stable  plant  case  which  [50]  implicitly  treats,  B  =  0  and  G  is  rational. 

In  [50],  the  authors  state  that  rank(A^)  <  2 N,  where  (translating  to  our  notation) 
N  =  order(G)  +  order(F’).  They  go  on  to  treat  in  their  calculations  the  “generic”  case  in 
which  the  poles  and  zeros  of  G,  G *,  F,  F*  and  M  are  all  distinct,  stating  that  in  this  case 
rank(A^)  =  2 N.  The  simplification  (4.16)  shows  that  for  the  mixed  sensitivity  problem 
the  “non-generic”  case  rank(A^)  <  2 N  is  in  fact  generic!  This  observation  shows  that 
the  equation  (18)  in  [50]  gives  insufficiently  many  equations  to  solve  the  problem.  Part  of 
the  solution  is  to  reduce  the  required  number  of  basis  elements,  and  this  is  the  point  of 
Proposition  4  below,  which  shows  that  the  Laurent  series  coefficients  of  at  each  pole  of 
Bw  are  zero. 

In  the  following  calculation,  we  shall  use  f?£(p)  to  denote  the  (— n)th  coefficient  in  the 
Laurent  expansion  of  the  meromorphic  function  /  at  its  pole  p.  We  shall  also  use  Rf(p )  to 
denote  the  residue  ((  — 1)*‘  coefficient  in  the  Laurent  expansion)  of  /  at  p.  Also,  V(f)  and 
2(f)  mean,  respectively,  the  poles  and  zeros  of  /.  (Since  we  assume  /  to  be  meromorphic 
here,  there  is  no  concern  about  non-ordinary  zeros  such  as  may  occur  on  the  boundary  of  a 
region  of  analyticity.) 

Remark  22  Since  W(s)  is  outer  by  construction,  we  have  Z(W*)  C  (right  half  plane). 

Lemma  6  The  poles  of  are  in  the  left  half  plane  and  V(^f)  =  V(^Bw)  =  V(Bw)- 

Proof:  We  know  that  V(^f;)  C  Z(W *)  U  V(W f).  If  p  E  £(W*),  since  Wi  £  H°°,  we 
have  p  £(right  half  plane).  Similarly,  Z(W*)  C  (right  half  plane),  since  W(s)  is  outer  by 
construction.  Furthermore,  since  W*W  =  WfWi  +  WfW2,  if  W*W-L(ju)  >  0,  we  have 
W*W(ju )  >  0,  and  so  W(s)  has  no  zeros  on  the  imaginary  axis.  So  'P(^b)  C  (right  half 
plane).  Noting  that  Bw  is  the  Blaschke  product  such  that  ^ Bw  G  H°°,  and  using  the 
notation  —V  =  {p  :  —p  €  V },  we  see  that  the  poles  of  yfcBw  are  in  the  left  half  plane  and 
V(^BW)  =  -V(%h).  Therefore  V(%)  =  -7>({£)  =  V(^BW)  =  P(BW )■  I 
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Remark  23  We  note  that  {W1/W2)*  Rf  =  ( W\jW )*  and  that  W\jW  is  an  outer  spectral 
factor  of  WiWf+WfWt  •  '^ie  P°^M  ofWi/W  are  the  zeros  of  W,  and  therefore  the  right-half¬ 
plane  poles  of  are  ^*e  zer03  0/  W*,  which  are  the  zeros  of  Bw-  Thus  the  poles 

of  Bw  are  the  poles  of  R.  IfWx  and  Wj  have  no  common  zeros,  zeros(W)  n  zeros(Wi)  =  0. 
In  general,  poles(Sw)  =  zeros(W)  \  (zeros(Wi)  fl  zeros(Wj)). 

Proposition  4  For  any  pole  rj  €  V(Bw),  if  V  &  A)  U  V(Ni),  then 

*%(v)  =  0. 

for  all  m  >  1 . 

Proof:  From  (4.23)  we  can  write 

( n 3  -  WZWxfa  =  -A^n+Ax^  +  A'b.  +  n+R*c_  +  o_, 

where 


a_  =  n.{-W*Wx  -  A'U-B  +  A'U+A)xm 

=  H-BXft 

c_  =  n_(yt  +  S)xM. 

From  Lemma  5  and  the  definition  of  A,  if  rj  is  an  nth  order  pole  of  Bw,  it  also  is  a  pole 
of  A*  of  order  £>n.  Then 

lim(j  -  vYA*h(s)  =  Rf{rj)h(r,)  VheH3_.  (4.24) 

Furthermore, 

lim(a  —  rj)l+1A*b-  =  lim(j  —  rj)Rf’(rj)b _  =  0 

M  ■  *1|  I  '  *t| 

since  is  analytic  in  the  left  half  plane.  We  have 

lim(j-J7)*-n+lII+B*c_  =  0 
*~*v 

since  II+B*c_  can  have  poles  only  at  the  poles  of  D{.  Now 

since  Wx  and  Bw  have  no  poles  in  common,  and  therefore 

lim(j  -  i?)^mA*n+y4xM  =  0  (4.25) 

for  every  m  >  1.  Since  rj  is  an  n**-order  pole  of  Bw,  and  Bw  has  no  poles  in  common  with 
A,  we  have 

n +A*.  =  £  fr^  +  “<*)■  <4  26> 

”  m—1  v  / 
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wWr  a(»)  is  analytic  at  s  =  77.  So  from  (4.24),  (4.25)  and  (4.26)  (taking  m  =  n  in  (4.25)) 
we  have  J 

W*W, 

Rt\l)-^bW(ri)D:{T,)R2>(T,)  =  0. 

Therefore  #*"(77)  =  0.  But  then 

=  ^(«iWWE 

and  =  0.  Similarly  we  see  R^{tj)  =  0  for  each  m  >  1.  I 

By  Lemma  4,  x ^  =  B^BwN^Dih^  for  some  6  H?..  Using  this,  in  Appendix  B  we 
explicitly  compute  the  terms  in  the  right  side  of  the  expansion  (4.23)  for  AM. 


4.5  Example:  The  Damped  Flexible  Beam 

We  next  illustrate  the  technique  developed  in  the  previous  sections  using  the  transfer  function 
model  of  a  damped  flexible  beam  presented  in  Section  2.1. 

for  simplicity,  we  shall  use  the  following  sensitivity  and  complement  ary  sensitivity  weight¬ 
ing  functions,  Wi  =  Wa  —  c(b  +  c),  where  a,  b  >  0.  We  note  that  these  weighting  func¬ 
tions  emphasize  low  and  high  frequency,  respectively,  as  good  design  practice  would  dictate 
(see,  e.g.,  [27],  [40]). 

Taking  (l  and  £a  to  be  the  zeros  of 

w'w  =  +  w;w2 

we  can  write 

W'W  =  ~  ~  ■*)(&  ^  ■s)(^s  ■») 

(a  -  s)(a  +  s) 

We  note  that  jj  +  (&*  —  sa)(oa  —  sa)  has  real  coefficients.  Since  it  cannot  have  pure  imaginary 
roots,  we  see  that  it  has  exactly  two  roots  with  positive  real  parts.  We  can  assume  that 
Re(i  >  0,  Re{ 3  >0.  If  &  is  complex,  then  (2  =  6  • 

Thus  for  this  example  we  can  explicitly  write 


FF 

F 

G 

M 

A 

B 


_ b3  -  s3 _ 

(6  -  J)(6  -  3)(6  +  3)(&  +  *) 

6  +  3 

(J  +  6)(J  +  6) 

_ 1 _ <s  +  &)(*  +  &)  N(s) 

6(j  +  a)(j-&)(j-{a)  (•»  +  <»)  JV(O) 

_ 1 _ 

e(j  +  a)(s  -  6 )(3  -  6) 

c(*  +  6)(*  +  6)  N(a) 

(s  +  a)  '  J\T(0) 


4.5.1  Essential  Spectrum 

NiDi  has  no  essential  singularity  on  j R,  so  by  Theorem  1,  we  find  that  the  essential  spectral 
radius  is  given  by 


P*{T)  =  sup  TOu,)|a 

=  [  nSS*.  if 

l  ^  i3  >  0 

following  a  simple  calculation  to  determine  the  zeros  of  ^  |  F(jui)  |a. 


4.5.2  Eigenvalues 

Assume  n3  >  cr,{T)  —  supw  |/’(jw)|a,  so  fi3  —  F*F(ju )  >  0.  This  means  that  /ia  —  F*F  =  0 
does  not  have  pure  imaginary  solutions,  and  so  it  has  exactly  2  roots  with  positive  real  parts 
which  we  shall  denote  as  3i  and  sj. 

_  (3-31)(3-3a) 

"  (j  +  *»X*  +  *l)’ 

where  we  can  write 


=  MJ*’( fi  -  »X6  -  •>)(«■  +  *)(6  +  >)  -  - »’)• 

In  this  case,  for  any  h  G  Ha  and  any  h-  G  H*,  Tl-Bh  =  0  and  U.+B*h.  =  0,  so  (4.23) 
becomes 

(n3  -  w*wx)x^  =  (n_A*n+A  -  A*n.+A  -  u.w*wx)xu,  (4.27) 

where  A  =  -frr-JV„  - 
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From  Appendix  B,  we  see  that  k  =  l,m  =  2  and  n  =  0,  and  the  corresponding  4>vj  = 
1, 2, 3, 4  are 

1  A*  A*  Am 

3  -  a’  3  +  a’  3  +  3i’  S  +  32 
So  we  can  write  down  the  matrix  A{y)  easily  as 


A(y)  = 


(  _L_ 

(-a  (+a 

_i_  *‘(-0 

-(-a  -C+a 

0  dfc 

c  Js 


*'«) 

C+*t 

*'(-0 

-<+«l 

1 

-»i+£i 

1 

-*i+£j 


£1S1  \ 

<+*j  ' 

-C+*a 
1 

-*j+£i 

_ L_  , 


Then  det[.4(/i)]  =  0  is  equivalent  to 

H-0 


(6  -  0(6  -  0 


det 


1 

1 

1 

a-( 

1 

1 

•J-< 

1 

a-(i 

1 

*1-6 

1 

*3 -{l 
1 

a-<a 

*1-6 

*i-6 

m) 


(6  +  0(6  +  0 


We  solve  the  above  equation  for  the  largest  y2  6  (oa(T),  ~y). 


Remark  24  In  the  stable  plant  case,  if  y2  =  4f,  then  C  =  0,  so  y 


,2-1 


is  a  solution  to 


(4-28).  The  open  loop  case  (i.e.,  with  compensator  C  =  0)  corresponds  to  this  solution. 

Remark  25  In  the  present  case,  Di  —  1.  Since  the  optimal  performance  of  the  weighted 
mixed  sensitivity  doesn't  depend  on  the  outer  part  of  the  denominator  of  the  plant,  so  we  see 
that  (4-28)  doesn’t  depend  on  the  instability  of  the  plant  at  all!  In  other  words,  with  a  stable 
plant  P,  =  Ni,  the  corresponding  equation  (4-28)  be  identical  with  that  of  the  present 

case,  so  the  optimal  performances  for  both  cases  will  be  the  same.  We  know  that  for  the 
stable  plant  case  the  optimal  performance  y,  £  j i ),  so  the  optimal  performance  for 

the  present  unstable  case  y  =  y,  €  (<r«(T),  pr). 

Remark  28  If  N*{—a)  =  0,  repeating  the  above  procedure  we  have 


Bi 


*■“(-*) 


where 


Si  = 


1 

1 

»i-£a 


=  0, 


1 

*3  -< 

1 

•i-(i 


Mfii  _  mi  Miii  _  Mill  , 

*i+Ci  *i+Ca  *a+Ci  *a+6  J 


In  order  to  have  nontrivial  solution,  we  require  that 

rank  (Si)  =  1. 


First 


det 


i  l 

•»-£  «j-£ 

l  l 

«i-£j  *i-6 


=  o  <=>  (6  -  6)(ji  -  si)  =  o 
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We  have  shown  that  the  results  and  techniques  of  [50]  can  be  modified  to  solve  the  H°° 
optimal  mixed  sensitivity  problem  for  a  very  general  SISO  unstable  plants.  Part  of  the 
extension  is  needed  to  treat  the  model  of  a  damped  flexible  beam  [11],  which  has  irrational 
outer  part.  It  is  also  the  case  that  when  N<  is  rational  and  Di  general  inner,  we  can  solve 
the  problem  by  the  same  techniques  given  in  the  previous  sections.  We  have  observed  that 
the  complementary  sensitivity  weighting  function  should  be  taken  to  be  as  improper  as  the 
plant  is  strictly  proper.  Because  of  this,  an  important  area  of  continuing  work  is  the  case  of 
irrational  complementary  sensitivity  weighting  function. 


54 


Chapter  5 

MIMO  Mixed  Sensitivity 


5.1  Introduction 

The  problem  of  computing  the  optimal  performance  for  infinite  dimensional  MIMO  systems 
is  different  from  that  of  the  SISO  case  mainly  because  we  no  longer  enjoy  the  commutativity 
of  scalars.  Under  some  commutation  conditions,  [35]  studied  the  computation  of  the  optimal 
performance  of  a  class  of  infinite  dimensional  MIMO  systems. 

In  this  chapter,  we  consider  H°°-optimal  mixed  sensitivity  design  for  an  idle  speed  control 
model  posed  in  [46].  This  is  a  three- input  four-output  system,  and  the  infinite  dimensional 
parts  of  the  system  are  time  delays  with  different  delay  parameters.  In  [46],  two  first  order 
lags  were  used  to  approximate  these  two  delays.  Following  now-standard  transformation 
of  the  problem  to  an  operator  norm  problem,  in  order  to  solve  this  design  problem  without 
using  rational  approximation  of  any  kind,  first  we  have  to  be  able  to  compute  the  inner-outer 
factorization  explicitly  for  an  irrational  H°°  matrix.  Also  in  this  case,  the  commutation 
conditions  of  [35]  can  not  be  satisfied  due  to  the  two  different  time  delays,  so  that  the 
the  techniques  and  results  developed  there  cannot  be  applied  to  this  design  problem.  We 
resolve  this  difficulty  by  a  more  detailed  characterization  of  the  eigenspace  of  a  certain 
operator.  This  characterization  enables  us  to  compute  the  eigenvalues  and  eigenfunctions  of 
the  operator  explicitly.  Also,  by  using  a  result  proved  in  [50],  we  decompose  the  operator  on 
three  orthogonal  subspaces  of  H3  and  compute  the  essential  spectra  explicitly. 

Although  in  this  chapter  we  restrict  ourselves  to  H°°-optimal  mixed  sensitivity  design  of 
this  particular  model,  the  method  itself  applies  more  generally. 


5.2  Problem  Description 

In  Section  2.4  we  presented  a  model  of  idle  speed  control  for  a  fuel  injected  engine.  This  model 
involved  a  system  having  two  incommensurate  delays,  and  multiple  inputs  and  outputs. 

Our  goal  is  to  design  a  feedback  compensator  controlling  the  valve  input  and  spark 
advance  to  minimize  the  effect  of  the  torque  load  disturbance,  and  at  the  same  time  guarantee 
a  good  stability  margin.  First  we  can  transform  the  problem  into  the  standard  compensation 
configuration  shown  if  Figure  5.1.  d  is  a  disturbance  at  the  plant  output,  u  is  the  control 
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Figure  5.1:  Standard  Compensation  Configuration 


input,  which  consists  of  two  components:  idle  valve  setting  and  ignition  timing  setting,  y 
is  the  speed  output,  e  combines  the  three  signals  (  the  two  control  inputs  u  and  the  speed 
output  y  )  to  be  constrained. 

We  want  to  design  a  compensator  K  =  (Jfej,  Jbj)T  for  the  closed-loop  shown  in  Figure  5.1 

so  that  ||  t vjfid  ||oo  achieves  a  minimum  under  the  constraint  of  the  closed-loop  internal 

Ws  JW 

stability.  We  shall  denote  this  minimum  by  and  call  it  the  optimal  performance,  wi ,  wj 
and  wj  are  weighting  functions  with  ti/i  a  high-pass  filter  and  wj,w3  low-pass  filters.  We 
shall  use  the  following  weighting  functions  as  in  [46]: 


C\S 

a  +  rj  * 


w?  =  c. 


a  +  $ 
a  +  rj  ’ 


u>3  — 


C3 

a  +  r3‘ 


5.3  Standard  Formulation  of  the  Problem 

In  this  section  we  transform  the  problem  into  a  problem  of  computing  an  operator  norm. 
The  basic  ideas  of  the  transformation  procedure  we  use  here  is  by  now  standard:  First 
we  parameterize  ([5],  [48])  all  stabilizing  compensators  and  transform  the  problem  into  an 
infimal  norm  problem  of  the  form: 

*  “  mSL  11  (*c  Z) 

where  X  6  Lfx,  and  C  €  H°*.  In  the  course  of  this  transformation,  an  inner-outer  factor¬ 
ization  of  a  H“  matrix  is  needed.  Although  the  existence  of  the  inner-outer  factorization 
has  been  proven  theoretically  [43],  explicity  writing  down  the  inner-outer  factorization  is  in 
general  difficult  even  in  the  scalar  case.  We  observe  that  for  a  quite  general  case[47],  we  only 
need  an  inner-outer  factorization  of  an  upper  (or  lower)  triangular  or  a  Hermitian  matrix. 
In  the  present  case,  for  example,  what  we  need  is  to  explicitly  calculate  the  inner  and  outer 
factors  of  a  two-by-two  upper  triangular  H*  matrix.  In  this  section  we  also  show  how  to 
get  this  factorization  explicitly. 

The  next  step  is  to  use  a  well-known  result  of  [36]  to  further  transform  this  minimization 
problem  into  an  operator  norm  problem.  Since  the  disturbance  d  is  a  scalar  in  this  case,  we 
shall  see  that  the  operator  we  get  is  a  scalar  operator  on  H3. 
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5.3.1  A  Preliminary  Transformation 

Since  P  is  stable,  using  [21,  pp. 822-824],  with  a  stabilizing  compensator  K  we  have 


fed  —  Pi  1  —  P12QP2I  — 


where  Z  €  H^r 

So  the  problem  becomes 


a0  =  min 
*€HS-Xl 


wipi  0 

0  W2P1 

,^3PlP2  U>2P\P3j 


) 


Z|U 


where  tui  is  a  high-pass  filter,  to3  and  W3  are  low-pass  filters,  Ho  is  the  optimal  performance. 

Without  loss  of  generality  we  can  assume  rw  >  r,  as  in  the  present  case,  and  let  r  =  tv—t,. 
Then 


u0  =  min 


=  mm 
*€Hr« 


^ipi  \ 

0  z\\m. 


W2P1 
W3P1P3  W3P1P2 ) 


/  0 
0 


Let  T\  — 

\W3P1 

B  =  (bij)2X2  satisfying 


and  Tj  =  pi 


We  wish  to  find  A  =  (a,j)3xi  and 


Solving  (5.1)  we  get 


/  A*B  =  2YTa> 

\  =  ra(,). 

4?  =  ifW1  +  W). 

bit  = 


(5.1) 


•)  _  „(•)/.„&?  .  ,J*U‘ 


4r1'=PrK,+«,)-@. 

a,  -  "Wfe 

°1  “  ‘u  ’ 
kn 


611  = 


6u  = 


&*«(s  +  fl,)(ei  -  s)  •  •  ■  (e4  -  s) 
(*«i  +  i)a(*.a  +  s)3(ra  +  a)(r3  +  s)  ’ 

<%P]PvKtKt{9]  -  s2)(ra  -  3) 


(r 3  +  s)(ir#1  +  s)(ir,a  +  s)(*wi  +  j)(t„3  -I-  s)(ej  +  s)  •  •  •  (e4  +  s) 
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£>21  =  0, 

. _ {01  +  •»)(£*!  +•»)•••  (£>8  4~  s) _ 

2  (*«i  +  s)(x,2  4-  a)(Tvi  4-  s)(irv  a  4-  s)(r!  4-  s)(r3  4-  a)(ei  4-  a)  •  •  •  (e4  4-  s)  ’ 

where  ei, ...,  e4  are  the  four  roots  of 

41*  -  -  »J)K\  -  **)(*,*,  -  *)  +  4fi4(*.  -  *X4  -  *)  =  o 

with  positive  real  parts,  and  hi, hg  are  the  eight  roots  of 

(cxv;  -  >’)  -  -  j'x^i  -  ok*!  -  o •  •  ■  («5  -  *’)- 

-4f>]0W.(*.  -  0(4  -  0(4  -  O  =  o 

with  positive  real  parts. 

Further  let  C <•>  =  T f*}  -  A<*>  =  u^Vi**  ~  A<*>. 

It  is  easy  to  check  that  (Ti  —  T3Z)^  =  (A  —  BZ )(*)  4-  C^.  Noting  that 


»-('¥•  ¥  )(lV-  L  ) 


and  €  H°°,  we  have 


where 


~  *Sfe. 11  ( ^  cBZ ) lu  "  zIb'R.  II  ( 14  c  ) lu' 

0  =  (4)=  (^  &  ). 


5.3.2  Inner-Outer  Factorization  of  D 

In  the  following,  we  explain  how  to  get  the  inner-outer  factorization  of  D.  The  key  here 
is  that  we  observe  that  each  element  of  the  outer  factor  for  this  matrix  takes  the  form 
Qi  +  Qae-*^  where  Q\  and  Q3  are  rational. 

Let 


D=fUl  Ua  )(Ul  Ua  )=:UV 

\u3  tt4  )  \v3  VK  ) 


be  the  inner-outer  factorization  of  D.  Since  det{D)  =  d\\  =  det(U)det(V),  and  dn  is  inner, 
we  know  [39]  det(U)  =  adu,det{V )  =  a,  where  |a|  =  1.  Of  course  we  can  set  a  =  1.  So 
from  (5.2)  and  =  I  we  have 


u=(dnvi  dnv ;  \ 

~  V  -U3  vi  /  ' 


v^  +  vP  =  1 
dn{v{v a  4-  V3V4)  =  d\3 
U1V4  —  v3v3  =  1 
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(5.3) 

(5.4) 

(5.5) 


Let  Vi  =  PVi  +  Qv.e  *r ,  t  =  1, . . . ,  4,  r  =  rv  -  rt.  Then  from  (5.3)  we  have 


+  (JW  + /-W  +  =  1  (5.6) 

^vi  Q*\  4"  PviQ«»  =  o  (5.7) 

From  (5.4)  we  have 

p*p*  +  qi,q^  +  +  q:,qv<  =  o  (5.8) 

Q^P*  +  Q'^Pv*  =  0  (5.9) 

<*ii  (P^Qv,  +  p^Qv<)  =  R  (5.io) 


where  R  is  defined  by  dia  =  Re~tT.  From  (5.5)  we  have 


Pv\Qv%  4*  Pv\Qv\  —  PvtQvt  4*  PvtQvt  (5.11) 

QviQv*  =  QvtQvt  (5.12) 

PvlPv.  -  P^P^  =  1  (5.13) 

Now  from  (5.7),  (5.9)  and  (5.13),  we  get  QVl  =<?„,=  0.  So  (5.6),  (5.8)  and  (5.13)  become 

f  ^ + pp  =  i, 

]  PvyP^  4"  Pv»PVt  =  0, 

{  PVlP„4  -  P^Pv,  =  1, 

and  we  get  =  §-,P„,  =  where  A  =  PH  +  PH.  But  \  ^  +  =  PH  +  PH  = 

1,  so  A  =  1.  Therefore  P^  =  P*t,Pv,  =  -P^.  From  (5.10)  and  (5.11)  we  have 

f  P^Qvt  4-  P^Qvt  =  dJjH, 

l  -Pv,Q^  4-PViQv4  =0. 

So  Q„,  =  d{iRPvx,Qvi  =  ^PPv,. 

Now  the  problem  becomes  to  find  Pn  and  P„,  such  that 

PW  +  PH  =  i, 

<  P'H'P**  *=  H<">  C5 

V4  =  P*  4-  d^P^Pe-"  6  H°°, 
va  =  Pj;  +dJxPvlPc-*T  €  H°°. 

Tf  J  _  (^v  *) _ j_e _ n«  /* n«  a*(«l  i 


If  =  irenfe:). we  ^  =  ^;by.  *>;  =  where 

fV.'l  -  t.  .*  -L  f.  _ Lf. 


In  (5.14)  we  have 


/*(•*)  —  £*■*  4"  fr-i-s  -  4"  •  •  •  4-  £o> 
?*(■*)  —  C*a*  4-  Cfc-i3*-1  4-  •  •  •  4-  Co- 

r  U4  = _ nil _ +  gwa«>-~  e  h» 

1  Ua  = _ n*i _ +  /(fi_*i«>-T  e  H°° 
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So  we  require 


nSi)  +  R(S,)e-s^g(St)  =  0 


*  =  1 . h.  Now  we  show  that  Pj-p  +  —  1  is  satisfied  provided  that 

*  -  l&l’  +  I  Cal3  =  1. 

In  fact,  multiply  (5.15)  by  /(5i),  (5.16)  by  —  y(6»)»  and  add  them  up,  we  see 

/W(«i)+SW(«i)-0- 

By  the  definition  of  involution  *,  we  also  see 

/(•)(_5,)  +  5C*)(_5.)  =  o. 

So  =  <^($i  —  s)  •  •  •  (fa  —  s)($i  +  •»)•••  (&k  +  •*)>  and  therefore 

/>(•)  +  p(-)  = _ /H(j)  +  g(,)(j)  _  _  x  =  i 

In  summary,  we  solve  (5.15),  (5.16)  and  (9)  to  get  f(s)  and  g(a),  and 


(5.15) 

(5.16) 


(5.17) 


y  —  (  (A +«)■••(*»+«) 
“  l  ri«)_ 

\  (^i +»)•••(<*+*) 


(«i 

OiRj^MihZ^ 

(«i -•)-(«»-») 


TJ  —  (  (Ji+*)"(L+«)  (4+*)"(L+*)  ^  ^  f  /  ( J)  5  (5)  ^ 

wsfSfcJ  «,+.)-»+*) U»(.)  /WJ' 

Since  (A  —  DZ)W  -f  =  (CfM  —  VZ)^  +  C ^  and  noting  that  V  is  invertible  in  H°°, 

we  have 

=  11  (^’V2)11- 

where  C  €  H°°. 

5.3.3  Transformation  to  an  Operator  Norm 

In  the  following,  the  above  infimal  norm  problem  is  transformed  into  an  operator  norm 
problem  in  a  standard  way  using  a  well-known  result  of  [36]. 

Let  T  :  H2  -  )  be  defined  as 

r/  -  (7/)  =  *  h3. 
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Here  (H^xl)_  =  L|xX  ©  Hfxl.  Then  we  have  [36,  Theorem  1] 

m.  =  \yr  Iko  ll, 

where  K0  = 

Let  C  =  nKoT  :  Ha  -  K0,  then  f  =  ||T*  |K#  ||a  =  ||CC*||  =  ||C*||*  =  ||C||a  =  ||C*C||. 
Since  C  =  (n(H!,^-C/M)  :  H>  -  K„,C-  =  ((n(H,  ^VAf.C-)  :  K„  -  H». 

For  any  x  €  (H|xl)_,y  €  Ha, 

<  (n(H| xl)-U*AYx,y  >=<  s,II(H*xl)_C/'My  >= 

=<  xtU*Ay  >=<  A'Ux,y  >=<  IIh *A*Ux,y  > 

and  for  any  y,  z  €  Ha 

<  C *y,z  >=<  IlHaC,*y,z  > 

so  =  IIh^-A*!/  and  C*  =  H.J&C*  and  thus 

=  nH»A*im(H2xi)_tfM  +  IlH>(7,C7  (5.18) 

We  conclude  that 

nl  =  sup  spec(C’C)  (5.19) 

Remark:  By  [30,  Theorem  2],  we  know  if  p0  >  |jC|joo  then 

nl  =  sup  spec(Hl.AHu'A  +  7c*c)- 

Here  'Hu- a  and  Tc-c  denote  Hankel  and  Toeplitz  operators.  It  is  easy  to  see  that  'H'U.A  = 
so  (5.19)  follows. 


5.4  Computation  of  the  Optimal  Performance  n0 

From  the  previous  section  we  know  that  the  optimal  performance  fi0  is  equal  to  the  square 
root  of  the  norm  of  the  operator  C*C.  In  this  section  we  shall  show  how  to  compute  the 
discrete  and  essential  spectra  of  this  operator.  From  (5.19),  these  two  computations  allow 
us  to  find  Ho  (as  in  [51]  and  [50]).  Note  that  all  the  operations  in  theorems  1  and  2  are 
explicitly  computable  using  the  expressions  in  section  4. 
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5.4.1  Discrete  Spectrum  of  C M 

/  g-***  o  \ 

Let  E  =  ^  q  J,  then  A  =  E*Ait  where  Ax  is  rational.  Since  in  general  t,  /  r„,  the 

multiplication  of  .£*  and  Ax  is  not  commutative,  so  the  techniques  used  in  scalar  case  [50] 
to  cancel  the  infinite  dimentional  parts  of  the  plant  no  longer  work  for  us.  [35]  assumes  the 
following  commutation  condition: 

e*a,  =  n  hM;, 

where  fl*  is  rational  and  Af*  is  inner  (possibly  infinite  dimensional),  so  that  the  cancellation 
of  the  infinite  dimensional  part  works  as  it  does  for  the  scalar  case.  It  is  easy  to  see  that  this 
commutation  condition  can  not  be  satisfied  in  the  present  case.  The  difficulty  of  the  present 
problem,  of  course,  is  mainly  due  to  the  noncommutativity  of  the  matrix  multiplication. 

Since  the  multiplication  of  E *  and  A\  is  not  commutative,  for  an  eigenvector  xM  of 
associated  with  the  eigenvalue  p2,  our  idea  here  is  to  characterize  AxXM  instead  of  xM  as 
[50]  does.  The  proof  of  the  following  lemma  is  basically  the  same  as  Lemma  1  in  [50]. 
The  rationality  of  Ax  and  U  allows  us  just  to  introduce  an  extra  finite  Blaschke  product  to 
“absorb”  the  unstable  poles  of  U* A\A\. 

Lemma  7  If  p3  is  an  eigenvalue  of  and  xM  is  an  associated  eigenfunction,  then 

AiXft  €  5j*5adt(H2Xi)_ , 

where 

bp  is  a  finite  Blaschke  product  such  that  5M[(pJ  —  C^)-1]*  £  H03; 
ba  is  a  finite  Blaschke  product  such  that  baU*  A\A\  €  H^2; 
dt  is  an  inner  function  such  that  d^E*  £  H^,2. 

Froof:  Since  /i2xM  =  we  have 

(/*’  -  c<->k  =  -  nH. 

Multiply  the  above  equation  by  (p3  -  C^*^)_16*6*dJAx,  we  get 

WWA^  =  (p3  -  c<*>)-,6X<CAx(nH»A^n(H|xl)_^A  -  nHlc^)xM  (5.20) 

For  any  h  £  H|xl, 

<  6*dJAx(nH*A*C/’II(a»)|j_C/*A  -  nH*_C7^^)xM,  h  >= 

=  <  {UwA9UTl{B*Ki)_irA  -  YllllC^))xtl,bad.A\h  >  = 

=<  A'Un{Hlxi)_U'Ax„bad.Alh  >=<  U{Wtxj_U*Ax„U'AbadtA\h  >= 

=<  Tl(Blxi)U-Axtl,bad.U-E-AlA\h  >=  0. 
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This  means  iadJAi^ipA'C/TIfHi^^l/'A  —  IIh*  C^)xm  G  (Hjx1)--  By  the  definition  of 
and  (5.20), 

Kb'.<r.AlX>e(  HS«)_, 

therefore  the  lemma.  ■ 


Remark,  b M,  ba  and  d«  can  be  chosen  as  follows: 

(1)  6m  can  be  chosen  as  the  Blaschke  product  whose  zeros  are  those  zeros  of  fj,2  — 
lying  in  iJe(s)  >  0. 

(2)  ba  can  be  chosen  as  the  Blaschke  product  whose  zeros  are  those  poles  of  U*AiA[  lying 
in  f2e(s)  >  0. 

(3)  d,  can  be  chosen  as  detE  or  e-*m“(r*>T0  = 

Note  that  we  can  rewrite  (5.18)  as 

C(,)  =  U^A9U(I -UKlJU*A  +  UBiC^ 

=  nm7(AW  +  cW)-nwA*unB,xu9A 
=  nH,(AM  +  cW)-{A'u -naiA'U)nH,xu*A  (5.21) 

It  is  not  difficult  to  see  that  Lemma  7  does  not  provide  enough  information  about  the 
eigenspace  of  C^m>  for  us  to  conclude  that  C^  —  (A^  +  C^)  is  of  finite  rank  on  the  eigenspace 
because  e_rr,H_  C  e_rr*H_.  The  following  lemma  provides  us  further  information  on  the 
eigenspace  so  that  we  can  show  that  C ^  —  ( A ^  4-  C^)  is  indeed  of  finite  rank  on  the 
eigenspace. 

Lemma  8  If  xM  is  an  eigenfunction  of  then 

*„  €  F(Ha)  +  flxe-^'H*, 
where  F  is  of  finite  rank,  R\  is  rational. 

Proof:  Since  from  (5.21)  we  have 

0  =  —  n2Xn 

=  (A<*>  +  C(,)  -  mj)*„  -  Bhi  (A<*>  +  C(,))xM  -  {A*U  -  nIP_A*U)Uv»nlU9Axl, 

=  ( A <*>  +  C(,)  -  /iJ)xM  -  nH»_(A(,)  +  C(,))xM  -  A*f/nH|xlC/*AxM  +  nw_A*UTlmlxU' Ax^ 

(5.22) 

Noting  that  both  IIh*  (A^  +  C^)  |h»  and  nH>  A*U  |jp  are  of  finite  rank,  we  rewrite 
(5.22)  as 

0  =  ( A .<*>  +  C(,)  -  mJ)xm  ~  A*t/nH|)<ll/*Ax^  +  Fi  (5.23) 

where  F\  is  of  finite  rank. 

In  Lemma  7,  if  we  choose  d,  =  e-rr* ,  we  have 


n bUu-Axu  =  nB,'U'E'A,x„  =  ne,nU-E-b^d,  (£- )  = 
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=  nH uM.tr  J)  (£)  =  ('T^-)  = 


=  n H.„  W  (nH- -"**-)  +  ne.„4At/-n(HUi)_  (*'^-)  = 


=  w  (“"V*1*)  -  n(H;„,)_W  (nH’7k-)  + 

■er  =  r„  -r„  (**')  €  (H5„i)-. 

Since  both  II(H lxl)-Kb<*U*  |h»xi  and  l(H|xl)_  Lave  finite  rank,  we  have 

nsj.^Mi,  =  W  (nH’ej'r'‘‘*)  +  F,  (5. 


Here  r  =  r. 


=  i„60£/-  (nH’ej"'‘1*)  +  F, 


(5.24) 


where  Fj  is  of  finite  rank. 

Combining  (5.23)  and  (5.24)  we  have 


0  =  (44(,)  +  C’(*)-M3)*M-M»^‘(nH,C0"Tfcl") +^3 

=  (A<*>  +  C<*>  -  M2)*m  -  M.  ^  (  nH,c0'T/l1-  )  +  F3 

=  (A(#)  +  C<*>  -  n2)xv  -  bjaa{  1e-*T*nH»e‘'r/ii-  +  F3  (5.25) 

where  i4J  =  (oJlt  oJa),  F3~F\  —  A*UF2  is  of  finite  rank. 

Now  from  Lemma  7  we  have  hx_  =  b*  6*  dJanxM,  so  (5.25)  becomes 

0  =  (A<*>  +  C<*>  -  ^2)xM  -  bM6aoJle_rr*IlH>e_*Tb*b*d*anXM  +  F3 
=  (A<*>  +  C(,)  -  /i2)x„  -  bMbaoJ1e-rr'nH*b*b*onerr*xM  +  F3 
=  (A<*>  +  C(,)  -  /i2)xM  -  +  6Mbaa1*1e"T-nHi^b;axle"'xM  +  F3  (5.26) 

x4  =  (m3  +  a};1  -  aw  -  c«)-1f3+ 

+(^3  +  aft  -  -  £<*>)- 

This  completes  the  proof.  I 

Now  from  the  second  equality  of  (5.26)  and  Lemma  8  we  get 

+  bM a* j e ~ II h j a 1 1>( /i2  +  aft  -  A<*}  -  CW)-lnH»  ^aaxxe^x,,  +  F3 

(5.27) 
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Noting  that  -  A ^  -  C^)~l  Ih*  is  °f  finite  rank,  we  see  that  actually  the 

right  hand  side  of  (5.27)  is  of  finite  rank.  So  we  have 

Theorem  4  Let  E  be  the  eigenspace  of  an  eigenvalue  p2  of  C^*\  then 

A  =  Ch  -  (A<*>  +  <?(*>) 

is  of  finite  rank  on  E. 

Proof:  From  the  above  derivation  we  have  that  —  C ^  equals  the  right  hand  side 

of  (19),  which  is  of  finite  rank  on  E.  I 

Now  if  p2  is  an  eigenvalue  of  C ^  and  xM  a  corresponding  eigenfunction,  then 

C^Xft  =  p2x„. 

From  Theorem  1,  we  have 

[p3  -  (A<*>  +  C(#))]xM  =  Ax„. 

Since  the  operator  A  is  of  finite  rank,  we  can  calculate  the  eigenvalue  p2  and  the  corre¬ 
sponding  eigenfunctions  as  in  [50]  and  [14]. 


5.4.2  Essential  Spectrum  of  C 'W 

In  this  section,  we  shall  show  how  to  compute  the  essential  spectrum  of  C^.  We  shall 
consider  the  following  more  general  form  of  A: 


fm{  0  \ 

\  0  m;  J 


Aj 


where  mj  and  m2  are  arbitrary  inner  functions  and  m i  divides  m2  in  H°°.  In  fact,  we  have 


Theorem  5  The  essential  spectrum  of  C ^  is 

<Ta(C(,))  =  { AP(jv)  +  :  ju  G  o,(mj )}  U  { inf  [a^O'w)  +  C{*\]w)\ , 

sup  [aSJO'.w)  +  <7w(jv)]}  U  {inf  Cw(j,w), sup  C(,)(ju/)}  . 

U*  W  tat 

Prior  to  proving  this  theorem,  we  give  some  preliminary  results.  We  make  the  following 
decompositions 


HJ  =  Hi  ©  H2, 

Hi  =  Hu  ©  Hu, 

where  H;  =  maH3,  Hu  =  (miHa)x  =  (miHl)  D  H3.  For  any  operators  X,  Y  in  a  Hilbert 
space,  by  X  ~  Y  we  mean  that  X  —  Y  is  compact.  First,  we  shall  prove  some  projection 
properties. 
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Proposition  5  For  any  rational  R  6  L°°,  any  inner  functions  m,  m3,  and  m  divides  m3  in 
H°°,  m3  divides  m3  in  H°° 

(a) 

(b)  11^1713  R(h> 

are  both  finite  rank  operators. 

Proof:  Let  m3  =  mmo,  since  IIh»-S|ep  is  a  finite  rank  operator,  we  have 

=  nH,7n2(IlHJ  +  ~  IltijmjIlHs =  0  . 

This  proves  (a). 

Similarly,  let  m3  =  Since  IIgj  i2|g*  is  of  finite  rank,  we  have 

UHlm3R\w  —  ^Bj  1712(113*  +  IlH^)f2m4|H*  ~  113,1713  nBiiZr^lB*  =  0  • 

This  proves  (b).  I 


Proposition  0  For  any  rational  R  €  L°°,  the  following  operators 

nsjAig^  iib^Ih^,  nHlJi2|BU 

are  of  finite  rank. 

Proof*  Since  H3  —  m;Ha  by  Proposition  5(b) 

nauftlH,  =  113^171213, 

is  of  finite  rank. 

Let  H  =  mx  H3.  Similarly  we  know  that 

nBu.R|B  =  IlBuf2m1|B, 

is  of  finite  rank.  Since  H13  =  H  fl  Hi  C  H,  so  IIb,,/?^,  is  of  finite  rank.  By  observing  that 

nHj-RlHt  =  (nHli2-|H,)* ,  nHl,iZ|Btl  =  (iib^’Ib,,)* 

we  also  see  that  and  nB„i2|Hii  are  of  finite  rank.  I 


The  following  lemma  is  a  result  in  the  proof  of  Theorem  1  of  [50]. 

Lemma  9  ([50])  For  a  self-adjoint  operator  X  on  a  Hilbert  space  H,  let  S  be  a  subspace  of 

H  If 

XUs  -  IIsA-IIs 


is  compact,  then 


<7e(nsXns)  =  at{X)  . 
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Proof  of  Theorem  2: 

Let 

Ai  4  Y[wA'[mw  s  U*A, 

'**1x1'- 

Aa  4  n^A-Upa^U-ATl^, 

A3  =  (^)anmiriHia- 

Then 

Ai  =  +  nHj)  ~  nHivr^n(H«xl)_£/MnHl  = 

=  (Hh!  +  nn,)>i*i7n(H*)<l)_£/'>inHl. 

By  Proposition  5(a), 

nHtA*tm(H3xt)_tfM~o 

and  therefore 


Ai  ~  UBiA-imiElxi)_U-ATlni 

=  nHlAtA!nHl  -  a2 

We  now  compute  Aj.  First  we  have 


(5.28) 


A'U  =  (t*u,  u31)aj1mi  +  (uu, tinja^ma  . 

It  is  easy  to  see  that  IIh*^  is  of  finite  rank  since  is  rational.  Also 

by  Proposition  5(b),  nHl(«ij,ua3)<x*2ma|H|xl  ~  0-  So  we  have 


Aa  ~  nHl(un,itai)aJ1m1IlH3)<l  ^“JoumJIlH, 

=  (IlHu  +  nHwXuu.UaO^mjIlHt^^^anm^nHu  +  nHlJ) 

~  A3  (5.29) 

The  last  equivalence  ~  is  because  that  (^l)anmJIlHu  is  of  finite  rank  and  also  by 

Proposition  5(b)  nHu(«ii,U3i)oJ1mi|H|)|i  ~  0. 

Further,  since  Hia  C  miH3,  we  see  that  II(H*xl)  (^l)anminHia  is  of  finite  rank,  so 

Aa  ~  nHl,(«u.  u2i)o*j  OuIIhij  -  nHl>(«u,  uai)a;im1n(H3>|l)_  ^  j  anm*nHlJ 

"  nHwOliOnllHu  (5.30) 
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here  y(,un  +  =  1  since  U*U  =  I. 

Combining  (5.18),  (5.28),  (5.29)  and  (5.30)  we  have 

Cw  =  -  nH,.o;1nH„  +  nH.c<-> 

By  Proposition  6,  we  further  have 

e<->  =  +  cw)nH„  +  nHll(^S-)  +  <?<•>  -  +  nH,c(->nH, 

and  therefore 

*.(£<•>)  =  +  cH)nHll)u  a.(nHli(a(:2]  +  c<*>) nHlJ)  u  a.(nH,c'(*)nHl)  . 

Noting  that 

(at?  +  cW)nH„  -  nH„(at?  +  cw)nHl,  =  nas>i?  +  c<->)na„+ 

+(nH.  -  nH„)(ot?  +  cf>)nH„  ~  (nH.  -  ne, ,)(«[?  +  c<->)nHl,  = 

=  nHu(at?  +  c<-))nHl,  +  nHl(at?  +  cM)na„  ~  nH,(at?  +  c«)nB„  . 

By  Proposition  6,  nBj(at?  +  C^)IIBl  ~  0  sad  since  Hia  C  H, ,  so 

nH,(°S?  +  c'->)nHl,  ~  o 

and  therefore 

(<•!?  +  cW)nHl.  -  nH„(aS?  +  cW)nB„  ~  o . 

So  by  Lemma  9,  e.(IIBls(at?  +  C(*')IIBll)  =  a.(at?  +  C(*i). 

Similarly,  by  Proposition  6, 

c^fin,  -  nH3c(*>nH,  = 

nHiC(*)nHj  +  (nH>  -  nHa)CHnHl  ~  nHlc(,)nHj  ~  o . 

So  by  Lemma  9,  ^(rinjC^IlH,)  = 

Note  for  any  rational  R  €  L00,  iZIIn*  =  IlHJ^nH*  +  ~  So 

cre(R)  =  cre(IlH*-ft)-  By  [38,  Theorem  6.2  (i),  p.  55]  we  have 

cre(4 V  +  ^(,))  =  { inffa^O’w)  +  C^\jui%  sup^j^ju;)  +  C{m)(ju;}}  , 

w  <y 

<7t(C^)  =  {inf C^(jw),  supC^(ju>)}  . 

w 

By  [33,  Corollary  1,  p.  125]  we  have 

^(IIh,^  +  C,(*))nHu)  =  {^(iw)  +  C^\ju)  :  ju>  6  cTe(mj)}  . 

I 
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5.5  Optimal  Compensators 

After  obtaining  the  optimal  performance  fi0,  one  can  compute  the  suboptimal  compensators 
using  the  rational  approximation  techniques.  In  the  multivariate  case,  even  for  rational 
plants,  there  is  more  than  one  optimal  compensator  in  general  [22].  However,  in  the  present 
case,  we  observe  that  our  problem  is  a  “vector”  problem  rather  than  a  general  “matrix” 
problem,  and  the  operator  we  get  in  the  previous  section  is  a  scalar  operator.  This  makes 
it  possible  for  us  to  find  an  eigenfunction  corresponding  to  the  largest  eigenvalue  fi0  of  the 
operator,  and  then  we  can  find  the  optimal  compensator  as  in  the  scalar  case. 

Assume  that  \ia  >  ||C'||00.  Using  the  technique  described  in  [30],  we  can  transfer  the 
original  mixed  sensitivity  problem  into  the  following  pure  sensitivity  problem  (with  the 
optimal  performance  yia  =  1): 

min  ||/f  —  Z||oo  =  1, 

where  =  n\  -  M,  €  H00,  H  =  WAM~l  6  L°°. 

This  is  a  infinite  dimensional  pure  sensitivity  problem  with  two  different  time  delays.  The 
method  developed  in  the  previous  section  can  be  used  to  find  the  maximal  eigenfunctions 
(corresponding  to  the  largest  eigenvalue  =  1)  of  the  corresponding  Hankel  operator  and 
the  minimal  sensitivity  can  be  found  as  described  in  [51]. 


5.6  Conclusion 

In  this  chapter  we  have  shown  how  to  compute  the  infimal  mixed  sensitivity  for  the  full 
infinite  dimensional  model  of  multivariable  pi  nt  considered  in  [46],  without  resorting  to 
finite  dimensional  approximation.  We  axe  presently  preparing  to  compute  numerical  values 
for  comparison  with  the  results  in  [46].  The  procedure  presented  here  is  readily  generalizable 
to  a  class  of  MIMO  plants  [47], 
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Chapter  6 

Numerical  Inner/Outer  Factorization 


6.1  Introduction 

The  motivation  for  the  developments  in  this  chapter  is  interest  in  explicitly  computing  solu¬ 
tions  to  if  “"-optimal  control  problems  for  distributed  parameter  systems,  such  as  the  damped 
flexible  beam  in  Section  2.1.  In  this  case  one  has  an  irrational  transfer  function 

P(s)  =  N(s)/D(*)  (6.1) 

where  N  and  D  are  both  functions  in  H°° .  The  solution  to  the  mixed  sensitivity  problem 
[50] ,[15]  involves  finding  the  largest  solution  fio  on  an  interval  [o,  6]  C  R  to  an  equation  of 
the  form 

MOCOO)S.OO  =  tt(-i<00W(0  (6.2) 

where  Ni  is  the  inner  factor  of  N  (see  Example  1  below  for  further  details),  and  (,  R\,  and 
Rj  axe  certain  functions  of  /i.  (/io  is  the  largest  eigenvalue  of  an  operator  which  is  related 
to  the  mixed  sensitivity  operator.  See  [50].)  If  we  search  for  /xo  by  numerically  computing 
solutions  to  this  equation,  computation  of  the  infimal  values  for  H' °° -control  criteria  will 
involve  the  numerical  evaluation  of  the  term  N%(jw)  for  different  values  of  uj. 

For  general  transfer  functions,  one  will  not  be  able  to  explicitly  write  down  an  inner-outer 
factorization,  and  so  the  evaluation  of  the  term  Ni(ju)  is  not  simply  a  matter  of  substituting 
a  number  into  a  formula.  In  this  chapter  we  present  a  numerical  technique  for  this  evaluation 
without  explicitly  finding  the  inner  factor. 

We  present  three  examples  of  irrational  inner  factors  on  which  we  will  test  our  numerical 
technique.  In  Example  1  below,  the  inner  factor  of  the  transfer  function  is  am  infinite 
Blaschke  product.  In  this  case  we  can  explicitly  find  the  right  half  plane  zeros,  but  there 
are  infinitely  many.  For  Example  2,  the  transfer  function  has  a  singular  inner  factor,  yet  an 
explicit  factorization  and  formula  is  available.  In  Example  3,  the  inner  factor  appears  to  be 
an  infinite  Blaschke  product,  but  the  alternative  of  taking  a  continued  product  of  separately 
computed  zeros  is  not  attractive  since  the  implicit  nature  of  the  expression  for  the  transfer 
function  and  lack  of  further  information  on  the  location  of  the  right  half  plane  zeros  make 
their  computation  difficult. 
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6.2  Mathematical  Framework 


The  outer  factor  of  an  H°°  function  is  given  in  the  open  right  half  plane  by  [26,  p.  133] 

,  ,  .  f  1  ,  ......  <3  + 1  dt 

On  the  imaginary  axis  the  boundary  value  function  of  f0  exists  almost  everywhere  and 
is  equal  to 

F 1  ,  ,  ,, .  v,  ts  +  t  dt 

-/  log|/(.OIT 


lim  exp 
o 


t  t  + 


+  ts  1  +  t2 


(u-il’  +  o-1  ) J  *  A  "l+t’j 


Taking  s  =  <r  4-  tu>,  with  ar,  u  6  R,  we  have 

=  Um«p {i /_”  [t^^)  - 

We  know  the  limit  of  the  real  part  of  the  integral,  since 

17-MI  =  I/Mi  a.e. 

The  limit  of  the  imaginary  part  of  the  integral  in  the  last  expression  is  defined  as  an  integral 
on  the  imaginary  axis  in  the  sense  of  a  Cauchy  principal  value,  which  we  can  simplify  as 
follows: 

=  (* + (r«_- o* +»,))  iogi/(i,)|Tf?} 


uit  *+■  1 


=  --4‘ 

TT  i-oo  (1  -f  t2)(u  —  t) 

■-'AC 


log  |/(*0I* 

/(**) 


(1  +  t3)(<jj  -  t ) 
Considering  the  last  term,  we  have: 
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Now  we  note  that 


U)t  4-  1 
u>  —  t 


log 


m 


/M 


has  a  removable  singularity  at  t  =  u>,  and  that  since 


log  !/(»*)! 

1  +  fa 


(because  f(s)  €  [26,  p.  51]), 


u/t  +  1 


Thus 


(l  +  t*)(cj-t) 
tdt  +  1 


log 


€  L1 


m 


r. 


oo  (1  -f-  fa)(w  —  t) 


/(*") 

/(*) 


€  I1- 


log 


/M 


is  finite  as  an  ordinary  (Lebesgue)  integral.  Therefore,  we  can  compute 


arg  (/>M)  ="J_ 


ut  -f-  1 


(i  +  i>)(«-<)log 

by  evaluating  the  integral  numerically  as  an  ordinary  integral 


m 


dt 


6.3  Discussion  of  Numerical  Techniques 

Four  problems  arise  in  the  numerical  computation: 

1.  the  integration  has  infinite  range, 

2.  the  integrand  has  a  removable  singularities  at  t  =  w,  and  possibly  other  points  which 
may  be  due  to  the  formula  for  the  transfer  function, 

3.  the  integrand  will  have  integrable  singularities  at  the  (possibly  non-rational)  zeros  of 
/(«*;),  and 

4.  the  presence  of  the  logarithm  in  the  kernel  may  lead  to  intrinsic  computer  word-length 
problems. 

The  first  three  are  standard  problems  in  numerical  quadrature,  and  our  approach  is  to 
assemble  standard  solutions,  as  we  discuss  below.  The  fourth  problem  is  more  serious,  as  we 
discuss  in  some  detail  below. 

To  complete  the  numerical  computation  of  the  inner  factor  of  /,  After  finding  the  nu¬ 
merical  value  of  arg(  we  can  simply  compute 

_ /M _ 

I/Ml  exp(»  '  arg(/0(w)) 


/.(«*>)  = 
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We  now  discuss  the  numerical  problems  mentioned  above,  and  describe  our  approach  in 
their  solution.  The  basic  reference  we  have  used  for  a  discussion  of  the  issues  and  techniques 
in  quadrature  is  [37].  (For  the  actual  computer  programs  we  obtained  updated  versions  of 
the  programs  in  [37]  via  electronic  mail  (“netlib”)  from  Oak  Ridge  National  Laboratories.) 

1.  Infinite  range  of  integration.  There  are  basically  two  standard  approaches:  either 
(a)  transform  the  range  to  a  finite  interval,  or  else  (6)  approximate  the  infinite  integral  by 
integration  over  a  finite  range.  Method  (6)  is  recommended  in  [37,  p.  80]  for  the  case  in 
which  the  integrand  decays  “rapidly”  to  zero  at  oo,  and  this  is  the  method  we  have  adopted. 
We  remark  that  we  have  also  tried  method  (a),  and  this  approach  can  generate  an  additional 
problem:  if  the  integrand  does  not  decay  sufficiently  quickly  at  oo,  then  when  transformed  to 
the  finite  interval  the  integrand  may  be  unbounded  at  the  image  of  oo,  leading  to  inefficient 
stepsize  selection  for  an  adaptive  quadrature  program.  We  encountered  this  problem  in 
Example  1  discussed  below,  and  the  difficulty  is  compounded  by  the  wordlength  problems 
in  evaluating  the  integrand  in  the  vicinity  of  oo,  which  we  discuss  further  below. 

2.  Removable  singularities.  For  finite  removable  singularities  this  is  easily  handled  by 
using  a  polynomial  approximation  in  a  neighborhood  of  the  singularity.  In  our  Example 
1,  we  have  a  removable  singularity  at  0  due  to  the  expression  for  the  transfer  function. 
Here  we  used  a  linear  approximation  to  the  integrand  over  a  neighborhood  which  yields  a 
variation  in  the  function  two  orders  of  magnitude  less  than  the  estimated  absolute  accuracy 
requested  of  the  adaptive  quadrature  program.  For  the  singularity  at  u ;  one  could  use  a 
similar  approach,  although  in  our  tests  we  simply  avoided  values  of  u  which  were  likely 
abscissae  in  the  quadrature  (by  taking  logarithmically  spaced  points  at  which  to  evaluate 
the  outer  factor). 

3.  Integrable  singularities.  Since  the  integrand  is  in  L1,  of  course  all  singularities  on 
the  imaginary  axis  will  be  integrable.  So  far  we  have  only  considered  an  example  with  an 
irrational  zero  at  oo.  We  expect  most  finite  zeros  to  be  easier  to  treat,  for  the  following  reason: 
A  finite  zero  of  the  transfer  function  will  lead  to  an  unbounded  but  integrable  integrand  at 
the  zero,  due  to  the  log  function.  Such  singularities  are  treated  in  [37].  The  idea  is  to  break 
up  the  interval  of  integration  into  subintervals,  with  the  origin  as  an  endpoint,  and  then  use 
a  specialized  routine. 

In  the  case  of  the  zero  at  oo,  integrand  must  vanish  at  infinity,  so  this  behaviour  occurs 
by  means  of  the  quotient  of  two  large  numbers  vanishing.  This  can,  and  does,  in  Example 
1,  lead  to  intrinsic  numerical  problems,  as  we  discuss  next. 

4.  For  large  t,  the  integrand  essentially  looks  like  ■  Assuming  /  has  a  zero  at 

infinity,  the  numerator  becomes  unbounded  (and  negative)  for  large  t,  and  the  integrand  is 
a  quotient  of  two  large  numbers.  Our  ability  to  integrate  this  numerically  depends  upon  a 
critical  balance  between  two  effects:  on  the  one  hand  the  quotient  is  getting  small,  while  on 
the  other  hand  the  denominator  is  getting  large.  If  the  rate  of  increase  of  the  numerator  is 
sufficiently  less  than  that  of  the  denominator,  then  we  can  start  neglecting  the  integrand  for 
numerical  purposes  (setting  it  to  be  zero  overall)  before  the  size  of  the  denominator  exceeds 
the  value  allowed  by  machine  wordlength.  In  Example  1,  we  see  a  case  in  which  this  balance 
is  at  a  marginally  acceptable  level.  In  the  other  examples,  there  is  no  problem. 
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6.4  Examples 

Example  1.  Our  first  example  is  a  transfer  function  model  for  the  damped  Euler- Bernoulli 
model  presented  in  Section  2.1.  We  performed  the  computations  for  this  model,  and  com¬ 
pared  the  accuracy  against  the  phase  of  the  inner  factor  computed  by  forming  the  Blaschke 
product  from  the  first  100  right-half-plane  zeros.  (By  adding  zeros,  it  is  easy  to  check  the 
latter  estimate  of  the  phase  has  essentially  no  error.)  In  Figure  6.1  we  plot  the  error  in 
computed  phase  of  the  inner  factor  as  a  function  of  frequency. 


Frequency  in  radians 

Figure  6.1:  Error  in  inner  factor  phase  for  beam. 

In  order  to  attempt  to  reduce  the  error  in  the  quadrature  procedure  it  is  necessary  to 
examine  the  source  of  the  error.  The  infinite  integral  is  approximated  using  a  finite  interval, 
so  there  is  an  error  contribution  from  the  “tails”  of  the  integral.  Since  the  integrand  is 
smooth,  the  adaptive  quadrature  routine  easily  computes  the  integral  over  the  finite  range, 
so  long  as  the  integrand  can  be  evaluated.  Using  the  asymptotic  behavior  (2.23)  above,  the 
tail  looks  like 


A  arg(/0(iu;)) 


dt  (6.3) 


(6.4) 


which  explains  the  linear  error  behavior  seen  in  Figure  6.1.  From  the  expression  (6.4),  it 
is  clear  (as  it  is  intuitively)  that  we  should  lengthen  the  interval  of  integration  in  order  to 
reduce  the  error. 

Unfortunately,  the  interval  [ — 10*,  10s]  used  in  this  computation  is  roughly  the  largest 
interval  we  can  use  with  the  double  precision  computer  program  employed  here,  without 
replacing  the  integrand  by  an  approximation.  The  reason  is  that  the  smallest  non-zero 
number  representable  on  the  computer  used  is  roughly  2  x  lO"30*.  Using  the  approximation 
(2.23),  we  see  that  at  roughly  u  =  3  x  107  the  magnitude  of  the  transfer  function  saturates 
machine  precision.  In  practice  the  problem  manifests  itself  not  as  an  underflow  but  as  an 
error  in  the  compiler’s  logarithm  computation,  when  the  program  is  unable  to  normalize  a 
floating  point  number  due  to  a  saturated  exponent. 
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Actually,  the  naive  approach  to  computing  the  integrand  has  problems  with  an  even 
smaller  quadrature  interval,  since,  in  the  formula  (6.1)  for  the  transfer  function,  the  denom¬ 
inator  is  roughly  the  square  of  reciprocal  of  the  transfer  function  at  high  frequency.  It  is 
only  by  carefull  decomposition  and  normalization  of  the  integrand  that  we  are  able  to  attain 
the  full  interval.  We  have  not  explored  the  use  of  the  asymptotic  behavior  (2.23)  as  a  direct 
approximation  in  order  to  extend  the  interval  of  integration. 

Example  2.  Consider  the  case  of  a  pure  time  delay  cascaded  with  a  rational  outer  factor 
which  is  non-zero  on  the  imaginary  axis. 


T(s) 


3  +  1 


In  this  case  we  can  directly  find  the  phase  of  the  inner  factor,  and  in  fact  the  quadrature 
technique  can  also  attain  an  essentially  perfect  result. 

Example  3.  We  consider  the  transfer  function  of  Section  2.2  of  a  flexible  beam  having 
a  tip  mass  at  one  end  and  a  motor  at  the  other,  taking  into  account  both  bending  and 
torsion.  Inspection  of  equations  (2.45  -  2.47)  quickly  leads  one  to  conclude  that  an  explicit 
inner-outer  factorization  will  not  be  possible.  We  have  not  yet  performed  computations  for 
this  model,  but  it  seems  clear  that  there  is  little  hope  checking  the  results  by  taking  a  finite 
Blacshke  product. 


6.5  Conclusions 

We  have  presented  a  numerical  technique  for  the  computation  of  the  inner  factor  of  a  stable 
transfer  function,  intended  for  use  in  solving  H°°  problems  when  an  explicit  factorization  is 
not  available.  Preliminary  computational  results  show  that  the  technique  works,  although 
it  is  also  clear  that  when  the  zeros  of  an  infinite  Blaschke  product  are  readily  computable, 
this  is  not  the  method  of  choice.  Work  on  applying  this  method  to  the  beam  in  Example 
3,  among  other  models,  is  in  progress,  along  with  the  integration  of  this  computation  into  a 
solution  for  the  overall  problem  addressed  by  (6.2). 
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Chapter  7 
Conclusions 


We  believe  that  our  work  to  date  has  shown  that  there  is  a  reasonable  prospect  that  direct 
H°°  design  of  feedback  controls  for  distributed  parameter  systems  will  develop  into  a  practical 
procedure.  The  work  reported  here  shows  some  of  the  potential,  both  for  computing  actual 
designs  and  for  finding  the  limits  of  attainable  performance.  We  have  also  been  able  to 
obtain  some  insight  into  the  nature  of  the  optimality  criteria  considered  by  examining  exact 
solutions.  Obviously  much  remains  to  be  done,  in  terms  of  increasing  generality  of  the  theory 
(to  cover,  for  example,  more  complex  SISO  systems,  MIMO  systems,  and  more  realistic 
and  comprehensive  design  criteria)  and  in  terms  of  developing  numerical  procedures  for 
computing  solutions. 
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Appendix  A 

Some  Projection  Formulas 


Here  we  shall  prove  some  projection  formulas  used  in  Chapter  4.  These  projection  formulas 
show  how  to  compute  the  projections  of  11+  and  II_  for  some  classes  of  L2  functions  by 
means  of  residue  theorem. 

We  shall  use  C  to  denote  the  complex  plane,  C+  and  C~  to  denote  the  open  right  and  left 
half  planes.  H3,  H3  ,  H°°  and  are  Hardy  spaces  on  the  half  plane. 

Lemma  10  For  a  rational  function  Q(s)  on  C  with  Q(ju> )  €  L°°(jR), 

|£J(x  +  jui)  —  Q(jw) |  — ►  0  uniformly  in  u> 


as  x  — ♦  0. 

Proof:  Since  Q  is  proper,  lim.-.a,  Q(s)  =  q  exists.  So  for  any  t  >  0,  there  exists  an  L  >  0 
such  that  for  js|  >  L, 

\Q(3)  <  \ 

So  for  |u/|  >  L,  we  have 

\Q(x  +  j(j)  -  Q{ju>) |  <  | Q(x  +  ;w)  -  g|  +  \Q(ju>)  -  g|  <  |  |  <  e. 

Since  Q  £  L°°(jR),Q  has  no  poles  on  j'R,  so  it  is  continuous  on  j R.  For  sufficiently  small 
V  >  0,  Q(s )  is  continuous  on  [ — 77, 77]  x  [—jL,jL\,  so  it  is  uniformly  continuous  on  this  closed 
rectangle.  Therefore  there  exists  a  6  >  0  such  that  for  |x|  =  |(x  4-  ju>)  —  ju>|  <6  we  have 

I Q(x  +  ju)  -  Q(ju) I  <  £. 


I 


Lemma  11  For  Q  as  in  Lemma  10  and  g  6  Hi, 

lim  f  | (Qg)(x  +  ju)  -  ( Qg)(juj)\du  =  0. 

*—0-  J-oo 


Proof:  Note 

/  l(<?0)(*  +  i")  -  (Qg)(ju)\<lu  <  f  \Q(x  +  ju)\\g(x  +  ju)  -  g(ju)\du+ 

J  —  OO  J  —  oo 


+ 


J. 


\g(]Uj)\\Q(x  +  ju)  -  Q(ju)\du>  <  \\Q\\l~  [  \g(x  +  ju)  -  g(jw)\du+ 

J  —  OO 


+ 


sup \Q{x  +  ju >)  -  Q{ju) |  /  \g(juj)\<L)  =  {1}  +  {2}. 

U l  J  —  oo 


The  first  term  tends  to  0  by  a  result  for  Hp  spaces  [26,  p.  128].  The  second  term  tends 
to  0  by  Lemma  10.  I 

Lemma  12  For  f(jui)  E  L2; 

J_  [°°  ,  -  J  (n +f)(s),3  6  C+; 

27T  7-00  jw-3  '  l  (n_/)(s),  5  €  C- . 

Proof:  See  [7,  p.  195] .  I 

Formula  Al:  For  a  function  /  =  Qh -  with  Q  E  RH°°  and  6  Hi, 

=  for  JSC+- 

where  7  is  a  contour  in  C"  and  encircles  the  poles  of  Q. 

Proof:  From  Lemma  12  we  only  need  to  show  §  ^-dz  =  j  /"  ^^dw. 

Let  A  be  the  region  encircled  by  7.  Since  is  analytic  on  C-  \  A  for  3  E  C+,  we  see 
for  any  R  >  0  sufficiently  large  and  5  >  0  sufficiently  small, 

jl 

j-r  3  —  z  J rRii  3  —  Z  J-R  3  —  (  —  0  +  JU>)  Jc,  3  —  He1” 

where  |  <  a  <  p  < 

For  R  large  enough  we  have  \^t\  <  5-  We  know  h^{Re^9)  — *  0  uniformly  as  72  — ►  00  for 
§  <  9  <  [26,  p.125].  So 


|J  J“  1  S  IWIIl-  J"  <  2||Q||l-  J°  - 0 

as  R  — ►  00.  So  we  have 

/  iw  * = ,  r  JtLtM-i,. 

Jy  3  —  Z  J- 00  3  —  [  —  0  +  ]<JJ) 

Also  notice  that  for  €  H2,  h-(ju>)  E  L2,  so  ^  L1.  By  Lemma  11, 


lim  r 


/(-5  +  ;u/) 

3  - (-6  +  jw) 


djuj 


-j: 


Md„. 

00  3  —  JCJ 
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Therefore 


Jy  3  —  Z  J- 


*  £M.t 

oo  3  —  ]L U 


This  proves  Formula  (Al). 


By  Formula  Al  and  the  residue  theorem,  we  have 
Formula  Al':  If  Q  has  poles  Pi,P2,--,Pn  in  C-,  then 


(n+/)W  =  h  /  ~iz  =  t 

ZTTJ  J-y  3  Z 


A  RQ(pt)h-(pt ) 


s  -  Pi 


where  RQ{pi)  denotes  the  residue  of  Q  at  p{. 

Similarly,  we  have  corresponding  formulas  for  a  function  f  =  Qh+  with  Q  £  /ZH“  and 
h+  £  H2.  We  state  them  explicitly  as  follows: 

Formula  A2:  For  a  function  /  =  Qhi  with  Q  £  and  k+  £  H2, 

rn  c\f  -\  _  ^  /  f(Z)  r _  - 


(II_/)(s)  =  —  — —dz  for  s  £  C  , 

J 7  3  —  Z 


where  7  is  a  contour  in  C+  and  encircles  the  poles  of  Q. 
Formula  A2':  If  Q  has  poles  pi,P2,-,Pn  in  C+,  then 


(n-/)M  = 


2irj  J-y  3  —  z 


”  R9(p.)h-(P.) 


3  ~  Pi 


where  RQ(pl)  denotes  the  residue  of  Q  at  p^. 


Appendix  B 

Computation  of  the  Eigenvalues  and 
Eigenfunctions  of  T 


In  this  appendix,  we  shall  use  the  projection  formulas  developed  in  Appendix  A  to  provide 
explicit  formulas  for  the  computation  of  the  eigenvalues  and  eigenfunctions  of  T. 

Recall  that  for  an  eigenvalue  /x2  and  an  associated  eigenfunction  we  have 


(/x2  -  W*Wx)x^  = 


(B.l) 


and 


=  n_A*(n+A  -  n_£)  -  a*u+a  -  n „w;wx  +  atls  +  n+s*[n_(A  +  #)]  (b.2) 

where 

W'W, 

_  wl  wx 

W  1  ” 

B  =  -WN0XD-, 

ww  =  w*wx  +  w;w2. 

Further  we  explicitly  write  out  Wx  and  Dx  as  follows 

F(s) 


Wi  = 


(s+Wi)---(s+  Wk) 
where  F(s )  is  a  polynomial  in  s,  G  C+,ji  =  1, fc,  and 

(dx  -  s)  - ■  ■  (dn  -  3 ) 


e  H°°, 


Dx  = 


(d\  +  s)  ■  •  •  (dn  +  3 ) 


where  d,  (E  C+ ,  j  =  1 , . . . ,  n. 


Remark  27  1.  By  the  definition  of  A,  it  is  easy  to  see  that  the  left  half  plane  poles  of  A 

are  {-wu  -u>fc}  =  V{WX). 
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2.  Since  Bw  is  a  Blaschke  product  such  that 


W* 

Bwt+  €  H°° 


W 


and  B ^  is  a  Blaschke  product  such  that 

B^jj.2 -F'F)~l  €  H°°, 

order(Svv)  =  number  of  poles  of  ^  on  the  right  half  plane  =  order (B^).  We  see  that 

order(Bw)  =  order(5;i). 

3.  By  the  definition  of  Bw,  we  know  that  the  right  half  plane  poles  of  jpj-  are  {—  pfw  : 
J  =  1,  ...,771},  or  equivalently,  the  left  half  plane  poles  of  are  {pfw  :  j  =  1, 

There  are  five  terms  in  (B.2).  By  Proposition  2,  Rx»(pfw)  =  0 ,j  =  In  the 

generic  case,  we  write  down  the  basis  functions  for  the  range  of  each  term  as  an  operator  on 
as  follows: 

1st  term: 


2nd  term:  (i) 

(n) 

3rd  term: 

4th  term: 

5th  term: 


s  —  wi  ’  ’  s  —  in* 


A*— A*  ‘ 


Am- 


3  +Wk 


=  1,2,..., 77i. 


3  +  WX 

1 

si 

3 -Pi 

1 


3  —  W\  '  '  3  —  Wk 


Am—Lr,...,A0  1 


s  —  di 

1 


3  -  dn 

1 


3  +  d\  3  +  dn 

We  shall  use  —  1, ...,  2k  +  m  -I-  2tx  to  denote  the  above  2k  +  m  +  2n  basis  functions 


of  A. 


So  we  see  that  the  rank  of  Au  is 


rank(  AM)  =  2k  4-  m  +  2n  =  2[order(  Wi)  +  order(Z?,)j  +  order  (B^). 
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Remark  28  In  [50,  Lemma  3],  the  authors  consider  the  stable  plant  case  and  state  that 
rank(A^)  —  2[order(G)  +  order(F)]  in  the  generic  case.  In  the  stable  plant  case  G  reduces 
to  - .  In  fact  in  [50]  G  is  also  assumed  to  be  in  H°°.  Therefore,  in  that  context  BwG 
is  substituted  for  G,  and  B\yM  for  M .  Of  course,  BwG  has  the  same  order  a3  G ,  which 
is  order (W~x )  =  order  (F).  By  he  definition  of  B Z?M(/i  —  FmF)_1  G  H°°,  so  order (F)  = 
order(i?M).  But  we  show  that  rank(A^)  =  2order(H^i)  +  order(B^).  (In  the  stable  plant  case 
Dt  =  1  so  order (Dt)  =  0 .)  Since  order (W~l)  =order(Wi^  +  order (W^)  in  the  generic  case, 
we  conclude  that  the  generic  order  of  the  mixed  sensitivity  problem  is  lower  than  that  m  the 
abstract  problem  considered  in  [50]  by  2  orde^W^)- 

Let 


+ 1  >  ^2 k)  ^2fc  +  l  j  •••)  Clk  +  mi 

c2  fc  +  m+2 
Clk  +  m  +  n  +  j 


From  (B.l)  we  have 


[(n+A-n_5)iJK) 

(II+i4xJ1)(u/i)  -  (U..Bx„)(wj),j  =  1, k, 


E(d])RD‘  (d])x^(dj),j  =  (B.3) 

E-Rd'[II4A  +  B)xtl]{-d,),j  =  1, ....  n.  (B.4) 


(p2  —  W]Wi  )xM  =  AMxM 

2A;+TTi-f-2n 

=  S  CA-  (B-5) 

2  =  1 


In  the  following,  we  shall  show  how  to  get  a  system  of  2k  -f  m  +  2n  homogeneous  linear 
equations  in  cx , c2fc4.n+jn. 

1.  Solve 

p2  -  w;wx  =  0 

to  get  2k  roots  r1,...r2fc.  So 


2k+rr.  +  2n 

5Z  c2^j(ri)  =  0,1  =  l,...,2fc. 

2  =  1 


2.  By  Proposition  2,  Rx>1{pfw)  =  0,  /  =  1, ...,  m.  So  we  have 

0  =  lim  (pfw  -  s)(p2  -  W(Wi)x^ 


lim  (pfv  -  s)AMxM,  for  l  — 

~*P{ 


which  are 

k 


T 

7=1  (“i  +  pT 


E; 


c2fc-t-., 


e7 


^2fc+m+ 1 


1  \Pj  ~  Pl  )  2  =  1  Vd2  Pi 


=  0  for  /  =  1, 


,  m. 
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3. 


2i-fHi+2n 

[(/i2  -  W'W^x^di)  =  c]<j>J(dl)>  for  /  =  1, ...,  n. 

i=i 


[(/x2  -  W-W^K-di)  =  lim  (d,  +  s)(M2  -  W*Wx)x, 

=  lim  (di  +  j)(Amxm) 

#— ►  —  di 


=  (  lim  (4 +  ,)*')( 

— *  V^s  «>  -  4  frpf’  +  d, 


-  Y1  T~7  |  +  c3fc+m*„*(,  for  I  -  1,  ,n 
,=  !  “i  +  d>  / 


From  (B.3)  we  get 


=  E(l)^(«i,)’  for  ,  =  1- 


Since  xM  6  B^BwNtDt H2  ,  we  see  that 


W7W,  tfx« 

Ax„(-d<)  =  (  for  /  =  1, ...,  n. 


(n.^x„)(-4)  =  (AxJi-d,)  -  (n+^x-d,) 


^1  v  C*+J  ™  c2k+, 


(n.Bx„)(-*)  =  -i;x~i,  for  1=1,  ...,n. 

;  =  1  +  “I 


From  (B.4)  we  have 


^  ^  ^  d')[c2fc+m+n+J  +  {E'R?'){  d')(E  ^  2d{  E  pBu  ^di 


•%—'  c2k+m+j 


)I.  ^  f=l, -,n. 


Summarizing  the  results  from  above,  we  have  the  following  procedure  to  find  eigenvalues: 
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1.  Solve 

M2  -  w;Wx  =  0 

to  get  2k  roots  ri,...,r2*. 

Let 

atj  =  4>j(rt),  for  i  =  1,  ...,2k,  j  =  1,  ...,2k  +  m  +  2n. 


2.  For  z  =  1, m,  let 


1 


a2fc+.,*+i  - 


w,  +p, 


for  j  =  l,...,k) 


1 


a2k+t,2k+]  —  B  B  i  f°r  J  —  1 ,  —  j  77T-; 

Pj  ~  P, 

f°r  J  = 

dj  —  px 


3.  For  z  =  1, n,  let 

a2*+m+.,j  =  <£;(£),  ^  j  =  1,  — >2A  +  m; 

/  ^fc+m+j  W  ~  (<*>),  for  t 

&2fc4-m+»,2fc+m+j  <  -  *  .... 

I  ^2fc+m+J(<i.),  for  = 

a2Jt+m+i,2Jlc+m+n+j  =  ^fc+m+n+j  ( )  j  for  J 

where  E  =  —WN0X. 

4.  For  i  =  1, n,  let 


^  2  Jk + m +n+ » +j 


m+n-M, 2Ac+j  — 


^2fc+m+n+*,2ik+m-f  i 


<^24s -*-m -f », 2* -4-m-fn+j  — 


(M»  -  a  ...  1 

{l  w-wJT-  R  K-*)-" 

for  j  = 

w;wxn:  n  ;  v  x)sdt  +  Pfu 

for  j  =  1, ...,  77i; 

{[  w;wxn;  ~  R  ]{~dx)  +  R  (_d’)} 4  +  4’ 

for  7  =  1, ... ,  7i; 

/  ~  for  1  =  hi  =  i.-.n; 

l  0,  for  7,7  =  1.  •■■»"■ 


5.  Let  all  other  atJ  =  0  for  1  <  z , 7  <  2*  +  m  +  2n. 
We  form  the  matrix 


84 


A{y.)  = 


/  an 

al,2*+m+3„ 

\ 

a2fc,l 

“2*,3*+m+3n 

0 

0 

a2*+l,*+l 

a3fc+l,2fc  +  m+n 

0 

r 

0 

0 

a1k+m,k+l 

a2fc-fm,2fc  +  m+n 

c 

0 

aj*fm+l,l 

a3*  +  m+l,3*  +  m+3Ti 

a2*+m+n,l 

a2fc-f  *n+w,2b-f  m+3n 

0 

0 

a3/k  +  r7»+n+l,*+l 

a2  k  +  m+ ni  ,2  k  +  m+  n 

a3*+m+n+l,3*+m+n+l 

0 

\  0 

0 

a3fc+m+ 3n,fc-f 1 

a2fc+m-f  2n,2fc+m+  n 

0 

a3fc+m+2n,2fc  +  m-f  2n 

/ 

and  we  can  then  write 

A{n)C  =  0,  (B.6) 

where  C  =  (clt ...,  c2t+m+2n)T. 

In  order  for  (B.6)  to  have  nontrivial  solutions,  we  set 

=  0,  (B.7) 

and  solve  (B.7)  for  /i2  to  be  am  eigenvalue. 
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Appendix  C 
Proof  of  Lemma  2 


Define 

X  =  W  -  MQ  . 

Let 

Xn(ju)  4  W(ju;)  -  M{juj)Tn{ju). 

We  now  consider  the  magnitude  squared  of  Xn(jLu), 

=  I W(jw)  -  M(ju)T„(,u)\'  =  \W(jw)  +  A.0'")^"1  - 

using  the  definitions  in  (3.20).  Noting  the  assumption  (3.13)  we  take 

Woo  =  sup  Rw{°o), 

and  so  for  any  v  >  0  we  can  take  n(v)  large  enough  so  that 

\W(ju>)\  <  Woo  +  v  a.e.  for  |u/|  >  n( u).  (C.l) 

Now  we  fix  u,  assume  that  n  satisfies  (C.l),  and  later  determine  how  to  pick  v  as  a 
function  of  e  to  satisfy 

\\Xn\\<in  +  e.  (C.2) 

Suppose  u>„  is  a  frequency  (possibly  oo)  at  which  ||Xn||  €  Rxn(vn)-  If  a>n  is  finite,  define 
V/n  =  sup  Rw(wn),  h  =  IMi^n)!,  S  =  arg[/in(;o>n)]  and  a  =  a(un). 

Otherwise  we  know  that  ||An||  =  /x,,  and  we  define 

Wn  =  nx,  h  =  0,  6  =  0  and  a  =  0. 

Note  that  these  are  functions  of  n,  as  is  u)n.  We  also  note  for  later  use  that  6  satisfies 
0  <  6  <  2ir/n. 

Now  we  show  that  the  norm  of  Xn  approaches  the  infimal  value  /x,  as  n  increases.  For 
e  >  0  let  Bc(ujn)  =  (u»n  -  e,ujn  +  e).  Then  for  all  u> 


86 


I^n(i^)|2  <  lim  ess  sup  \W{jQ  +  hn(jQ(X(jC)  -  W(jO)|2 

e^°  c  e  *.(<*.) 

=  lim  ess  sup  (|W(iC)|2  +  2Re  {  W”(»<)A»(*0(*(*0  -  ^(*0)}  + 

(  €  N'(vrx) 

\hn(i()(x(iO  -  W(i0)\2 ) 

=  lim  ess  sup  (\W(i()\2  -  2Re  (An(zC)} \W(i(,)\2+ 

c~*°  (ew.W 

2Jfe{w(*)M«0*(*0}  +  IMOI’l-fW)  -  ^(*01’ ) 

=  lim  ess  sup  (|W(iC)l*  -  2Be  <M*0}  |H^«)|2  +  l*»(iOIJ|W(>0|J 

+  M>OI’l*«)|a  +  2B«  {»"«)*»(*<)■*(*<)}  - 

2|/in(iO|Jfie  {X(i()W(i«}  ) . 

For  n  >  4,  0  <  a.ig(hn(i( ))  <  tt/2,  and  we  write 

l*-(j«)|2  <  lim  ess  sup  (|^K)|3  -  2|An(»C)|  ■  \W{iQ\2  +  M*C)lW>C)|2+ 
e~*°  (a,W  v 

2(\hn(i()\  -  Re  {^n(*0})|Wr(*C)l2  +  IM*  )I2I*(<C)I2  + 

2 Re  (M*C)}  ^  {W(i( )X(«C)}  -  2/m  {hn(i()}  Im  {W(»0*(*0}  - 

2|M<0r^{^C)wr*(iC)}) 

=  lim  ess  sup  (iW^tOI*  [(1  ~  IM*C)I)2  +  2(|M*C)I  ~  Re  {M*0»] 

e_*  C  e 

+IM*C)l1!-?(>C)l’  +  2  (Be  {»”(*<)*«)}  [Be  {*.(>()}  - 
IM*C)|J  ]  -  Im  {/>„(. C)}  Im  {W’(i <)*«)}  ) 

Using  continuity  of  hn  and  the  definition  of  we  get 

\Xn(ju>)\3  <  W^-  [(1  -fi)a  +  2/i(l  -coatf)]  +  +  2Wnnt{h(\  -h)  +  hsin6} 

<  [Wn{  1  -  h)  +  +  2 hS2Wl  +  2 WnfuhS 

Given  n,  there  are  two  possibilities: 

case  (i).  (cj„  >  n)  In  this  case  we  shall  use  the  fact  that  Wn  — >  W, x  as  wn  -*  oc. 
|*n(;u;)|2  <  {{W00  +  v)(l-h)  +  hnt}2  +  2h62{W00tv)2  +  2h>it8(Woo  +  i') 
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<  [(m  +  1/)(1  -  h)  +  hm}3  +  2 hS2(m  +  v)3  +  2 hfj,i6(m  +  u) 

—  [m*  +  v]2  +  2  +  v)  +  263(h*  +  i/)2 

~  Hi  +  v  ■  2/Xi(l  +  6)  +  2 fi2S  +  v3  +  2 63(ni  +  ^)2 

<  /i2  +  v  •  2/Xi(l  +  — )  +  47T/x?/n  +  v3  +  2  (— )  (fi,  +  v)2. 

ti  In/ 

Now  we  simply  state  how  to  pick  ne  and  t/(e)  such  that  n  >  ne  and  0  <  v  <  t/(e) 

guarantees  that  the  right  hand  side  of  the  last  inequality  gives 

\Xn{ju)\3  <  n2  +  e. 


Assuming  0  <  e  <  1,  we  take 

0  <  v(e)  <  — r — ,  and  ne  >  2r  such  that  \W(ju)\  <  \W H - —  for  |w|  >  ne, 

Id  Hi  16  Hx 

which  gives  v  ■  2/i,(l  +  2s )  <  e/4.  Taking 


(C.3) 


ne  >  - 

e 

gives  7 r/i2/n  <  e/4  for  n  >  nc.  Taking 

s  I 

gives  i/2  <  e/4.  Finally,  using  (22)  with  0  <  v  <  i/(e)  and  taking 


(C.4) 


(C.5) 


^  Sn(Hi  + 1) 

7lt  >  - 

e 


(C.6) 


allows  us  to  conclude  that  8  ( Hi  4-  v)3  <  e/4, 

Thus  (C.3)  —  (C.6)  together  will  ensure  ||^n||L  <  Hi  +  c  in  this  case, 
case  (ii).  ( ujn  <  n)  The  idea  in  this  case  is  that  h  — *  1  since  u/n  <  n.  Since  KWHo,,  >  i^n! 
and  h  <  1, 


\Xn(jujn)\3  <  [W^l-^  +  hH,]2 +  2hS3Wl)  +  2W00Hth6 

<  h-+(  1-h)-  2\\W\\00h>  +  II^IIL(1  -  h)3  +  4irW||W'||0o/» 

+2(^)Vj|i,  (C.7) 

We  proceed  to  satisfy  (C.2)  by  bounding  each  of  the  last  four  terms  of  the  right  hand 
side  of  (C.7)  to  be  individually  bounded  by  e/4.  Now  since  u/n  <  n  in  this  case, 


1  >  fcJ  =  IM;ey)l! 

=  W  +  uir 
>  [oV(7J  +  n’))1'" 
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It  is  easy  to  check  that  limfl-<00[a2/(72  +  n2)]1/"  =  1,  so  h2  ->  1  as  n  -*•  oo,  and  therefore 
(I  -  h)  — ►  0.  In  order  to  estimate  the  required  n€l  we  want  (1  —  /i)  •  2||W||/ij  <  e/4.  It  is 
sufficient  to  take 


log(1-ip^;)<slog[“1/(7l+n,)I' 


For  n  >  7 


a2/(72  +  n2)>a2/2n2, 


so  a  sufficient  condition  is 


W1  _  £  n  .  21og(q)  -  2  log(n)  -  log(2) 

ZWWWfiJ  2  n 

Assuming  n  >  max(2,  1/a2),  a  sufficient  condition  is  then 

Similarly,  ||W||^0(1  —  h )2  <  e/4  can  be  assured  by  requiring 

,  V*  s  21og(n) 

105(1  2||IV||  <  ~  n  ' 

47T/Z, ||  W\\oo/n  <  e/4  will  be  obtained  by  assuming 

„  ^  l^ia\\W\\ 


(C.8) 


(C.9) 


and2(^)2||^||L<6/4  by 


«.  >**\\wy-  ■ 


Thus  using  (C.3)-(C.10)  we  pick  nt  to  satisfy: 


ne  >  max  ^2ir,  !=£,  J,,  8,||^||^) 

\W(ju)\  <  I  Woo  I  +  g£-  for  |w|  >  ne 

>  max  (log(l  -  sp^R),  log(l  -  jjjgjj)) 


to  obtain  ||A„||2  —  /j,2  <  e.  | 


(C.  10) 
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